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ABSTRACT 

A  frequency  dependent  hybrid  mode  analysis  of  general 
coplanar  transmission  lines  is  presented.   A  Fourier  trans- 
form technique  is  applied  and  the  resulting  expressions  are 
evaluated  numerically  using  the  method  of  moments.   Shielded 
slot  line  resonators  as  well  as  dispersion  characteristic  and 
characteristic  impedance  of  open  and  closed  boundary  slot 
lines,  coupled  slot  lines  and  coplanar  strip  lines  are  obtained 
An  extension  of  this  technique  is  used  in  a  perturbational 
analysis  of  slot  lines  on  a  ferrite  substrate. 

The  theoretical  performances  of  various  slot  line  direc- 
tional couplers  are  predicted.   Numerical  results  are  compared 
with  experiments  in  various  cases,  where  excellent  agreement 
has  been  obtained. 
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1 . 0   INTRODUCTION 

The  study  of  electromagnetic  energy  transmission  is 
but  one  important  area  in  microwave  engineering,  where  the 
electromagnetic  waves  are  travelling  through  some  trans- 
mission medium,  which  provides  the  link  between  the  trans- 
mitting and  the  receiving  part  of  a  transmission  system. 
The  transmission  medium  may  be  the  free  space  or  it  may  be 
a  guidance  structure  which  will  direct  the  energy  flow  into 
a  certain  prescribed  direction.   This  structure  will  be  in 
general  a  boundary  where  the  surrounding  fields  interface 
with  the  surface  currents  and  charges  on  this  boundary  which 
will  be  in  most  cases  an  electric  conductor,  although  in 
some  cases  a  dielectric  material  may  be  used  instead.   Some 
standard  examples  for  guided  wave  structures  for  microwave 
applications  are  the  open  wire  transmission  lines,  coaxial 
cables  and  rectangular  or  circular  waveguides. 

1.1  TRANSMISSION  LINES  FOR  INTEGRATED  CIRCUIT  APPLICATIONS 

In  most  cases  of  designing  a  transmission  link  the  deci- 
sion of  which  guidance  structure  should  be  used  will  be  based 
upon  various  determining  factors  which  are  sought  to  be  opti- 
mal in  some  sense.   Some  examples  of  these  factors  are  the 
cost  of  design  and  maintenance,  the  transmitted  power  level 
and  the  frequency  band  of  operation,  to  name  only  a  few. 
While  most  of  these  examples  given  so  far  apply  mainly  to 
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large  scale  transmission  systems,  recent  advances  and 
increasing  use  of  integrated  circuit  technology  at  microwave 
frequencies  set  different  standards  where  new  transmission 
lines  are  required  that  are  planar  in  nature.   Several  rea- 
sons may  have  led  to  this  ever  increased  interest  in  planar 
transmission  lines;  e.g.  the  use  of  microwave  circuits  in 
satellite  communications  requires  light  weight  and  small 
sizes;  another  reason  is  that  the  etching  of  a  planar  trans- 
mission line  on  a  dielectric  or  ferrite  substrate  is  just 
another  production  step  within  the  fabrication  of  a  wholly 
automatically  fabricated  integrated  circuit  module.   The  most 
commonly  used  transmission  line  to  date  in  microwave  IC's 
is  the  microstrip  line  shown  in  Figure  1.1/1.   It  was  first 
proposed  in  1955  by  M.  Arditi  [Ref .  1]  and  has  since  then 
received  considerable  attention  by  innumerable  investigators. 

Although  the  investigations  of  single  and  coupled  micro- 
strip  lines  continue  to  be  the  subject  of  more  sophisticated 
research  today  only  recently  the  interest  shifted  towards 
investigations  of  coplanar  transmission  lines.   Some  examples 
of  these  lines  are  shown  in  Figures  1.1/2-4.   These  lines 
offer  the  advantage  of  being  more  easily  adaptable  to  shunt 
element  connections  without  the  need  to  penetrate  the  die- 
lectric substrate  as  in  the  case  of  microstrip  lines.   Active 
devices  as  well  as  passive  elements  may  be  fabricated  in 
either  series  or  shunt  configuration  directly  on  to  these 
lines  in  either  hybrid  or  monolithic  fabrication  techniques 
[Ref.  2]. 
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Figure    1.1/1      Microstrip  Transmission  Line 
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Fiaure   1.1/2      Slot  Transmission  Line 
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Figure    1.1/3      Parallel   Coupled  Slot   Transmission   Line 
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Figure  1.1/4   Parallel  Coplanar  Strip  Transmission  Line 
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The  slot  line  in  Figure  1.1/2  consists  of  a  narrow  gap 
between  two  conducting  surfaces  which  are  on  one  side  of  a 
dielectric  substrate.   A  fairly  high  permittivity  of  this 
substrate  confines  the  waves  travelling  on  this  structure  to 
the  region  of  this  slot  such  that  the  main  mode  of  propaga- 
tion is  similar  to  a  TE..  -.  rectangular  waveguide  mode  [Ref.  3]. 

Figure  1.1/3  shows  two  adjacent  slot  lines  parallel  to 
each  other  and  of  equal  width,  where  the  coupling  effect 
between  these  lines  may  be  advantageous  in  the  design  of 
microwave  filters  and  couplers.   If  this  configuration  is 
used  as  a  single  transmission  line  such  that  both  outer  con- 
ductors are  the  ground  electrodes,  this  line  is  called  a 
coplanar  waveguide  (CPW)  according  to  Wen  [Ref.  4], 

Both  lines  use  a  fairly  broad  area  of  conductor  surface 
which  may  serve  as  the  ground  plane  for  a  micros trip  placed 
onto  the  opposite  side  of  the  substrate,  thus  coupling  from 
one  transmission  line  to  the  other  can  be  easily  achieved 
[Ref.  5] . 

Another  promising  property  of  these  two  lines  is  the  fact 
that  both  provide  an  area  where  the  RF  magnetic  field  is 
elliptically  polarized  which  makes  these  lines  suitable  for 
nonreciprocal  applications  when  the  dielectric  is  replaced  by 
a  ferrite  substrate.   This  elliptic  polarization  is  in  con- 
trast to  microstrip  lines  where  a  similar  effect  can  be  achieved 
only  by  using  meander  type  lines  [Ref.  2].   The  coplanar  strip 
line  shown  in  Fig.  1.1/4  consists  of  two  strip  conductors  on 
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a  dielectric  substrate  with  no  ground  plane  and  can  be 
thought  of  as  a  dual  to  the  coplanar  waveguide  [Ref.  4], 

1.2  MATHEMATICAL  ANALYSES  OF  COPLANAR  TRANSMISSION  LINES 

In  the  design  of  a  transmission  line  the  two  important 
line  parameters  are  the  characteristic  impedance  and  the 
propagation  or  phase  velocity. 

A  standard  approach  in  the  analysis  of  these  lines  is  to 
consider  the  line  as  having  distributed  parameters  of  capaci- 
tance ar.d  inductance  and  to  solve  for  these  in  terms  of  the 
physical  dimensions.   This  solution  assumes  a  TEM  mode  of 
propagation  since  static  fields  are  assumed  for  the  calcula- 
tion of  the  line  reactances .   However  it  can  be  easily  proved 
that  the  inhomogeneous  guidance  structures  shown  in  Figures 
1.1/1-4  strictly  cannot  support  a  TEM  mode.   Assume  waves 
are  propagating  along  a  non  homogeneous  guiding  system  such 
as  a  microstrip  or  a  slot  line  with  a  propagating  term 
e  •*  where  co  is  the  radian  frequency,  t  the  time  and 

Y  =  a  +  j6  the  propagation  constant  assumed  from  now  on  to 
be  pure  imaginary  (lossless  case,  a=0) ,  j  =  v-1  . 

The  transverse  field  components  can  be  related  to  the 
longitudinal  field  components  through  an  algebraic  manipula- 
tion of  Maxwell's  curl  equations  to  yield 


,  3E  .      3H  . 

xi    k  2  „2  VJ  &i  3y      Jwdx 
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.  3E  .       3H  . 

Hyi  "  ^2  <^i  3^  +  *»  IT"' 

,  3E  .         3H  . 

Exi  =  ~T^  ^  33T1  +  **H  IT* 
ki  "B 

,  3E  .          3H  . 

Eyi  "  j^2  MB  37^  +  3»"i  3^ 

2     2 

where  k.   =  w  e.y.  i  =  1,2,3,...   defining  the  spatial 

region. 


If  the  TEM  mode  is  assumed  E   and  H   are  non  existing 

z       z  3 

throughout  the  guidance  structure.   A  solution  other  than 

2   2 

the  trivial  solution  clearly  can  only  exist  if  k.  -p   =0. 

Since  for  the  inhomogeneous  guidance  system  3  is  the  same 

for  all  regions  but  for  certain  i^j ,  k.^k.    the  TEM  mode 

i   J 

condition  cannot  be  met. 

But  although  a  planar  transmission  line  such  as  the 
microstrip  line  has  an  inhomogeneous  cross-section  a  TEM 
or  often  called  "quasi-static"  mode  approximation  can  be  used 
since  most  of  the  fields  are  confined  to  the  dielectric  be- 
tween the  two  conductors .   Recent  investigations  show  however 
that  this  is  only  valid  for  frequencies  up  to  approximately 
5  GHz  [Ref .  6] .   For  coplanar  transmission  lines  a  quasi- 
static  approximation  seems  no  longer  suitable  for  microwave 
frequencies  in  the  GHz  region  since  the  waves  are  not  as 
closely  bound  in  the  dielectric.   Investigations  of  coplanar 
transmission  lines  thus  should  start  from  the  general  field 
theory  to  include  the  hybrid  mode  which  provides  a  frequency 
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dependent  analysis.   This  means  that  in  general  a  solution 
to  Helmholtz's  equation  rather  than  Poisson's  equation  is 
to  be  found  with  inhomogeneous  boundary  conditions  due  to 
the  interaction  between  fields  and  guiding  structure.   This 
is  not  an  easy  task  in  a  mathematical  sense.   A  Green's 
function  solution  e.g.  should  be  feasible  at  least  theore- 
tically but  has  not  been  reported  in  the  literature  so  far, 
probably  due  to  the  immense  complexity  of  the  problem. 

In  his  analysis  of  the  slot  line  Cohn  [Ref.  3]  introduced 
electric  and  magentic  boundary  walls  which  permitted  the 
slot  line  configuration  to  be  treated  as  a  rectangular  wave- 
guide problem.   The  set  of  rectangular  waveguide  modes  was 
used  as  a  set  of  orthogonal  functions;  in  the  final  analysis 
the  fields  were  approximated  by  the  TE,Q  rectangular  waveguide 
mode. 

Wen  [Ref.  4]  investigated  the  coplanar  waveguide  and  the 
coplanar  strips  with  the  assumption  that  the  dielectric 
substrate  is  thick  enough  to  be  considered  infinite  for 
conformal  mapping  purposes.   For  large  values  of  the  permit- 
tivity this  may  be  a  practical  assumption  but  it  appears  to 
be  impractical  for  relatively  small  values  or  for  thin 
substrates.   An  alteration  of  this  method  used  by  Wen  is 
given  by  Davis  et  al .  [Ref.  7]  and  takes  the  finite  thickness 
of  the  substrate  into  account.   The  coplanar  strip  line, 
however,  is  not  treated.   Furthermore,  both  methods  lack  any 
frequency  dependent  information  on  the  behavior  of  the  line 
parameters  since  a  quasi-static  approximation  is  used. 
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Recently,  Itoh  and  Mittra  [Ref.  8]  proposed  a  new  method 
for  calculating  the  dispersion  characteristic  of  a  single 
slot  line.   Here  Helmholtz's  equations  for  the  scalar 
potential  functions  [Ref.  9]  are  transformed  into  the  Fourier 
domain  and  are  solved  for  in  accordance  with  the  given 
boundary  conditions .   The  method  of  moments  is  then  used  to 
solve  for  the  propagation  constant.   This  approach  showed 
good  agreement  with  results  from  Cohn  [Ref.  3] .   However  no 
analysis  of  the  characteristic  impedance  has  been  reported 
in  Reference  8. 

This  spectral  domain  approach  was  used  by  Luna  [Ref.  10] 
to  compute  the  dispersion  characteristic  of  coplanar  strips. 
No  further  analysis  of  this  transmission  line  has  appeared 
in  the  literature. 

1.3  UNSOLVED  COPLANAR  TRANSMISSION  LINE  PROBLEMS 

Besides  the  known  parameters  of  coplanar  transmission 
lines  given  in  Section  1.2  there  remain  several  problems  open 
for  analysis  which  will  be  shortly  listed  here  and  are  the 
subject  of  this  thesis.   Although  the  slot  line  was  studied 
in  19  6  8  by  Cohn  [Ref.  3]  with  a  frequency  dependent  analysis 
and  the  coplanar  waveguide  structure  was  analyzed  by  Wen 
and  Davis  et  al.  [Refs.  4,7]  by  a  quasi-static  method,  no 
complete  analysis  of  two  coupled  slot  lines  has  appeared  yet, 
although  originally  planned  by  Cohn.   Line  parameters  which 
describe  the  coupling  effect  between  these  lines  are  important 
in  the  design  of  directional  couplers  and  in  filter  design. 
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Moreover,  in  the  design  of  slot  line  filters,  the  application 
of  resonant  slot  structures  may  be  desirable  and  these  have 
not  yet  been  studied  analytically  at  all. 

Another  unsolved  problem  is  the  shielded  slot  line 
structure.   Since  waves  traveling  on  a  slot  line  are  not 
closely  bound  to  the  dielectric  next  to  the  slot  but  extend 
somewhat  into  the  air  regions  it  seems  important  to  study 
the  behavior  of  the  line  parameters  under  the  influence  of 
nearby  obstacles  such  as  an  enclosing  or  shielding  structure. 
These  results  could  be  of  some  significance  in  the  design  of 
microwave  IC  modules. 

A  further  investigation  is  needed  for  a  slot  line  on  a 
ferrite  substrate.   Although  recently  a  publication  by 
Kitasawa  et  al .  [Ref.  11]  discussed  this  problem,  the  analysis 
presented  is  restricted  since  the  operating  frequency  is 
assumed  to  be  far  above  the  resonance  frequency  of  the 
ferrite.   No  analysis  about  attenuation  of  the  waves  on  the 
ferrite  is  mentioned. 

Finally,  although  the  frequency  dependent  dispersion 
characteristics  of  coplanar  strips  were  analyzed  by  Luna 
[Ref.  10]  a  frequency  dependent  formulation  of  the  charac- 
teristic impedance  is  still  missing. 

1.4  CONTENTS  AND  SUMMARY 

Inhomogeneous  waveguide  structures  such  as  coplanar 
transmission  lines  will  not  support  a  TEM-mode  of  propagation 
at  high  frequencies.   Therefore  any  analysis  will  assume  a 
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full  or  hybrid  mode  which  may  be  thought  of  as  being  a 
superposition  of  TE  and  TM  modes.   This  formulation  will 
yield  a  frequency  dependent  phase  velocity  and  characteristic 
impedance . 

This  thesis  is  devoted  to  a  hybrid  mode  analysis  of 
coplanar  transmission  lines.   In  Section  2.0  open  boundary 
structures  are  analyzed.   Following  the  proposed  method  of 
Itoh  and  Mittra  [Ref.  8]  the  dispersion  characteristic  of 
a  single  slot  line  is  formulated.   The  investigations  are 
extended  by  comparing  various  possible  approximations.   Using 
Parseval's  relation  it  is  shown  that  the  power  flow  on  this 
waveguide  structure  can  be  computed  in  the  spectral  domain. 

The  resultant  calculations  for  the  characteristic 
impedance  are  compared  with  results  from  Reference  12 .   This 
is  the  first  verification  of  Conn's  method  by  an  entirely 
different  numerical  analysis  and  shows  very  close  agreement. 
The  theory  is  extended  for  an  analysis  of  two  coupled  slot 
lines  of  equal  width,  the  results  are  new  and  very  important 
for  the  microwave  industry  in  the  design  of  slot  line  couplers 
and. filters .   A  frequency  dependent  analysis  of  coplanar 
waveguide  is  included. 

An  easy  extension  of  the  theory  is  then  used  to  analyze 
the  line  parameters  of  a  coplanar  strip  transmission  line. 
The  final  part  of  Section  2.0  shows  a  perturbational  expression 
for  the  propagation  constant  of  slot  lines  on  a  ferrite 
substrate. 
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In  Section  3.0  various  closed  boundary  structures  are 
analyzed.   The  shielded  slot  transmission  line  is  compared 
with  the  open  boundary  structure.   Resonating  slot  configura- 
tions are  investigated  such  as  a  single  resonant  slot  and 
two  end  coupled  resonating  slots.   A  byproduct  of  the 
resonant  slot  investigation  is  the  end  effect  analysis  of 
slot  lines.   Theoretical  results  are  compared  with  experi- 
mental data  from  open  boundary  experiments  and  show  good 
agreement. 

Section  4.0  is  devoted  to  an  experimental  verification 
of  the  numerical  data  for  coupled  slot  line  parameters. 
Using  basic  network  theory  the  theoretical  response  of  a 
slot  line  directional  coupler  can  be  computed.   Due  to  the 
non-equal  phase  velocities  in  the  two  natural  modes  co-  or 
contra-directional  coupling  can  be  achieved.   Two  slot  line 
couplers,  a  co-  and  a  contra-directional,  were  built  and 
their  performance  measured. 

In  Section  5.0  remarks  are  made  about  the  numerical  analysis 
and  the  computer  programs.   It  is  believed  that  results  of 
this  thesis  are  of  interest  to  the  microwave  industry. 
Hence  two  of  the  most  important  computer  programs  used  in 
this  work  are  attached  to  the  thesis  and  their  use  explained 
in  this  section. 

Section  6.0  shortly  outlines  an  attempt  to  analyze  a 
circular  resonating  aperture  on  a  dielectric  substrate  and 
shows  the  difficulties  encountered  there.   Suggestions  for 
expansions  of  the  presented  work  conclude  this  section. 
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2.0   OPEN  BOUNDARY  STRUCTURES 

In  this  section  open  boundary  coplanar  transmission  lines 
are  analysed  where  it  is  assumed  that  the  lines  extend  infi- 
nitely in  the  z-direction  in  an  infinite  x-z  plane.   It  is 
further  assumed  that  conductors  are  negligibly  thick  and  that 
dielectric  and  conductor  are  perfect. 

The  electromagnetic  fields  surrounding  the  guidance  struc- 
ture are  derived  in  a  most  general  way  from  hybrid  modes  where 
use  is  made  of  the  fact  that  any  hybrid  mode  can  be  thought 
of  being  a  superposition  of  TE  and  TM  modes  by  the  orthogon- 
ality properties  of  open  boundary  waveguide  modes ,  proved  by 
Collins  [Ref .  13] .   The  derivation  of  all  equations  is  general 
and  not  limited  to  any  specific  mode  of  propagation.   Approxi- 
mations for  the  lowest  order  mode  are  made  in  the  very  end 
and  can  be  easily  modified  to  yield  results  for  higher  order 
modes . 

It  is  further  assumed  that  the  operating  frequency  is  kept 
below  the  frequency  for  the  onset  of  surface  wave  modes  which 
could  propagate  in  the  dielectric-air  medium  for  y  <  D  (see 
Figure  1.1/2-4) . 

This  onset  of  surface  wave  modes  for  the  slot  line  con- 
figuration e.g.  can  approximately  be  determined.   Consider  the 
region  y  <  D  and  assume  that  the  dielectric  has  the  permit- 
tivity e2  while  for  y  <  0  the  permittivity  is  that  of  free 
space  e  .   All  further  assumptions  stated  above  hold  and  the 
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existence  of  a  slot  in  the  metal  is  irrelevant.   Attwood 
in  Reference  14  has  derived  the  determinantal  equation  for 
TM  surface  wave  modes  which  yield  the  lowest  order  mode. 


This  equation  is  given  as   k  ~  tan (k  ~D)  =  -  —  k    where 
^ ^  c_2 *  c2        e    co 

k  .  =\/k02-B2   and  k  n   =  V62-k  2   ,   k  2  -  0)2e  P   and 
c2     2  co         o       o       o  o 

2     2 

k2   =  w  e_y   .   A  bound  for  the  lowest  order  cutoff  wave- 
length can  be  found  by  approximating  this  transcendental 
equation.   For  k  2  tan(k  2D)  <  0  ,   tan(k  2D)  <  0   or 
k  2D  >  j  .   But   k  2  =Vk2  "^        so  that  the  following 


inequality  holds   k2D  >\fk22-$2    D  >  J  .   If   A2  = 


OJ 


V^o' 

the  bound  for  the  onset  of  surface  waves  is 


X2  <  4-D 


2.1   SINGLE  SLOT-LINE 

The  schematic  of  a  single  slot  line  and  its  orientation 
with  respect  to  a  three-dimensional  rectangular  coordinate 
system  is  shown  in  Figure  1.1/2.   The  substrate  is  assumed 
to  be  dielectric.   Let  this  transmission  line  be  divided  into 
three  regions  in  space  and  define  these  regions  as 


gion    1 

for 

D   <  y 

2 

0    <   y    <_  D 

3 

y    <    0 
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Since  the  lossless  case  is  assumed,  the  propagation  constant 
Y  -  a  +  j$   is  pure  imaginary  with  a   =  0.   The  following 
derivation,  however,  is  perfectly  general  and  the  above 
simplification  is  made  only  at  the  end. 

2.1.1   Dispersion  Characteristic 

In  the  construction  of  solutions  to  Maxwell's  field 
equations  for  a  transmission  line  problem  it  is  sufficient 
to  assume  two  scalar  potential  functions,  a  scalar  electric 
potential  <j>  and  a  scalar  magnetic  potential  ^  which  may  be 
both  related  to  TM  and  TE  modes  respectively  [Refs.  9,13]. 
The  axial  components  of  TM  and  TE  modes  are  then 


E  -=  k  .  2cj)(x,y)e±YZ  (2-1) 


and     H   =  k  .  2ip  (x,y)  e±YZ  (2-2) 

£»  CI 


2      2     2  .  . 

where   k  .  =  k.   +  y  i      i=l,2,3   defining  the  spatial  region 

of  the  transmission  line  as  stated  in  Section  2.1. 

Through  Maxwell's  curl  equations  the  transverse  field 

components  are  then  determined  by  these  axial  components  and 

can  be  given  as 


EXi  =  ^  33T  -  J»»  *r>  eY  •      (2"3) 

Eyi  -  <Y  W-  +  JWP  33T}  eY"  (2"4) 
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Hxi  =  (Y  3iT  +  ^Ei  3^'  eY  (2-5> 

Hyi  =  ^  JT  '   JW£i  ^  e  (2"6) 

1=1/2,3 

where  propagation  in  the  z  direction  is  assumed  and  the 

dependence  of  the  scalar  potentials  on  the  variables  x  and 

y  is  understood. 

The  scalar  potential  functions  satisfy  Helmholtz's 

equations  in  the  three  spatial  regions ,  thus 


(vxy2  +  kci2>  *i  "  °  (2"7) 


and      (V   2  +  k  .2)  if..  =  0  (2-8) 

-Xj         CI       ■*■ 


where  V    denotes  the  two-dimensional  Laplacian  operator 

Xy  cr  tr 

in  the  transverse  (x,y)  direction.   So  equations  (2-7)  and 
(2-8)  are  second  order  partial  differential  equations  (P.D.E.). 
A  common  practice  in  solving  these  P.D.E.  which  essentially 
present  a  boundary  value  problem  is  to  transform  these  into 
ordinary  differential  equations  (O.D.E.)  by  an  integral 
transform  technique  [e.g.  Ref.  15]  which  may  then  be  solved 
analytically.   The  boundary  conditions  for  (2-7)  and  (2-8) 
are  intuitively  obvious,  for  instance   <+> .  (x,y)  ■*■  0   for 

00 

2 

x,y  ->  ±°°  and  furthermore    /  j(J>.  (x,y)  |  dx   is  bounded  for 

—  00 

finite  energy  reasons  —  the  same  holds  if  <f>  is  replaced  by  \p   — 
since  the  time  average  power  flow  is  related  to  the  scalar 
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potentials  as  shown  in  Reference  13  and  is  assumed  to  be 
finite.   Although  no  rigorous  proof  will  be  given  it  is 
therefore  assumed  that  the  scalar  potential  functions  can 
be  transformed  into  the  Fourier  domain  via  a  Fourier 
transformation  defined  as 

CO 

F{<|>(x,y)}  =  $U,y)  =  /   <f>  (x,y)  ejax  dx         (2-9) 

—  00 

If  this  transform  is  applied  to  equations  (2-7)  and  (2-8) 
the  following  O.D.E.  are  obtained: 


32  2 

*--■   <*>•  (a,y)  =  y.  4*.  (a,y)  (2-10) 

8y   x  1 

32    ,..  ,;. .,    2 


9y 


j   H'i(a,y)  =  Y±  ^(Oiy)  (2-11) 


.        2     2     2.2 
where  y.      =  a   -  y      -  k. 


Again  for  energy  reasons  certain  solutions  to  these 
O.D.E.  may  be  eliminated  by  boundary  conditions  in  the  three 
spatial  regions  since  for  y  -*■  °° ,   <j)(x,y)  ■*•  0   and  so   <Ha,y)  ->■  0 

The  solutions  to  (2-10)  and  (2-11)  may  then  be  given 


as 


-Y1(y-D) 

$-|_(a,y)  -   A  (a)e 

$2(a'Y)    =  Be(a)sinh  y?Y  +  C  (a)cosh  y ~y     (2-12) 

e     YTy 
$->  (a,y)  =  D  (a)  e 
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and  .     -y, (y-D) 

H'1(a/y)  =  An(a)e 

Y2(a,y)  =  Bh(a)sinh  y2Y  +  Ch(a)cosh  y2Y     (2-13) 

h  Yly 
*,(a,y)  -  D  e 


The  following  observation  about  the  solutions  in 

region  2  must  be  emphasized.   Any  wave  on  this  inhomogeneous 

waveguide  structure  will  partly  travel  through  air  and  partly 

through  the  dielectric  slab.   It  is  then  obvious  that 

k   <  8  <  k~   in  the  lossless  case  where   kn  =  k-,  =  k  . 
o         2  1     3     o 


,  „  v/7  2  Q2*    .,       2 

a   <  Vk0  -8- 


Since   -°°  <  a  <  °°  ,   for   |a|  <  V  k2  -6     then   y    <  0  . 
For  these  values  of  the  transform  variable  a  the  solutions 
in  (2-12)  and  (2-13)  for  region  2  will  change  from  hyperbolic 


I 2 — T 

to  trigonometric  functions  since   y?  =  JY2"  /   Yo"  =  ^^  2  ~a 

e  h 

The  eight  unknown  coefficients  A  (a)  through  D  (a) 

are  not  independent  but  can  be  related  to  each  other  through 
the  continuity  conditions  of  the  field  components  at  the 
interfaces  between  the  three  spatial  regions .   These  continuity 
conditions  in  the  spatial  domain  are  for 


y  =  0 


Et2(x,0,z)  =  Et3(x,0,z) 


Ht2(x,0,z)  =  Ht3(x,0,z) 


(2-14) 
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and  for  y  =  D 


£tl(x,D,z)  =  Et2(x,D,z) 


Ezl(x/D/z) 


Exl(xfD,z) 


*,<*>•*"      |x|  <J 

0  elsewhere 


j?x(x)e'       |x|  <  j 
0  elsewhere 


(2-15) 


Hzl(x,D,z)  -Hz2(x,D,z)  = 


W 


jv(x)eYZ 


Hxl(x,D,z)  -  Hx2(x,Dfz)  = 


0  elsewhere 


Dz(x)e'      |x|  >  j 
0  elsewhere 


The  subscript  t  is  an  abbreviation  for  tangential 

x  and  z  components  to  avoid  writing  essentially  identical 

equations  twice.  E         (x)  and  j    (x)  are  the  electric  field 
^  x,  z        Jx, z 

components  and  the  current  densities  in  the  x-z  plane  at 
y  =  D  in  region  1. 

Denoting  the  Fourier  transforms  of  the  x-  and  z- 
directed  current  density  and  electric  field  components  by 


Ev(a)  =  F{£  (x)}  E  (a)  =  F{2?  (x)} 

<«*■  J*  z  z 


Jv(a)  =  F{i  (x)}  J  (a)  =  F{j(x)} 

a  X  z  z 
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and  applying  these  transforms  and  the  transforms  of  the 

continuity  conditions  to  equations  (2-12)  and  (2-13)  an 

e  h 

elimination  of  the  unknown  coefficients  A   through  D  by 

substitution  is  possible.   The  details  of  this  lengthy 

algebraic  manipulation  are  given  in  Appendix  A.   The  results 

yield  a  set  of  coupled  equations  of  the  form 


N1(a/&)      N2(a,3) 
N3(a,B)      N4(a,S) 


E  (a) 
x 


Ez(a) 


Jx(a) 


Jz(a) 


(2-16) 


where  the  elements  of  the  N  matrix  can  be  thought  of  being 
the  Fourier  transforms  of  dyadic  Green's  function  components 

[Ref.  13]  by  an  application  of  the  convolution  integral  in 
the  space  domain.   The  lossless  case  y   =  jB   is  assumed  in 

(2-16)  . 

Nothing  seems  to  be  gained  so  far  since  neither  the 
electric  field  and  current  density  components  nor  their 
transforms  are  known.   The  coupled  equations  in  (2-16)  can 
however  be  further  simplified  by  the  methods  of  moments 

[Ref.  16],  applied  in  the  spectral  domain.   Define  the  scalar 
product  on  the  space  of  complex-valued  functions  of  the  real 
variable  a  over  the  domain   -°°  <  a  <  °°   according  to 
Reference  15  as 

CO 

<f(a),g(a)>  =   /   f(a)g*(a)  da 
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and  choose  as  a  suitable  set  of  weighting  functions  the 
Fourier  transforms  of  the  electric  field  components,  which 
means  Galerkin's  method  [Refs.  9,16].   The  moment  method 
applied  to  equation  (2-16)  leads  to 

<Nx(a,3)Ex(a) ,Ex(a)>  +  <N2 (a, 6) Ez (a) , Ex(a) >  =  0 

(2-17) 
<N3(a,6)Ex(a)  ,t"z(a)>  +  <N4  (a,B)  Ez  (a)  ,Ez(a)  >  =  0 

That  the  right  hand  side  is  zero  can  be  easily  verified  by 
Parseval's  relation  [Ref.  15],  observing  that  at  y  =  D  the 
presence  of  an  electric  field  component  excludes  the  presence 
of  a  current  density  component  and  vice  versa,  thus  electric 
field  and  current  density  are  orthogonal  to  each  other. 

Up  to  this  point  the  formulation  of  the  problem  is 
exact  since  no  approximation  has  been  made  yet.   Note  that 
the  electric  field  components  and  their  transform  although 
unknown  can  be  expanded  in  a  set  of  complete  basis  functions 
so  that  in  general 


a 

Vx)  ~     I        *n  exn(x) 


(2-18) 


and     E    (x)  =  I        be   (x) 

z        ,    n   zn 
n=l 


The  following  analysis  is  restricted  to  a  one  term  approxima- 
tion for  each  of  the  series  in  (2-18)  since  it  was  found  that 
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sufficient  accuracy  could  be  obtained  for  all  practical 
purposes.   So  n  =  1  in  equation  (2-18)  and  this  subscript 
will  be  omitted  from  here  on.  , 

With   E  (a)  z    F{e  (x)  }   and   E  (a)  s  F{e(x)} 

X  X  Z  £* 

(2-17)  then  becomes 

oo  oo 

a     f    N,  (a, 3)  |E    (a)  1 2da  +  b      /     N0(a,3)E    (a)E    *(a)da  =   0 

—  oo  — oo 

(2-19) 

00  CO 

a/     N  3(a,3)E    (a)E*(a)da  +  b      /     N    (a,  3)  |  Ez  (a)  |  2da  =   0 

—  GO  — OO 

The  dispersion  characteristic  of  the  slot  line  can  now  be 
found  by  varying  3  such  that  the  determinant  of  the  coefficient 
matrix  of  (2-19)  is  zero  for  a  given  set  of  physical  parameters 
at  a  desired  frequency  of  operation. 

The  remaining  question  is  what  kind  of  basis  functions 
to  choose.   The  choice  of  this  complete  set  of  basis  functions 
is  arbitrary  in  a  mathematical  sense,  since  as  long  as  this 
set  is  complete,  any  closed  form  of  the  field  components  can 
be  represented  by  it.   However,  the  rate  of  convergence  of 
the  series  representation  will  depend  on  how  well  the  first 
few  terms  approximate  the  closed  form.   This  in  general 
requires  some  a  priori  knowledge  of  the  true  distribution. 

In  order  to  determine  the  sensitivity  of  this  analysis 
to  the  choice  of  basis  functions,  an  investigation  of  various 
one  term  approximations  was  made.   One  choice  for  the  electric 
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field  components  which  requires  almost  no  insight  into  the 
field  problem  is  of  the  form 


•i"  « 


x   < 


W 


(2-20a) 


elsewhere 


.«"w    = 


•  2  <  x  <  0 
0  <  x  <  j 

elsev/here 


(2-20L) 


The  Fourier  transforms  of  these  basis  functions  are  then 


*<»<* 


2   .  ,aWs 


(2-21a) 


EzX)  (ot)  "  j  I  [1  "  cos(-r>i 


(2-21b) 


Another  choice  which  certainly  approximates  the 
field  more  closely  is 


42)  w 


tf 


,W,  2     2 
(y)   -  x 


x   < 


W 


(2-22a) 


elsewhere 


<2)  <*> 


x  V^f 


v2     2 

)   -  x 


i    W 
x|  <  2 


elsewhere 


(2-22b) 
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The  reason  that  (2-2  2)  are  superior  basis  functions  to 
(2-20)  is  that  (2-22b)  approaches  zero  smoothly  as 
|x|  -^y  and   |x|  -►  0  .   Furthermore  (2-22a)  approximates  the 
fields  more  closely  near  the  edge  of  the  slot  than  (2-20a) . 
The  Fourier  transform  of  (2-22)  can  now  be  found  in 
a  standard  table  of  Fourier  transforms  (e.g.  Ref.  17], 


E^2)  (a)  =  it  BQ(^)  (2-23a) 

E<2)(a)  =  j  ^,^B2(f)  (2-23b) 

where  BQ  and  B~  are  Bessel  functions  of  the  first  kind  of 
order  zero  and  two  respectively. 

Call  the  use  of  x  and  z  directed  field  components 
the  second  order  approximation  while  the  use  of  an  x  directed 
component  only  is  called  a  first  order  approximation.   Then 
a  reduction  in  computation  time  is  possible  using  this  first 
order  approximation  since  only  one  integral  in  (2-18)  has 
to  be  evaluated  instead  of  four  for  the  second  order 
approximation.   These  two  formulations  of  the  dispersion 
problem  were  programmed  into  a  digital  computer  language  and 
the  root  3  which  rendered  the  determinant  of  the  coefficient 
matrix  to  zero  was  found  numerically  by  an  iteration  scheme. 
Details  of  the  programming  are  given  in  Section  5.   First  and 
second  order  approximations  using  equations  (2-21)  and  (2-23) 
are  compared  in  Figure  2.1/1  with  results  from  Conn's  method 
in  Reference  12.   Although  no  indication  about  the  rate  of 
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cohn's  result 


PRESENT   METHOD 

—  Is1   ORDER,  eq.(2-20a) 
nd  ORDER,  eq.(2-20  ) 

* 1°'   ORDER,  eq.(2-22a) 

2nd0RDER,eq.(2"22) 


€r=  !6.0 
W/D  =  I.O 


OOI  0.03  0.05  0.07 

NORMALIZED     FREQUENCY     D/^ 


Figure  2.1/1   Comparison  of  Approximations  for  Dispersion 
Characteristic  of  a  Single  Slot  Line 
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convergence  and  hence  the  accuracy  of  the  two  series  for  a 
general  problem  can  be  obtained,  a  comparison  for  this 
particular  example  shows  that  the  largest  deviation  between 
the  different  approximations  is  less  than  4%.   As  already 
noted  above,  the  basis  set  (2-22)  is  vastly  superior  to 
basis  set  (2-20)  in  a  one  term  approximation  for  the  field 
components  due  to  the  physically  impossible  discontinuities 
of  equation  (2-20b) . 

A  numerical  investigation  with  regard  to  the  magnitude 
of  the  x  and  z  components  of  the  elctric  field  was  made,  using 
the  field  equations  (2-2  3)  in  the  second  order  approximation 
for  an  arbitrary  choice  of  W/D  =  0.2  and  e   =  20.   Up  to  a 
value  of  normalized  frequency   D/A  =  0.03   the  x  directed 
electric  field  component  was  found  to  have  a  magnitude  greater 
than  ten  times  that  of  the  z  directed  component.   This  provides 
justification  for  the  use  of  the  more  efficient  first  order 
approximation.   As  one  might  expect,  the  magnitude  of  the 
z  component  of  the  electric  field  increases  for  larger  values 
of  normalized  frequency.   However,  the  agreement  between  this 
and  Conn's  [Refs.  3,12]  method  is  excellent  for  a  first  order 
approximation.   The  theory  presented  here  is  a  verification 
of  Conn's  theory  and  in  a  sense  a  check  on  the  numerical 
analysis  by  Itoh  and  Mittra  [Ref.  8]  since  it  was  derived 
independently  with  only  their  very  brief  letter  for  guidance. 
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2.1.2   Characteristic  Impedance 

The  definition  of  the  characteristic  impedance  for 
an  ideal  transmission  line  is  uniquely  given  by  static 
quantities  because  these  lines  support  a  TEM  mode  wave.   On 
the  other  hand,  in  investigations  of  pipe-like  waveguides 
the  wave  impedance  is  widely  used  [Ref .  13] ,  however  the 
impedances  are  different  for  TM  and  TE  waves.   Since  the 
slot  line  is  neither  a  TEM  transmission  line  nor  does  it 
support  a  pure  TE  or  TM  waveguide  mode  but  rather  supports 
hybrid  modes,  no  unique  definition  of  the  characteristic 
impedance  can  be  found. 

Since  the  choice  of  definition  is  arbitrary  Cohn ' s 
very  practical  definition  in  Reference  3  is  used  in  the 
following  derivation  of  the  characteristic  impedance  for  a 
slot  line.   Define 

V2   . 
Z0  =  j-pS-  (2-24) 

avg 

where  V  is  the  slot  voltage  and  P    is  the  time  average 
power  flow  on  the  slot  line  which  according  to  Reference  13 
may  be  computed  as 


P 


=  i  Re{  //  E\  xh  *  .  a    dxdy }  (2-25) 

2      _    t    t     z     J 


avg         s 


Observe  that  although  the  hybrid  mode  can  be  fully  described 
as  a  superposition  of  TE  and  TM  modes,  the  evaluation  of  the 
time  average  power  flow  is  not  just  a  superposition  of  the 
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single  time  average  power  flow  carried  by  each  single  mode. 
A  simple  manipulation  unveils  that  in  addition  to  the  single 
TE  and  TM  power  terms  there  are  cross  terms  as  shown  below 
in  (2-26).   Assume  y  =   j3  f      which  means  for  z  dependence 
eJ    a  wave  traveling  in  the  -z   direction. 

The  transverse  field  components  of  the  integrand 
in  (2-2  5)  can  be  expressed  in  terms  of  the  scalar  potential 
functions  as  shown  in  equations  (2-3)  through  (2-6) .   Equation 
(2-2  5)  becomes  with  this  substitution  and  dependence  of  the 
scalar  potentials  on  x  and  y  understood 


•    -  |  Re  /  /   {.,,BI(H>2+  0h    ♦  u£g[(M)2+  (||)2i  + 

-*  —  O J   —  CO  "*  "* 


w2+Jt2)(M-|i_  |1||);  dy  dx    (2_26) 


3y  3x    3x  '6y 

The  limits  of  integration  are  infinite  since  the  slot  line 
is  an  open  boundary  structure,  which  makes  the  use  of 
Parseval's  relation  [Ref.  15]  feasible.   By  this  method 
(2-26)  is  transformed  into  the  spectral  domain  where  again  a 
double  integral  is  obtained  with  the  variables  of  integration 
a  and  y  instead  of  x  and  y. 


00    oo 


p 


-     1     o«        r       r     r       2o       I  *  1 2  ,        |3V|2         D        , 3  $ , 2 

=  7—  Re     J      J    l"a   3(o£|<J>|      -   3wy|-r~|      -   Puel-srrl 


*vg        4tt   -      '      '     '   "   —  «'i  —  ir.yi  —  i3y 


—  CO      «»oo 


-v 


2Q ivi2    _     ... l02,.    3J'      A    3* 


* 


^  «             I   in  t  ^  •        r-  *■    t   ,        0  1             ,        0  4-  .,,      \ 

-a    3a)y  |  V  J       -  3aB    ($  j—  +  |y  ¥    ) 

+    j«k2(4>*   |I  +   |i-  H')  ]    da   dv         (2-2  7) 

J                  3y  3y 
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where  the  dependence  of  the  scalar  potentials  on  variables 
a  and  y  is  implicitly  understood.   Besides  any  dependence 
on  physical  parameters  and  the  frequency  of  operations,  these 
scalar  potentials  are  furthermore  dependent  on  the  propagation 
constant  3,  which  is  known  from  Section  2.1.1,  and  on  the 
assumed  field  distributions  as  given  by  equations  (2-12) , 
(2-13)  and  the  detailed  derivation  in  Appendix  A.   Note  that 
the  functional  dependence  on  the  variable  y  may  be  removed 
by  analytical  integration,  which  changes  the  double  integral 
to  a  single  integral.   For  a  numerical  integration  on  a 
digital  computer  the  integral  (2-2  7)  has  to  be  evaluated 
separately  for  the  three  spatial  regions  since  according  to 
equations  (2-12)  and  (2-13)  the  integrand  assumes  different 
forms  there.   The  integrands  for  the  three  regions  are  given 
in  Appendix  B  in  equations  (B-l)  through  (B-4) . 

It  should  be  noted  that  the  amount  of  algebraic 
complexity  in  these  equations  although  still  quite  large 
could  be  decreased  for  numerical  purposes  by  neglecting  the 
z-directed  electric  field  component  which  means  that  the 
spatial  phase  shift  between  E   and  E   required  no  attention 
in  this  algebraic  derivation.   It  is  shown  in  the  comparison 
with  Cohn ' s  method  [Refs.  3,12]  that  this  use  of  a  first 
order  solution  results  in  good  agreement. 

A  computer  program  was  developed  which  first  computes 
the  dispersion  characteristic  and  then  the  average  power  flow 
where  the  electric  field  was  approximated  by  (2-22a) .   To 
solve  for  the  characteristic  impedance  in  equation  (2-24)  , 
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N 


the  slot  voltage  has  to  be  computed.   Since  the  slot  voltage 

W 
2 

V  =    /   E  (x)  dx   ,   this  integration  may  be  accomplished 

O  X 

w 

~2 

analytically;  depending  on  the  assumed  electric  field 

approximations   V  =  tt   for  (2-22a)  for  instance.   A  comparison 
of  results  obtained  by  this  formulation  of  the  characteristic 
impedance  with  values  obtained  by  Cohn's  method  in  Reference 
12  is  given  in  Figure  2.1/2.   The  very  close  agreement  for 
the  two  arbitrarily  chosen  sets  of  parameters  substantiates 
the  accuracy  of  both  methods  and  is  the  first  known  verifica- 
tion of  Cohn's  method. 

For  comparative  purposes  it  is  noted  that  the  field 
distribution  (2-20a)  resulted  in  lower  values  for  the 
characteristic  impedance  by  not  more  than  10%  from  those 
obtained  from  (2-22a) . 

A  byproduct  of  this  analysis  is  a  consideration  of 
the  power  distribution.   The  spatial  division  of  the  slot 
line  in  this  numerical  formulation  allows  an  analysis  of  the 
amount  of  power  flowing  in  each  of  these  three  regions.   In 
Figure  2.1/3  the  ratio  of  the  power  in  each  region  over  the 
total  power  carried  on  the  structure  is  shown  as  a  function 
of  normalized  frequency  for  two  different  permittivities. 
Several  facts  can  be  learned  from  these  graphs : 

a)  The  amount  of  power  carried  in  each  region  is  frequency 
dependent.   As  the  frequency  increases  more  power  is 
carried  through  region  2,    the  dielectric  region. 

b)  Only  a  small  amount  of  the  total  power  is  found  in 
region  3. 

c)  More  power  is  carried  through  the  air  regions  if  the 
permittivity  is  lowered. 
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Figure   2.1/2      Comparison    for   Characteristic   Impedances 
of   a  Single    Slot   Line 
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Figure  2.1/3   Time  Average  Power  Ratio  in  Region  1,  2  and  3 

vs.  Normalized  Frequency  for  a  Single  Slot  Lino 
with  W/D  =  1.0. 
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These  results  may  have  been  known  previously  but  for  the 
first  time  are  here  quantitatively  analyzed  by  the  means  of 
the  spectral  domain  formulation. 

The  following  investigation  indicates,  however, 
that  a  slight  error  due  to  the  approximation  in  the  first 
order  exists.   This  example  serves  more  illustrative  purposes 
but  might  perhaps  stimulate  further  research.   With  the  above 
formulation  of  the  time  average  power  flow  a  graphical  display 
of  the  power  density  for  the  slot  line  may  be  obtained.   For 
an  arbitrarily  chosen  set  of  parameters  the  power  density 
in  the  spectral  domain  is  shown  in  Figure  2.1/4a  versus  the 
normalized  transform  variable  aD  for  the  spatial  region  1 
above  the  substrate.   The  unit  on  the  ordinate  is  omitted 
since  it  is  insignificant  in  this  illustration. 

From  Cohn ' s  zero  order  analysis  in  Reference  3  it  is 
obvious  that  the  power  density  is  the  largest  close  to  the 
slot  and  it  is  also  obvious  that  it  does  not  change  its  sign. 
To  verify  this  the  Fourier  inversion  of  the  spectral  domain 
power  density  has  to  be  taken  to  obtain  the  power  density  in 
the  space  domain.   Since  the  analysis  is  numerical  in  nature 
inversion  by  a  discrete  Fourier  transform  technique  (DFT)  on 
a  digital  computer  seemed  appropriate.   The  result  is 
displayed  in  Figure  2.1/4b.   Clearly  the  power  is  the  largest 
close  to  the  slot  and  decays  rapidly  but  contrary  to  any 
expectation  changes  its  sign.   It  can  be  easily  verified  that 
this  behavior  is  due  to  the  small  dip  at  aD  =  0  for  the  power 
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Figure  2.1/4   Power  Density  Distribution  in   a)  Spectral, 

b)  Space  Domain  for  Region  1  of  a  Single  Slot 
Line  with  e   =20,  W/D  =1.0  and  D/A  =  0.0  2 
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density  distribution  in  the  spectral  domain.   This  dip  is 
believed  to  originate  from  the  first  order  approximation  of 
the  electric  field.   Although  a  similar  phenomenon  can  be 
observed  in  an  investigation  of  the  power  densities  in  regions 
2  and  3  it  is  seen  that  the  error  introduced  is  quite  small. 

2.2   COUPLED  SLOT  LINES 

With  two  transmission  lines  placed  adjacent  to  each 
other  such  that  electromagnetic  coupling  effects  the  fields 
on  both  line  structures,  the  design  of  directional  couplers 
and  filters  is  possible.   Figure  1.1/2  shows  two  parallel 
slot  lines  whose  analysis  is  the  subject  of  this  section. 
There  are  many  ways  to  describe  the  coupling  effect  between 
two  transmission  lines.   According  to  Reference  18,  which  is 
a  compilation  of  some  twenty-nine  papers  on  coupled  trans- 
mission lines,  the  most  common  formulation  is  to  define  the 
natural  modes  of  the  guided  wave  structure.   Since  in  this 
analysis  two  slots  of  equal  width  and  constant  separation  are 
considered,  define  the  two  natural  modes  as  the  even  and  odd 
mode  such  that 


E  (x)  =  E  (-x)       for  the  even 


and     E  (x)  =  -E  (-x)      for  the  odd  mode, 
x        x 


The  existing  theory  developed  so  far  in  Sections  2.1.1  and 
2.1.2  can  now  be  easily  extended  for  the  analysis  of  line 
parameters  for  these  two  modes  of  the  coupled  slot  lines. 
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It  is  observed  that  the  dyadic  Green's  function  transforms 
in  equation  (2-16)  are  functions  of  a,  (3  ,  c    ,    D  and  the 
operating  frequency  but  are  independent  of  the  slot  configura- 
tion.  The  physical  dimensions  of  the  slot  enter  into  the 
calculation  only  through  the  assumed  field  distributions  or 
basis  functions.   Thus,  it  is  only  necessary  to  modify 
equations  (2-20)  and  (2-22)  such  that  the  mathematical  formula- 
tion of  the  basis  functions  corresponds  to  the  physical 
configuration  and  field  distributions  of  the  coupled  slots. 
As  an  example  Figure  2.2/1  shows  the  assumed  electric  field 
distribution 
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0  elsewhere 

The  modification  for  the  z-directed  electric  field 
component  is  similar.   The  following  short  manipulation  leads 
to  the  Fourier  transform  of  this  distribution: 
The  shifting  property  of  Fourier  transforms  states  that  if 

00 

F{<{>(x)}  =  »(a)  =   /   <j>(x)e3ax  dx  then   F{<J>(x-xn)}  =  eDax0$(a) 
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Figure  2.2/1   Assumed  Electric  Field  Component  in  x-direction 

vs.  x  for  Even  and  Odd  Mode  on  Coupled  Slot  Lines 
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The  Fourier  transform  of  (2-22a)  is  given  in  (2-2 3a)  as 

(2) 
E    (a) .   Then  the  Fourier  transforms  of  the  electric  field 

in  (2-28)  become  for  the  even  mode 

•   S+W       .   S+W 
E^e  (a)  "  (e       +   e     *  )  E^'  (a) 


=  2  cos(a  ~)  •  E(2)(a)  (2-28a) 


and   for   the   odd  mode 


E(2)(a)    =   +J2    sin(a  —-)-  E(2)  (a)  (2-28b) 

xo  J  2  x 


Another  change  to  adapt  the  formulation  of  Section  2.1 
for  the  coupled  line  analysis  is  necessary  in  equation  (2-24)  ; 
since  the  total  time  average  power  surrounding  the  transmission 
lines  is  now  due  to  two  lines 

V  2 

0 
— (2-29) 


0o,e   P 

avg  o , e 

for  the  odd  and  even  mode  respectively. 

Beyond  these  no  further  major  modifications  are  necessary 
to  use  the  previously  developed  theory  and  computer  programs. 
In  the  following  discussion  only  results  using  a  field 
approximation  (2-2 8a)  and  (2-2 8b)  are  presented.   There  are 
two  reasons  for  this:   a)  the  approximation  of  equations 
(2-22)  was  found  to  be  superior  to  (2-20)  in  Section  2.1 


53 


b)  the  use  of  a  second  order  approximation  became  prohibitive 
for  the  dispersion  characteristic  of  coupled  slots  since  the 
amount  of  computer  time  increased  tremendously  due  to  the 
necessity  of  using  a  more  sophisticated  numerical  integration 
scheme  (details  in  Section  5).   Besides,  the  time  average 
power  computation  uses  only  a  first  order  approximation 
anyway . 

Figures  2.2/2a-b  show  new  an  example  of  the  numerical 
results  for  the  dispersion  characteristics  and  the  charac- 
teristic impedances.   These  results  are  basically  new  in  the 
literature  and  only  a  partial  comparison  with  other  results 
is  possible.   A  quantitative  comparision  e.g.  between  these 
odd  mode  line  parameters  and  the  quasi-static  parameters 
calculated  by  Davis  et  al .  in  Reference  7  was  made  and  is 
shown  in  Figure  2.2/3  for  W/D  =  l.C  —  or  in  the  notation  of 
Reference  7,  t  =  S  —  at  frequencies  of  1,  3  and  5  GHz.   Good 
agreement  is  found  for  the  lower  frequency  range.   Note  that 
the  characteristic  impedance  of  a  CPW  is  one  half  of  Z^   (odd 
mode  characteristic  impedance)  for  equal  slot  widths. 

A  qualitative  check  on  these  results  can  be  made  by 

investigating  the  even  mode  characteristic  impedance  in  the 

limiting  case  as  the  separation  between  the  slots  becomes 

very  small.   One  expects  that  the  even  mode  characteristic 

impedance  Z~   will  be  close  to  one  half  of  Zn,  where  Zn  is 
Ue  u         u 

the  characteristic  impedance  of  a  single  slot,  whose  width 
is  twice  the  width  of  the  slot  in  the  coupled  structure. 
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Figure  2.2/2a  Even  and  Odd  Mode  Dispersion  Characteristics 

for  Coupled  Slots  with  W/D  =0.4    and  e   =  16 
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Figure  2.2/2b  Even  and  Odd  Mode  Characteristic  Impedances 

for  Coupled  Slots  with  W/D  =0.4  and  e  =  16 
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and  e   =10 


57 


The  dispersion  characteristic  for  both  structures,  however, 

should  be  the  same.   This  comparison  for  the  impedances  is 

indicated  in  Figure  2.2/2b  and  again  shows  an  excellent 

agreement. 

It  is  interesting  to  observe  the  coupling  and  decoupling 

between  the  waves  in  the  two  slots  as  the  frequency  varies 

for  larger  values  of  the  separation  constant  S/D.   With 

increasing  frequency  the  waves  become  more  closely  bound 

to  the  slot  which  means  that  there  is  less  interaction 

between  the  two  waves.   In  this  case  ZA   and  ZA   converge 

ue      Oo       ^ 

to  Z0,  the  characteristic  impedance  of  a  single  slot  with 
no  coupling  effect. 

Another  interesting  phenomenon  is  the  fact  that  for  a 
fixed  frequency  the  wavelength  ratio  \'/\   in  the  even  mode 
first  increases  and  then  decreases  as  S/D  is  increased 
from  very  small  to  larger  values  as  shown  in  Figure  2.2/2a. 
An  explanation  may  be  given  as  follows.   For  small  enough 
values  of  S/D  the  metal  strip  between  the  slots  has  little 
effect  on  the  propagating  wave  and  the  wave  propagates  as 
if  it  were  in  a  slot  of  width  2  •  W/D  +  S/D.   Increasing  the 
separation  between  the  slots  effectively  increases  the  slot 
width  so  that  the  metal  strip  still  has  little  effect  and 
the  wavelength  ratio  increases.   As  the  separation  continues 
to  increase,  the  two  waves  start  to  decouple  and  behave  more 
as  two  waves  on  two  slot  lines  which  will  finally  be  totally 
decoupled.   Each  wave  then  propagates  on  a  slot  line  with 
width  W/D,  hence  \*/\    decreases. 
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However,  again  an  error  due  to  the  first  order  approxima- 
tion can  be  observed.   Note  that  for  larger  values  of  S/D 
the  wavelength  ratios  and  the  characteristic  impedances  for 
the  two  natural  modes  should  merge  evenly  together  towards 
the  respective  line  parameters  for  a  single  slot  line.   It 
is  observed,  however,  that  for  still  distinct  wavelength 
ratios  for  the  two  modes  the  characteristic  impedances  for 
these  modes  merge  not  only  together  but  overcross  each  other 
by  a  small  amount.   This  fact  is  not  shown  in  Figures 
2.2/2a-b  for  clarity  purposes  but  is  presented  in  Figure  2.2/4 
where  for  one  arbitrarily  chosen  value  of  normalized  frequency 
and  slot  width  the  two  impedances  for  the  two  modes  are 
displayed  as  a  function  of  separation.   Obviously  the  error 
introduced  is  not  very  large  in  this  example  but  tends  to 
increase  somewhat  for  higher  frequencies.   Since  the  numerical 
methods  used  in  this  analysis  were  carefully  checked  it  is 
believed  that  this  error  is  caused  by  neglecting  the 
z-directed  component  of  the  electric  field  in  the  power 
calculation.   As  noted  in  Section  2.1.1  this  component 
increases  with  frequency  and  the  approximation  caused  a 
slight  error  in  the  power  calculation  on  a  single  slot  line 
in  Section  2.1.2. 

However  the  partial  comparison  of  the  results  with  other 
data  given  above  and  an  experimental  verification  described 
in  Section  4  confirms  the  accuracy  of  the  first  order 
approximation  for  the  frequency  ranges  considered  here. 


Z    For  single  slot 
o 
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Figure   2,2/4      Even   and  Odd  Mode   Characteristic   Impedance 
vs.    Separation    for   D/X   =    0.03,    W/D   ==0.4 


and   e      =   16 
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2.3   COPLANAR  STRIP  LINE 

The  configuration  of  a  coplanar  strip  line  is  shown  in 
Figure  1.1/4  and  consists  of  two  metal  strip  conductors 
placed  on  one  side  of  the  dielectric  substrate.   The  spacing 
between  these  two  strips  is  assumed  to  be  constant  in  the 
z  direction.   The  permittivity  of  the  substrate  confines 
the  fields  to  regions  around  the  conductors  to  minimize 
radiation  losses  and  to  make  this  line  suitable  for 
integrated  circuit  applications.   A  hybrid  mode  analysis  in 
the  spectral  domain  is  again  used  to  obtain  the  two  important 
line  parameters.   However  the  key  for  a  solution  by  this  method 
is  that  the  investigator  has  to  have  some  idea  of  how  the 
fields  on  this  structure  behave.   In  the  study  of  the  slot 
line  an  approximation  of  the  electric  field  across  the  slot 
was  required  and  was  not  difficult  to  obtain  by  physical 
reasoning. 

However  for  the  coplanar  strip  line  an  approximation  of 
the  current  density  across  each  strip  is  more  feasible. 
Hence  a  set  of  coupled  equations  similar  to  (2-16)  is 
desirable  with  the  roles  of  current  density  and  electric 
field  interchanged.   Note  that  this  inversion  was  already 
obtained  in  the  development  of  (2-1G)  and  is  the  equation 
(A-22)  in  Appendix  A  give  here  as 


M-^a^)       M2(cc,3) 


M3 ( a , 3 ) 


M4(a,3) 


J  (a) 
x 


Va). 


E  (a) 
x 


Ez(a) 


(2-30) 
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According  to  the  scalar  product  definition  and  Galerkin's 
method  given  in  Section  2.1.1  a  suitable  set  of  weighting 
functions  are  the  transforms  of  the  current  density  distribu- 
tions.  The  methods  of  moments  applied  to  (2-30)  and  an 
expansion  of  the  current  densities  in  basis  functions  yields 
the  second  order  approximations  similar  to  (2-19) 


GO  WW 

c   /  M,  (a,  3)  I J  (a)  I  da  +  d   /  M9(a,3)J  (a)  J,  *  (a)  da  =  0 


(2-31) 


c   /  M_(a,3)J  (a) J  *(a)da  +  d   /  M.(a,B) |J  (a) |2da  =  0 


In  the  detailed  analysis  of  different  approximations  in 
Section  2.1.1  it  was  found  that  a  first  order  approximation 
using  only  E  (a)  was  sufficient.   Since  the  necessary  power 
calculation  for  the  computation  of  the  characteristic 
impedance,  will  only  use  a  first  order  approximation,  this 
approach  is  used  throughout  the  analysis  of  coplanar   strips 
Thus  J  (x)  ~  0   and  (2-31)  reduces  to  the  problem  of  finding 
a  propagation  constant  3  such  that 


/   M4(a,B) | J2(a) |2da  =  0  (2-32) 


Again  one  has  the  task  here  to  choose  a  proper  basis 
function.   The  following  functions  are  reasonable  basis 
functions  for  the  current  density  component  j_(x)  over  the 
strip  conductors  and  have  been  used  extensively  in  the 


62 


analysis  of  microstrip  transmission  lines.   According  to 
Figure  1.1/4  the  left  strip  extends  for  x  such  that 
-S/2  -  W  <  x  <  -S/2   and  the  right  strip  for 
S/2  <  x  <  S/2  +  W  respectively.   Then  the  three  possible 
basis  functions  are 


-1 
D^1}(x)  -j+1 

0 


left 

right 

elsewhere 


(2-33a) 


-1 


3<2!<*> 


,W, 2    ,   ,  S+W. 2 
(_)   -  (X  +  --J-) 


,W.  2    ,     S+W,  2 


left 


right 


(2-34a) 


elsewhere 


Jj3)w 


_1  - 


1  + 


X 

+  i 

>+W 

2 

W 
2" 

x  - 


S+W 


W 
2" 


left 


right 


(2-35a) 


0  elsewhere 

A  schematic  presentation  of  these  basis  functions  is  given 
in  Figure  2. 3/1. 
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Figure  2.3/1 


Assumed  Current  Density  Component  in  z~Direction 
vx.  x  for  Coplanar  Strip  Line 

a)  equation  (2-33a) 

b)  equation  (2-34a) 

c)  equation  2~35a) 
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The  Fourier  transforms  can  now  be  easily  found  by  using 
again  the  shifting  property  of  the  Fourier  transformation 
indicated  in  equations  (2-28) .   The  transforms  are  given  as 

K      (<*)  =  3  77  sin  [a  -~— ]  •  sm(-^-)  (2-33b) 


z        ■*  a 


T' 


t(2)  /   x         o      •   r    S+W,      „   ,aWx  ,„   0„iv 

Jv  '(a)  =3  2tt  sin  [a  — ~— ]  •  Bn  (-*-)  (2-34b) 

T(3)  ,  ,    .8   .  ,   S+W,    T    24      (aW)2-12   .  ,aW. 
J2   («)  =  D  77  sin  [a  -^—]  •  [ ~  +  5—  sin  (-5-) 

+  3  (aW)2-8  C0S(^}]   (2_35b) 
(aW)J 


respectively.   (BQ  is  the  zero  order  Bessel  function  of  the 
first  kind. ) 

It  is  obvious  now  that  the  theory  and  the  numerical 
analysis  of  the  slot  line  dispersion  characteristic  needs 
only  minor  changes  to  adapt  a  solution  scheme  for  a  coplanar 
strip  line  dispersion  characteristic. 

To  compute  the  characteristic  impedance  of  this  trans- 
mission line  a  reformulation  of  (2-24)  is  necessary  since 
the  voltage  between  the  two  strips  is  not  explicitly  known. 
Again  this  definition  is  somewhat  arbitrary  because  of  the 
non-TEM  nature  of  the  problem. 

Define 

2  P 
ZQ  =     a^g  (2-36) 

I0 
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where  Ifl  is  the  total  current  on  one  strip  and  can  be  found 
analytically  through  an  integration  of  the  current  densities 
(2-33a)  through  (2-35a)  with  respect  to  the  variable  x, 
yielding 


I^    =  W  (2-37) 

I(j2)  =  7T  (2-38) 

I0(3)  =  |  W  (2-39) 


The  computation  of  the  time  average  power  flow  follows  in 
general  the  method  outlined  in  Section  2.1.2.   Some 
modifications  are  necessary,  however,  and  are  shown  in 
detail  in  Appendix  B. 

The  results  for  the  dispersion  characteristic  and 
characteristic  impedance  with  these  three  different 
approximations  are  shown  in  Figure  2.3/2a-b. 

These  line  parameters  are  basically  new.   The  wavelength 
ratio  A1 /X   was  measured  by  Luna  [Ref.  10]  and  good  agreement 
is  found  between  the  present  theory  and  the  experimental 
data.   Preliminary  measurements  on  the  characteristic 
impedance  with  a  time  domain  ref lectometer  indicated  good 
agreement  with  the  theory  using  the  current  density  distribu- 
tion (2-34a) .   This  work  should  be  extended  in  the  future. 

Reasonable  agreement  for  the  impedances  may  be  found  by 
comparing  the  present  values  with  the  quasi-static  results 
by  Wen  [Ref.  4].   As  one  might  expect,  however,  the  present 
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method  yields  somewhat  larger  values  due  to  the  finite 
thickness  of  the  dielectric  substrate. 

2.4   PERTURBATIONAL  ANALYSIS  OF  A  SINGLE  SLOT  LINE 
ON  A  FERRITE  SUB ST FATE 

The  hybrid  mode  analysis  presented  so  far  has  been  applied 

to  coplanar  transmission  lines  on  a  dielectric  substrate. 

Since  these  lines  are  new  and  only  one  complete  frequency 

dependent  analysis  has  been  proposed  it  is  not  surprising 

that  almost  no  analysis  has  appeared  yet  for  coplanar 

transmission  lines  on  ferrite  substrates. 

The  reasons  why  these  lines  offer  interesting  properties 

are  obvious  and  will  be  shortly  given  here.   Ferrites  are 

magetic  oxides  with  high  resistivity,  a  scalar  dielectric 

constant  and  a  tensor  permeability  at  microwave  frequencies. 

Electromagnetic  fields  within  these  ferrites  will  interact 

with  the  precession  of  electron  spins,  where  this  precession 

is  related  to  the  mage fixation  vector  which  also  precesses 

about  the  constant  biasing  magnetic  field.   The  magnetic  bias 

field  determines  the  orientation  and  the  frequency  of  the 

precession  and  thus  can  be  adjusted  to  any  desired  microwave 

frequency.   Depending  on  the  relative  polarization  and 

orientation  of  the  RF  magnetic  field  and  the  magnetization 

precession  to  each  other,  stronger  or  weaker  interaction 

is  obtained.   Since  the  amplitude  of  the  precession  decreases 

due  to  damping  the  amount  of  energy  absorption  can  be  quite 

large  for  resonance.   A  simple  example  of  the  use  of  this 


68 


phenomenon  in  fer rites  is  the  isolator.   If  the  RF  magnetic 
field  is  circularly  polarized  and  has  the  same  orientation 
as  the  electron  spin,  strong  interaction  is  obtained  and  a 
maximum  of  energy  is  absorbed.   This  attenuation  would  be  a 
maximum  for  a  wave  travelling  in  one  direction  —  where  the 
orientation  of  RF  field  and  electron  spin  is  aligned  —  and 
smaller  for  a  wave  travelling  in  the  opposite  direction. 
Thus,  the  amount  of  absorption  is  dependent  on  the  propagation 
direction  which  provides  non  reciprocal  attenuation.   Other 
non-reciprocal  devices  are  switches,  phase  shifters,  filters, 
and  circulators  [Ref.  19,  20]. 

It  is  then  seen  that  coplanar  transmission  such  as  single 
or  coupled  slot  lines  —  or  CPW  —  are  bound  to  be  used  in 
ferrite  applications  since  these  lines  posess  a  region  of 
elliptically  or  circularly  polarized  magnetic  RF  field,  much 
more  pronounced  than  the  quasi  TEM  microstrip  structure.   A 
schematic  representation  of  the  fields  on  a  slot  line  are 
shown  in  Figure  2.4/1.   For  these  reasons  and  the  fact  that 
coplanar  transmission  lines  offer  advantages  over  the  micro- 
strip  line  as  stated  in  Section  1,  it  can  be  foreseen  that 
nonreciprocal  devices  for  integrated  circuits  will  be  demanded, 
The  need  for  an  analysis  of  these  structures  becomes  apparent. 
Furthermore,  since  the  coplanar  transmission  lines  as  well 
as  the  properties  of  the  ferrites  are  highly  frequency  depen- 
dent any  analysis  of  these  structures  must  be  also  frequency 
dependent.   In  this  section  sn   attempt  is  made  to  analyze  the 
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Figure  2.4/1   Field  Configuration  on  Slot  Line 


70 


slot  line  on  a  ferrite  substrate  by  a  perturbational  method 
for  the  propagation  constant  in  the  spectral  domain.   This 
method  can  then  be  easily  extended  for  the  analysis  of 
coupled  slots  and  coplanar  strip  lines  by  the  means  described 
in  Sections  2.1-2.3. 

2.4.1   Perturbation  Expression  for  the  Propagation  Constant 

Perturbational  analyses  are  useful  for  the  calcula- 
tion of  changes  in  some  quantity,  if  the  problem  itself  is 
slightly  changed  or  perturbed.   Usually  a  solution  for  the 
unperturbed  problem  is  known  and  the  solution  of  the  perturbed 
problem  is  found  by  assuming  that  these  disturbances  are 
small.   Examples  for  the  use  of  the  perturbation  theory  are 
given  in  References  19  and  21. 

In  the  present  application  the  unperturbed  problem 
is  the  computation  of  the  wavelength  ratio  A  * /A  —  or  the 
propagation  constant  —  of  the  slot  line  on  a  dielectric  sub- 
strate.  The  solutions  to  this  problem  are  known  and  given 
in  Section  2.1.1.   For  the  perturbed  problem  the  dielectric 
substrate  is  replaced  by  a  ferrite  substrate  with  the  same 
dielectric  constant.   It  is  assumed  that  this  change  in  the 
structure  is  small  enough  to  allow  a  perturbational  analysis 
of  the  change  in  the  propagation  constant. 

The  perturbational  expression  for  the  propagation 
constant  due  to  small  changes  in  the  type  of  material  within 
the  guidance  structure  is  given  by  Helszajn  in  Reference  19 
as 
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ju  //  [(en [AX®] •£•£*)  +  (yn[AXe] -H-H*)]ds 
ac     u    e 

y'+Y*  =  — - (2-40) 

//   (E*  x  H  +  E  x  H*) -a   ds 
S  Z 


where  the  following  nomenclature  is  used 

S       total  cross  section  of  waveguide 

AS      cross  section  of  waveguide  which  is  perturbed  by 
a  change  in  type  of  material 

E  and  H  are  the  electric  and  magnetic  field  of  the  unperturbed 
structure 

[Axe]    change  in  magnetic  susceptibility  tensor 

[Aye]    change  in  electric  susceptibility  tensor 

Y1  =  a'+j3'   propagation  constant  of  perturbed  problem 

Y  =  a+ j 3     propagation  constant  of  unperturbed  problem. 

Although  (2-40)  applies  in  Reference  19  strictly  to 

a  closed  boundary  waveguide  it  can  be  shown  that  this 

expression  also  holds  for  an  open  boundary  structure  such 

as  the  slot  line.   This  derivation  is  shown  in  Appendix  C. 

A  particular  application  of  (2-40)  for  the  slot  line  is  now 

possible  and  is  derived  in  the  next  section. 

2.4.2   Computation  of  Perturbation  Expression 
in  the  Spectral  Domain 

From  the  field  configurations  on  the  slot  line  it  is 

seen  that  the  magnetic  field  is  elliptically  polarized  within 

the  y-z  plane.   The  strongest  interaction  between  the  RE 

magnetic  field  and  the  magnetization  of  the  ferrite  can  then 

be  expected  for  a  magnetic  biasing  field  in  the  x-direction. 
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From  the  discussion  of  the  susceptibility  tensor  in  an 
ellipsoidal  medium  in  Reference  19  it  is  seen  that  the 
present  configuration  is  similar  to  that  of  a  thin  disk. 
Since  the  biasing  magnetic  field  is  along  the  x-axis  no 
demagnetization  for  this  field  has  to  be  considered. 
However,  a  demagnetization  for  the  y-component  of  the  RF 
magnetic  field  is  to  be  noted.   This  was  the  deciding  factor 
in  using  the  external  rather  than  internal  susceptibility 
tensor,  where  the  demagnetization  enters  into  the  formulation 
of  the  tensor  rather  than  the  fields  and  was  found  to  be  much 
easier  to  use  in  the  following  derivation.   According  to  the 
small  signal  approximation  of  the  equation  of  motion  in 
Reference  19  two  radian  frequencies  can  be  defined  where 


to„  =    tort    -  N    to      +  N,u 
y  0  x  m         y  m 

w      =   io~    -  N    to      +  N   to 
z  0  x  m  z  m 

MS 
with      coA    =   -Y   H  ,         to      =   -   Y  -— 

0  '       xo       '  m  ^o 

Y   =    -1.4    •    g    •    106     [Hz/Oe] 

g    Lande"  factor 

H    magnetic  biasing  field  in  x  direction 

M_    saturation  magnetization 

N    demagnetization  factor  (in  this  application  N   =  N   =  0, 
N   -  1)  x     z 
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Furthermore ,  define  the  ellipticity  of  the  normal  modes  of 

/to — ' 


the  uniform  precession  as   e   = 


V 


and  note  that 


b) 


Voo  •  co     is  the  Kittel  resonance  relation.   According 
Y   z 

to  the  Landau-Lif shitz  phenomenological  damping  formulation 

Y  /  H 
[Refs.  19,22]  the  damping  can  be  defined  as   acu  =  -  -*-£- 

where  AH  is  the  linewidth  of  the  ferrite  material.   The 

external  susceptibility  tensor  is  then 


lxel 


0 

- 

0 

Xyy 

*yz 

xzy 

xzz 

(2-41) 


where   x  =  x'~"Jx"    and   Xyz  =  "Xzy 


2  2S  2  2 

co  co  (co„  -co  )  +  co  o)  .co  a 


Ayy     O 


m  r  r 


m  r 


DN 


(2-41a) 


tt„w  [co„  +  co  (1+a  )  ] 

y"   =  e         — 

*yy    o         DN 


(2-41b) 
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2  xyy 
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(2-41c) 


(2-41d) 


(2~41e) 


7z 


3    zVrw  a 

~dn 


(2-41f) 


?  ?  7        ?        ?  ?  ? 

with     DN   =    [cor      -   co    (l+a^)3      +    4co     co  a 
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Equation  (2-40)  can  now  be  evaluated.   Note  that 

the  denominator  can  be  expressed  as  4-P  (see  also 

c  avg 

Section  2.1.2)  and  that   [Ax']  =  0   due  to  the  assumption 

that  the  dielectric  constant  remains  unchanged,   [Ax  ] 

is  obviously  identical  to  equation  (2-41) .   To  evaluate  the 

remainder  of  the  numerator  note  that  the  product 

[Ax  ]»H«H*   has  to  be  integrated  in  the  crosssectional  area 

AS  which  extends  such  that   -00  <  x  <  °°   and   0  <_  y  <_  D  . 

The  RF  magnetic  field  is  given  in  terms  of  the  scalar 

potentials  in  equations  (2-2) ,  (2-5)  and  (2-6) .   With  these 

substitutions  (2-40)  becomes 


avg   -°°  0      JJ  -r 


+  2x   k  ?ii>^(+j3  ~ +  joj£.,  ■* — ) 

Ayz  c2^2   J   3y     J   2  3x 


+  X   k  n  {$<*)     3  dvdx 

AZZ  CZ        2 

(2-42) 


The  upper  and  lower  signs  denote  waves  traveling  in  the 
negative  and  positive  z-direction  respectively. 

This  double  integral  is  transformed  via  Parseval's 
relation  [Ref.  15]  into  the  spectral  domain,  the  integration 
with  respect  to  y  computed  analytically  and  the  remaining 
single  integration  is  performed  numerically  on  a  digital 
computer.   Here  again  use  is  made  of  the  fact  that  the 
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transforms  of  the  scalar  potentials  are  known  after  one  has 
solved  for  the  propagation  constant  6  according  to  Section 
2.1.1.   The  time  average  power  flow  P     is  computed  as 
explained  in  Section  2.1.2.   Some  details  of  the  derivation 
leading  to  this  single  integral  are  shown  in  Appendix  D,  the 
method  is  the  same  as  outlined  in  Sections  2.1  through  2.3. 
This  perturbational  analysis  was  carried  out  for  an 
arbitrarily  chosen  set  of  parameters,  such  that  resonance 
occurred  at  4  GHz.   Figure  2.4/2  shows  the  result  of  this 
calculation  where  3„/$  represents  the  ferrite  phase  constant 

r 

normalized  by  the  pure  dielectric  phase  constant.   The 
magnetic  biasing  field  of  H   =  683.6  Oe  was  directed  into 
the  positive  x  direction.   The  material  constants  are 
4ttMc  =  2  300  Gauss,  AH  =  100  Oe  and  e   =20. 

This  result  in  Figure  2.4/2  can  be  only  qualitatively 
checked  with  results  for  similar  structures  for  which  a  theory 
has  been  verified  by  experiment.   Similar  to  the  propagation 
constant  for  a  ferrite  filled  parallel  plane  waveguide  with 
transverse  magnetization  [Ref.  2]  it  is  seen  that  there  are 
two  regions  of  propagation,  below  and  above  the.  resonance, 
frequency,  separated  by  a  region  of  high  attenuation.   It  is 
noted  that  this  attenuation  is  dependent  on  the  direction  of 
propagation.   Another  qualitative  comparison  can  be  made  here 
with  the  resonance  isolator  described  in  Reference  19.   Both 
waves  traveling  in  either  direction  are  attenuated  at 
resonance,  although  one  not  as  much  as  the  other.   According 
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to  Helszajn  [Ref.  19]  optimum  forward  attenuation  dees  not 

correspond  exactly  to  optimum  reverse  attenuation  for 

X   4   x   ,  since  the  normal  modes  are  in  general  elliptically 

yy   zz 

polarized.   Due  to  the  demagnetization  described  above  this 
is  the  case  here.   Furthermore  it  must  be  noted  that  for  the 
resonance  isolator  not  to  attenuate  one  of  the  propagating 
modes,  the  ferrite  slab  should  be  placed  close  to  the  center 
at  some  optimum  position.   But  here  the  physical  nature  of 
the  slot  line  determines  an  offcenter  position  of  the  ferrite 
if  the  conductor  plane  in  y  =  D  is  thought  of  as  a  line  of 
symmetry. 

2 .  5   SUMMARY 

A  new  method  for  a  frequency  dependent  analysis  of  open 
boundary  coplanar  transmission  lines  has  been  formulated. 
The  two  important  line  parameters,  phase  velocity  and 
characteristic  impedance,  are  derived  from  a  complete  hybrid 
mode  analysis  through  scalar  potential  functions  which 
satisfy  Helmholtz's  equation.   This  P.D.E.  is  solved  by  a 
Fourier  transform  technique  and  a  matching  of  boundary 
conditions  in  terms  of  the  electric  field  and  current 
density  components.   The  application  of  Galerkin's  method  in 
the  spectral  domain  and  an  approximation  of  the  electric 
field  or  current  density  by  basis  functions  permits  a  solution 
for  the  phase  velocity  to  any  desired  degree  of  accuracy. 
Through  the  application  of  Parseval's  theorem  to  the  integral 
of  the  complex  Poynting  vector  in  the  space  domain  the 
evaluation  of  time  average  power  flow  on  the  waveguide 
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structure  in  the  spectral  domain  is  feasible  and  allows  the 
computation  of  the  characteristic  impedance  of  the  line. 
The  often  cumbersome  inversion  of  a  Fourier  transform  is 
thus  avoided. 

The  only  other  existing  analysis  of  a  single  slot  line, 
Cohn_[Ref.  3],  has  been  numerically  verified  by  a  surprisingly 
close  agreement  between  these  two  theories  which  are  radically 
different.   The  present  theory  also  provides  new  and  more 
detailed  information  about  the  power  density  on  the  slot  line 
structure. 

The  formulation  of  the  problem  is  so  general  that  a 
simple  change  in  the  assumed  field  distribution  yielded  a 
complete  analysis  of  two  coupled  slot  lines.   The  data 
obtained  is  new  and  will  be  of  major  importance  for  the 
microwave  industry  in  the  design  of  slot  line  filters  and 
couplers.   Wherever  a  comparison  with  related  theories  was 
possible  excellent  agreement  was  found.   Another  modification 
of  the  existing  theory  provided  the  complete  analysis  of  the 
coplanar  strip  transmission  line.   To  some  extent  a  comparison 
with  experimental  data  was  possible  and  showed  a  very  good 
agreement. 

A  perturbational  spectral  domain  method  for  the  analysis 
of  the  propagation  constant  of  a  coplanar  transmission  line 
on  a  ferrite  substrate  was  developed  and  applied  to  the 
specific  example  of  a  slot.  line.   An  experimental  verification 
was  not  performed. 
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The  spectral  domain  technique  is  relatively  straight- 
forward in  concept,  but  extensive  algebraic  manipulation  is 
required  to  achieve  computational  efficiency.   However,  once 
the  boundary  value  problems  of  Helmholtz's  equation  is 
solved  in  terms  of  electric  fields  and  current  densities  at 
the  air-metal-dielectric  interface,  any  coplanar  waveguide 
structure  can  be  easily  analyzed. 
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3.   CLOSED  BOUNDARY  STRUCTURES 

Lumped  circuit  elements  are  no  longer  of  any  practical 
use  at  microwave  frequencies.   Losses  due  to  the  increased 
resistance  by  the  skin  effect  or  due  to  radiation  effectively 
terminates  the  use  of  such  elements.   A  common  alternative  is 
the  construction  of  circuit  elements  by  sections  of  trans- 
mission lines  of  which  the  quarter  wave  transformer  is  one 
of  the  most  important  examples.   Using  standard  transmission 
lines  such  as  coaxial  cables  or  rectangular  waveguide  these 
lines  are  mostly  short  circuited  at  one  end  to  keep  radiation 
losses  to  a  minimum.   The  short  circuit  can  be  assumed  in 
most  cases  to  be  ideal.   This  may  be  true  and  practical  for 
closed  boundary  structures  but  is  not  easy  to  implement  with 
an  open  boundary  structure  such  as  the  slot-line.   If  a  slot- 
line  is  simply  terminated  within  the  conducting  plane  by 
ending  the  slot  it  is  found  that  a  considerable  amount  of 
energy  is  stored  behind  this  short  circuit  due  to  surface 
currents  flowing  around  this  end.   This  indicates  that 
magnetic  energy  will  be  stored  in  this  region  so  that  the 
short  circuit  is  not  ideal  but  rather  appears  as  an  inductive 
reactance  at  a  reference  plane  coincident  with  the  end  of  the 
slot.   This  effect  was  first  observed  by  Mariani  and  Agrios 
[Ref.  23]  during  their  experimental  investigation  of  slot 
line  filters  and  later  experimentally  investigated  by  Knorr 
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and  Saenz  [Ref.  24]  who  also  provided  data  curves  of  this 
so  called  "end-effect"  for  design  purposes.   No  theory  has 
appeared  yet  which  describes  this  effect. 

Another  example  of  how  lumped  circuits  are  replaced  by 
alternative  circuits  at  microwave  frequencies  is  the 
resonating  circuit.   Resonating  structures  are  important  in 
the  design  of  filters  and  oscillators  and  have  been  constructed 
as  cavities  for  rectangular  or  circular  waveguides.   For  the 
slot-line  to  be  an  alternative  transmission  line  to  micro- 
strip  line  the  use  of  resonating  slots  becomes  evident  for 
previously  mentioned  applications.   Quarter  wave  coupled 
resonating  slots  and  end  coupled  resonating  slots  have  been 
presented  on  an  experimental  basis  by  Mariani  and  Agrios 
[Ref.  23] ,  but  again  no  theory  exists  to  provide  means  for 
the  design.   The  problem  here  is  how  to  determine  the 
dimensions  of  a  resonating  slot  such  that  resonance  at  some 
desired  frequency  occurs  which  is  often  the  lowest  possible 
frequency.   Using  the  wavelength  ratio  of  the  slot  line 
computed  in  Section  2  for  this  half  wavelength  resonator 
will  not  yield  this  information  because  of  the  above 
described  end-effect.   The  wavelength  ratio  was  derived  for 
a  transmission  line,  which  extends  infinitely  in  the 
z-direction  (see  Figure  1.1/2). 

An  analysis  of  these  kinds  of  problems  is  now  proposed 
in  this  section.   The  hybrid  mode  analysis  in  the  Fourier 
domain  is  used  again  as  in  Section  2  to  provide  a  complete 
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frequency  dependent  solution.   To  minimize  the  losses  due 
to  radiation  the  resonating  slot  section  is  placed  inside  a 
shielding  structure;  this  configuration  is  schematically 
shown  in  Figure  3.1/1. 

As  a  byproduct  of  this  analysis  the  above  cited  end 
effect  is  solved  theoretically.   Some  further  modifications 
in  the  derived  theory  of  the  resonating  slot  structure 
unveils  that  line  parameters  such  as  the  dispersion 
characteristic  and  the  characteristic  impedance  for  shielded 
coplanar  transmission  lines  can  be  obtained.   As  an  example 
results  for  a  single  shielded  slot  line  are  provided. 

An  important  result  of  Section  2  is  that  once  the  boundary 
value  problem  of  Helmholtz's  equation  is  solved  in  terms  of 
the  electric  fields  and  surface  currents  in  the  plane  between 
the  air-conductor-dielectric  interface  any  similar  structure 
can  be  easily  analyzed  by  minor  modifications.   This  result 
holds  also  for  the  closed  boundary  structures  in  this 
section  and  will  be  demonstrated  with  the  example  of  end 
coupled  resonating  slot-lines.   This  structure  can  be 
described  by  a  circuit  model,  where  the  elements  of  this 
model  can  be  solved  for.   The  result  is  a  solution  to  an 
inductive  coupling  between  two  serial  slot-lines. 

Throughout  the  following  analyses  in  this  section  it  is 
assumed  that  the  two  rectangular  waveguides  formed  by  the 
shielding  structures  have  dimensions  such  that  rectangular 
waveguide  modes  are  kept  below  cutoff.   The  cutoff  condition 
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Figure  3.1/1   Shielded  Slot  Line  Resonator 
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for  the  waveguide  formed  by  region  1  in  Figure  3.1/1  is 
well  known.   Besides,  if  dimensions  h-,  and  h2  are  kept  such 
that  only  a  TE,0  limit  mode  (E  field  in  y  direction)  could 
be  excited,  the  lower  waveguide  iiormed  by  regions  2  and  3 
will  have  the  lower  cutoff  frequency  due  to  the  dielectic 
loading  and  has  to  be  considered  only.   A  cutoff  condition 
was  derived  by  Harrington  [Ref.  21]  for  instance  for  a 
partially  dielectric-filled  rectangular  waveguide.   A 
solution  is  required  for  a  transcendental  equation  given 
according  to  Reference  21  for  the  lowest  possible  mode  as 


"ST  tan<kx2D)  ■  "  TT  tan(kx3h2>  <3-1' 


where   k  -  =  to  e2y2  ~  ^T~^  (3-2a) 


kx3  "  "\v3.r    <£>2  (3-2b) 


The  subscripts  2  and  3  refer  to  the  regions  2  and  3  in 
Figure  3.1/1.   This  equation  (3-1)  is  given  here  only  for 
completeness  and  is  solved  numerically  within  the  computer 
program  described  in  Section  5.2. 

3.1   SINGLE  RESONANT  SLOT  LINE 

The  schematic  of  the  resonating  structure  to  be  analyzed 
is  shown  in  Figure  3.1/1  where  it  is  assumed  that  the 
shielding  structure  extends  for  z  such  that   -°°  <  z  <  °°. 
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The  conducting  material  on  top  of  the  dielectric  and  the 
surrounding  shield  is  assumed  to  be  perfect  and  the  dielectric 
material  is  lossless  and  isotropic. 

Similar  to  the  beginning  in  Section  2.1  a  full  mode 
analysis  will  require  two  linearly  independent  scalar 
potential  functions.   The  z  directed  electric  and  magnetic 
field  components  in  terms  of  these  potentials  may  then  be 
written  according  to  References  13  or  21  as 


Ezi(x,y,z)  =  (-^-2  +  ki2)4)i(x,y,z)  (3-3) 

8z 

2 
Hzi(x,y,z)  =  (l-y  +  k±2)^i(x/y/z)  (3-4) 

8z 


1=1,2,3  defining  the  spatial  region  according  to  Figure 
3.1/1 


k32  =  kl2  =  k02  =  w2y0e0    and   k22  =  Vk02 


The  transverse  field  components  of  the  hybrid  mode  can  then 
be  derived  through  an  application  of  Maxwell's  curl  equations 
to  (3-3)  and  (3-4) . 
Then 

E  .  =  •= k—  ~   JwVn  ^ (3-5a) 

xi    3x  oz    J   0  3y 

82cf).  8\Jj. 

Eyi  =  37^7  +  j^o   s3T  (3"5b) 
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9*i    824>i 

H   .   =  JU)G.  vr +  ■* *r-~  (3~6a) 

xi   J   -i  3y    8x  9z 

2 

H  .  =  -jcoe.  r— =■  +  t J—  (3-6b) 

yi     J   1  dx  dy    3z 


The  scalar  potential  functions  satisfy  Helmholtz's  equation 
in  the  three  spatial  regions  [Ref.  21],  thus 


V2^  +  ki24>i  =  0  (3-7) 


V2^  +  ki2ipi  =  0  (3-8) 


2 

where  V   denotes  the  three  dimensional  Laplacian  operator 

in  x,y,z  coordinates.   Two  Fourier  transforms  are  now  used 
to  transform  these  three-dimensional  partial  differential 
equations  into  ordinary  differential  equations  which  can  be 
readily  solved  in  terms  of  the  boundary  conditions. 
Define  a  two  dimensional  Fourier  transform  as 


00   a  j(a  x+£z) 

G(a  ,y,£)  =     f     I        g(x,y,z)  e   n      dxdz    (3-9) 
_oo  _a 


Note  that  the  transform  with  respect  to  x  is  a  finite 
transform  [Ref.  25] ,  also  widely  known  as  a  complex  Fourier 
series  coefficient. 

Some  physical  reasoning  must  be  used  here  to  determine 
the  discrete  transform  variable  a  .   For  the  lowest  order 
mode  in  this  resonant  slot  structure  the  fields  of  the 
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standing  waves  inside  the  slot  in  the  z  direction  are  similar 
to  the  z  directed  fields  of  a  traveling  wave  shown  in 
Figure  2.4/1.   Thus,  the  z  directed  magnetic  field  component 
H  will  be  an  even  function  with  respect  to  x  and  from 
Section  2.1.1  and  the  second  order  approximation  it  is  known 
that  the  z  directed  electric  field  component  E   is  an  odd 
function  with  respect  to  x.   A  distinction  in  equation  (3-9) 
in  sine-  and  cosine-transforms  aids  in  the  following 
discussion.   It  is  easily  seen  that  if 


IT  r 

F(n)  =   /   f(x)  \°°*    (nx)  dx 


then 


/  f"  (x)cos  (nx)dx  =  ~nzF(n)  +  (-l)nf'(Tr)  -  f'(-Tr) 


-TT 


and 


/  g"  (x)sin(nx)dx  =  -n2G(n)  +  n[g(-ir)  -  (-3fg(ir)] 

-TT 

So  these  transforms  are  expedient  for  problems  in  which 

either  the  function  to  be  transformed  or  its  derivative 

is  prescribed  at  the  endpoints  of  the  interval.   But  note 

that  for  the  above  problem  at  the  endpoints  x  =  -a  and  x  =  a 

the  tangential  electric  fields  must  be  zero  due  to  the 

conducting  shield,  hence  E  .  and  E  . .  or  as  seen  from 
3  zi      yi 

equations  (3-3)  and  (3-5)  <f> .  and  ■* vanish  at  these 

1        o  X 

boundaries  since  E  .  -  <b.  and  H  .  ~  ib.     .   The  transform 

zi   Yi        zi    ri 
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variable  a   in  (3-9)  is  thus  determined  to  be 


aR  =  ~L    for  n  =  0,  ±1,  ±2,  ... 


With  these  considerations  the  two  dimensional  transform  of 
equations  (3-3)  and  (3-4)  via  (3-9)  can  be  accomplished 
to  obtain 


Ezi(cVy'C)  =  (ki2  "  ?2)  *i(an'y'^  (3"10) 


Hzi(ctn/y,£)  =  (k±2  -  £2)  Vcxn,y,S)  ^-11) 


Similarly  equations  (3-5a)  and  (3-6a)  are  transformed  into 


Exi  -  "an^i  "  ^v0   W-  (3-12) 

Hxi  =  i"£i  W"  -   Vi  (3'13) 


The  y  directed  components  are  of  no  concern  in  the  development 

Define 

2     9  2  ? 

an   +  r  -  k±  =   y±Z  (3-14) 

then  the  two  Helmholtz's  equations  in  (3-7)  and  (3-8) 
are  transformed  into 


2 

(•^-j  -  Y±  )  *i(anfy#0  -  0  (3-15) 

3y 

2 

(•^  -  Y,2)  V.   (an#y,C)  -  0  (3-16) 

9yz     a.     in 
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It  is  observed,  that  the  boundary  conditions  for  the  tangen- 
tial electric  fields  of  y  =  D  +  h.  and  at  y  =  -h„  require 
E   =  0  and  E   =  0.  so  that  certain  solutions  to  (3-15)  and 

X  Z      ' 

(3-16)  are  eliminated.  These  conditions  are  independent  of 
x  and  z  and  may  be  applied  in  the  transform  domain  as  well. 
The  solutions  to  (3-15)  and  (3-16)  are  then 


VVY'^  =  Ae(cxn,£)  sinh  Y1(D+h1-y) 


*2(an,y,£)  =  Be(an,0  sinh  y2y  +  Ce(an,£)  cosh  y2Y  (3-17) 


&3(a  ,y,£)  =  De(an,£)  sinh  y^h^+y) 


1- 

yi(an/y,C)  =  A  (un,0    cosh  y^D-Hi^y) 


H/2(an,y/C)  =  Bh(an/^)  sinh  Y2y  +  Ch(an,0  cosh  y 2y  (3-18) 


Wy'5)  =  E)h^an'^)  cosh  Yl(h2+Y) 


It  is  to  be  noted  at  this  point  that  n  may  assume  the 

2 

value  zero  and  -»  <  £  <  «>.   y.   may  therefore  be  less  than 

2      2 
zero  for  certain  values  of  £.   Furthermore,  since  k«   >  k, 

there  will  be  three  sets  of  solutions  in  equations  (3-17) 

2         2 
and  (3-18)  for  n  =  0,  namely  y-.   <  0/  Y?   <  0  where  trigo- 

2 

nometric  functions  replace  all  hyperbolic  functions;  y,   >  0, 

2 

y9   <  0  where  trigonometric  functions  replace  the  hyperbolic 

2         2 
functions  in  spatial  region  2  and  y.   >  0,  y.   >  0  for  which 
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the  solutions  are  as  given  above.   It  can  furthermore  be 


seen  that  these  three  types  of  solutions  will  exist  only  for 

X 


2 

n  =  0.   Assume  n  =  1  and  £   assumes  its  lowest  value  of   £=  0 . 


2  7T 

For  Y.   <  0  then  —  <  k.  or  -^-  <  a.   However,  to  avoid  rectan- 
'  l  a     l     2 

gular  waveguide  modes  it  is  assumed  that  y  >  2a,  thus 

2 

y.       >  0  for  n  ±1,  ±2,  ...   at  all  frequencies  of  interest. 

According  to  the  method  of  solution  outlined  in  Section  2  the 
next  step  is  the  application  of  the  continuity  or  boundary 
conditions  at  the  interfaces  between  the  three  spatial  regions 
These  conditions  are  given  in  equations  (2-14)  and  (2-15) 
and  apply  for  this  closed  boundary  structure  as  well  with  a 
minor  modification.   The  last  four  equations  of  (2-15)  must 
now  include  a  general  dependence  on  the  variable  z  for  the 
electric  field  and  current  density  components.  A  z  dependence 

Y  Z 

of  e    is  sufficient  for  a  transmission  line  only.   For 
instance,  the  current  density  component  in  the  x  direction 
would  now  be  of  the  form  j  (x,z),  the  other  components  are 
changed  similarly  in  (2-15) .   Denote  the  transforms  of  these 
components  at  y  =  D  as 


E  (a  ,£)  =  ¥{E    (x,z)}         E  (a  ,£)  =  V{E    (x,z)} 
x  n         x  z   n  ^       z   ' 


J  (a  ,K)    =  F{j  (x,z)>         Van'^  =  f(j(x,z)}     (3-19 


If  these  components  and  the  transforms  of  (2-14)  and  (2-15) 
are  used  with  the  substitution  of  the  scalar  potential 
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transforms  through  the  solutions  to  Helmholtz's  equations 

(3-17)  and  (3-18)  it  is  possible  to  eliminate  again  the 

e  h 

unknown  coefficients  A   through  D  .   Details  of  this  exten- 
sive algebraic  manipulation  are  shown  in  Appendix  E,  the 
result  of  which  is 


M, (a  ,£)     M0 (a  ,£) 
In        2   n 


M, (a  ,£)     M„ (a  ,£) 

3   n's       4   n 


E  (a  ,0 
x  n'  *' 


E  (a  ,0 
z   n  ^ 


J  (a  ,0 
x  n' s 


J  (a  ,£) 
z  n 


(3-20) 


To  apply  here  again  the  method  of  moments  by  Galerkin's 
method  [Ref.  16]  in  the  spectral  domain  the  following 
definition  of  the  scalar  product  is  used  [Ref.  15] 


<f(an,0  ,g(a  ,£)   =   S     /  f  (a  ,  U  g*  ( aR  ,£)  d£ 


(3-21) 


For  the  application  of  Galerkin's  method  the  suitable  set  of 
weighting  functions  is  (5  (a  ,£)/E  (a  ,£)).   With  this  set 
and  through  equations  (3-21) ,  (3-20)  becomes 


<M1(an,C)Ex(an,C),Ex(an,C)>  +  <M2  (an'C)  Ez  (an'?)  ,Ex(otn'C)  >  =  ° 


(3-22) 


<M3(VOEx(ahrOfEB(ahf5)>   +    <M4  <V  V  *z  (a^O  fE,  <V«>  >   =    0 
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The  transforms  of  the  field  components  are  not  known  but 

can  be  expanded  in  terms  of  basis  functions  such  that 

00 

E  (<*.£)  -   E   a   E  .  (a  ,U  (3-23a) 


E  (a  ,£)  =   Z   b.  E  .  (a  ,£)  (3-23b) 

z  n  ^     ._,   l   zi   n  ^ 


Equation  (3-22)  can  then  be  written  explicitly  as 


00  00  oo 

Z   a.      Z  /  M,  (o    ,C)E    ,  (a    ,£)E       (a    ,£)    d£ 

i=l    1   n=_co   _oo         1      ^  xi      n's      xm      n's 

00  CO  CO 

+      Z   b.       Z  /  M.(a    ,£)E    .(a    ,?)E*    (a    ,£)    d£   =    0 

i=l   1   n=_oo   >oo        2      n  zl      n  xm     n 


m  =    1,2,3,. . .  (3-24) 


4ai   JL  J  M3(V«>Exi<Ve>C(V5>    d^ 

00  00  00 

+      E    b.       Z  /   M„(a    ,£)E    .(a    ,£)    E*    (a    ,£)    d£    =    0 

i=l    x   n=_co   _oo        4V    n's/     ziv    n's/       zm     n's/ 


m  —   1,2,3,... 

Equation  (3-24)  is  still  an  exact  formulation  of  the  resonant 
problem  to  which  a  solution  could  be  found  in  the  following 
fashion.   The  coefficients  M,  to  M.  are  besides  other  factors 
functions  of  the  operating  frequency,  while  the  geometry  of 
the  structure  and  therefore  the  slot  length  enters  through  the 
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basis  functions  into  the  calculations.   Given  a  desired  fre- 
quency of  operation  the  slot  length  L  is  changed  until  the 
determinant  of  the  coefficient  matrix  of  (3-24)  is  equal  to 
zero.   This  still  is  a  rather  theoretical  way  of  solving  this 
problem  since  in  (3-24)  two  infinite  sets  of  basis  functions 
are  implied.   Any  solution  will  involve  the  numerical  evalua- 
tions of  the  various  coefficients,  the  integrations  and  other 
algebraic  manipulations  so  that  an  approximate  solution  seems 
desirable.   In  Section  2  it  was  found  for  the  open  boundary 
structure  that  good  agreement  between  this  theory  and  others 
or  experiments  could  be  obtained  if  a  good  choice  for  a 
one  term  expansion  of  the  electric  field  components  was  made. 
If  the  slot  is  relatively  narrow  the  experience  with  the 
different  approximations  in  Section  2.1.1  suggests  the  z 
directed  electric  field  component  at  y  =  D  may  be  neglected 
while  the  x  directed  field  component  at  y  =  D  may  be  approxi- 
mated such  that 


Ex(x,D,z)  =  2x(x,z)  z    exl(x,z)  =  e^  (x)  -e^p  (z)     (3-25) 


where  1  i  ■    W 

x   < 


\f 


,W.  2     T 
(^r)   -  X 


2 


e.(i}  (x)  =  (3-26a) 


'xl 


elsewhere 


,  2T  cos  t  izl  K  * 

e.;f;  (z)  =  (3-26b) 


'xl 


elsewhere 
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The  Fourier  transform  of  (3-26a)  was  already  given  in 
equation  (2-2 3a)  but  is  presented  here  for  completeness  for 
the  finite  transform.   The  transform  of  (3-2 6b)  can  be  found 
analytically. 

The  transforms  of  (3-26)  are  then 


/ ,  >  a  W 

Exl  (an>  "  *  V-lH  (3"27a) 


/9v         TT  COS  (^-) 

Z\V  (g)  =  —  (3-27b) 

2        2 
*      -    (£L)Z 


where  B_  denotes  the  zero  order  Bessel  function  of  the  first 
kind. 

It  is  now  possible  to  obtain  a  solution  for  the  resonant 
slot  length  by  varying  L  such  that 


I  f      Mx(an,0  |Exi)(V'Exl)(0  I*  d?  =  °    (3"28) 

j}=— 00   —00 


This  root  finding  can  be  accomplished  now  by  a  numerical 
solution  on  a  digital  computer,  the  details  of  which  are 
given  in  Section  5. 

A  representative  graph  of  a  result  of  this  computation 
is  shown  in  Figure  3.1/2  where  the  resonant  length  versus 
frequency  for  the  lowest  resonant  mode  is  displayed.   All 
quantities  are  scaled  by  the  thickness  D  of  the  dielectric 
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Figure  3.1/2 


Resonant  Slot  Length  vs.  Normalized  Frequency 
for  a  Halfwave  Slot  Line  Resonator  with 


a/D  =  5.5, 


h1/D  =  h2/D 


8.0  and  er  =  20 
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substrate  which  makes  these  results  more  versatile  for  any 
applications.   These  results  are  new  and  cannot  be  compared 
with  any  values  resulting  from  ether  theories.   However, 
values  for  the  resonant  length  by  this  method  are  used  to 
compute,  theoretically,  the  end  effect  of  a  shorted  slot 
line.   Very  good  agreement  with  experiments  is  found  which 
substantiates  the  validity  of  the  theory  and  suffices  as 
an  experimental  verification. 

3.2   DISPERSION  CHARACTERISTIC  AND  CHARACTERISTIC  IMPEDANCE 

After  the  results  for  the  resonating  slots  were  obtained 
a  thorough  review  of  the  involved  formulas  led  to  the 
following  investigation  which  analyzes  the  line  parameters 
for  a  shielded  slot  line.   These  parameters  are  of  importance 
for  mainly  two  reasons. 

The  first  reason  is  that  for  a  theoretical  analysis  of 
the  end  effect  of  a  single  shorted  slot-line  the  propagation 
constant  of  this  line  must  be  known.   But  so  far  only  the 
open  boundary  slot  line  has  been  analyzed.   An  approximation 
could  certainly  be  made  in  assuming  that  the  boundary  walls 
are  far  enough  removed  so  that  their  effect  on  the  slot  line 
is  negligible  allowing  the  calculation  of  the  end  effect  with 
the  open  boundary  propagation  constant.   But  for  reasons  of 
completeness  line  parameters  for  the  shielded  line  should  be 
involved  to  provide  for  studies  on  the  effect  of  shields  on 
the  line  parameters. 
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This  effect  now  provides  the  second  reason.   If  the 
slot-line  is  to  be  used  in  integrated  microwave  circuits 
physical  dimensions  will  in  general  be  small.   Cohn  noted 
in  Reference  3  that  the  permittivity  of  the  dielectric 
substrate  must  be  sufficiently  large  to  confine  the  wave  to 
the  slot  and  to  minimize  radiation  losses.   However,  from 
the  graph  of  the  power  density  in  the  neighborhood  of  the 
slot  shown  in  Figure  2.1/4b  it  is  seen  that  even  for  e   =  16 
a  substantial  amount  of  power  is  found  for  distances  of  some 
multiples  of  slot  width  away  from  the  center.   Hence,  it  can 
be  anticipated  that  a  shielding  structure  around  the  slot 
line  seems  desirable  in  some  applications  to  comply  with  the 
miniaturization  of  microwave  circuits.   Since  no  analysis 
of  a  shielded  slot-line  has  appeared  yet  in  the  literature 
the  following  study  is  another  contribution  towards  a 
complete  investigation  of  coplanar  transmission  lines.   The 
method  used  follows  exactly  the  one  described  in  Sections 
2.1.1  and  2.1.2. 

A  careful  comparison  of  the  equations  used  in  Section  3.1 
for  the  resonant  slot  and  those  used  in  Section  2.1  for  the 
open  boundary  slot  line  reveals  that  only  slight  modifica- 
tions are  necessary  to  adapt  the  theory  in  Section  3.1  for 
an  analysis  of  shielded  coplanar  transmission  lines,  shown 
by  the  example  of  a  single  slot-line. 

If  the  general  scalar  potential  function  <j>(x,y,z) 

•  +YZ 

in  equation  (3-3)  is  replaced  by  c|>(x,y)e    which  accounts 
for  a  traveling  wave  in  the  negative  z  direction  it  is  found 
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that  equations  (2-1)  to  (2-8)  are  equivalent  to  (3-3)  to 
(3-8)  if  the  partial  differentiation  with  respect  to  z  is 
accomplished,  provided  a  similar  replacement  for  tj>(x,y,z) 
of  course.   Now  define  the  finite  Fourier  transform 


a        3  ax 
Fn{f(x,y)  }  =  $(an-,y)  =   /   <J>(x,y)e     dx       (3-29) 

-a 


where  a   is  determined  by  the  same  considerations  stated  in 
Section  3.1. 

The  Helmholtz's  equations  for  the  scalar  potential 
functions  are  then 


<7T"  Yi2)  *i(cVy)  "    °  <3"30> 


{^~2  "  Yi2)  Vvy)  =  °  (3"31) 

9y 


1=1,2,3 

2      2    2     2 
where  for  the  lossless  case   y  =  j$   and   y.   =  a   +  3  -k. 

Hence  equations  (2-10)  and  (2-11)  are  obtained  from  (3-15) 

and  (3-16)  if  £  is  replaced  by  3.   Similarly  equations  (3-10) 

through  (3-13)  hold  for  the  transmission  line  problem  if  this 

replacement  is  carried  out  there.   But  then  it  is  observed 

that  an  extensive  algebraic  manipulation  towards  a  solution 

like  that  of  Section  2.1  can  be  avoided  since  all  equations 

derived  in  Section  3.1  up  to  equation  (3-20)  and  all  of 

Appendix  E  hold. 
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The  final  part  of  the  formulation  is  then  obvious.   To 
solve  for  B  Galerkin's  method  is  applied  again,  which 
requires  the  scalar  product  to  be  defined  as 


<f (an) ,g(an)>  =   Z    f (aj -g* (aR)  .  (3-32) 


As  in  Section  2.1,  the  electric  field  components  at  y  =  D 
are  approximated  by  basis  functions  and  used  in  connection 
with  the  above  scalar  product  definition  on  equation  (3-20) 
to  obtain  the  set  of  simultaneous  equations.   The  zero  of 
the  determinant  of  the  associated  coefficient  matrix  will 
provide  the  desired  propagation  constant  3.   The  details 
are  by  now  familiar  and  given  in  Section  2.1.1. 

A  first  order  approximation  was  used  to  obtain  a  solution 
where  E  (x,D)  s  0   and  E  (x,D)  =  E    (x)   was  given  by 
equation  (2-22a) .   A  representative  result  and  comparison 
with  an  open  boundary  slot  line  is  given  at  the  end  of  this 
section. 

Inspired  by  the  experience  that  only  minor  modifications 
were  necessary  to  adapt  equations  from  Section  3.1  for  a 
solution  of  the  shielded  slot  line,  the  following  investiga- 
tion is  determined  to  evaluate  the  second  line  parameter,  the 
characteristic  impedance,  whose  definition  is  given  in 
equation  (2-2  4) .   The  time  average  power  flow  on  the  shielded 
slot  line  can  be  computed  according  to  equation  (2-26)  as 


10  0 


Pavg  "  I  Re  J1+°  _{      t-U0B[(||)2+  (f)2l  +  -B[(f|)2  +  (ff)2] 

+  (62+k2)  (||  ||  -  ||  f|)  }   dx  dy  .   (3-33) 


This  formulation  applies  to  the  spatial  domain  and  cannot  be 
computed  since  the  scalar  potential  functions  are  not  known 
in  this  domain.   Through  the  use  of  Parseval's  relation 
[Ref.  15]  it  is  possible,  however,  to  transform  this 
expression  into  the  spectral  domain  where  the  scalar  potentials 
are  known  after  the  propagation  constant  3  is  solved  for. 

Due  to  the  finite  transform  in  equation  (3-29)  the  integra- 
tion with  respect  to  x  is  transformed  via  Parseval's  relation 
into  an  infinite  summation,  thus  (3-33)  becomes 


oo  h   +D 

Pavg  =  ^Re      l  „•'       f-ar/WM2    "    %2Soniom 

3  n=-°°   -h2 


|9M'.2  ,3$, 2 


3anB    (^F"        9Y   *   } 


9  Pi  W  aw* 

+    jankz(<D*  |i  +  |i-  y)}      dy.     (3-34) 


According  to  the  general  procedure  outlined  in  Section 
2.1.2  each  coefficient  of  this  series  can  be  evaluated.   The 
transforms  of  the  scalar  potential  functions  are  given  by 
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equations  (3-17)  and  (3-18)  where  again  the  transform 
variable  £  is  replaced  by  the  propagation  constant  3.   The 
simple  functional  dependence  on  the  variable  y  allows  one 
to  perform  the  integration  in  equation  (3-34)  analytically. 
Note  that  each  coefficient  of  the  series  in  (3-34)  is  actually 
the  sum  of  three  further  coefficients  which  present  the  power 
carried  in  each  of  the  three  spatial  regions.   The  details 
about  these  are  shown  in  Appendix  F. 

The  final  part  in  the  evaluation  of  the  characteristic 
impedance  is  shown  in  Section  2.1.2  and  involves  only  the 
computation  of  the  slot  voltage  to  obtain  the  characteristic 
impedance  according  to  equation  (2-2  4) .   For  programming 
purposes  it  should  be  noted  that  the  computation  for  the 
shielded  line  is  easier  to  perform  since  a  numerical  integra- 
tion is  avoided  and  replaced  by  a  summation.   The  coefficients 
in  the  infinite  series  of  equation  (3-34)  were  found  to  decay 
rapidly  so  that  a  finite  approximation  yields  good  results. 

A  representative  set  of  curves  from  this  calculation  is 

displayed  in  Figure  3.2/la-b  where  the  wavelength  ratio  \x /\ 

and  the  characteristic  impedance  ZQ  of  the  shielded  slot-line 

is  compared  with  the  line  parameters  of  an  equivalent  open 

boundary  slot  line  for  two  different  slot  widths.   The  cutoff 

frequency  for  the  lowest  order  waveguide  mode  was  computed 

according  to  equation  (3-1)  and  is  in  this  example  in 

normalized  form  D/X   =  0.0421. 

c 

Several  facts  can  be  learned  from  this  comparison. 
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a)  Dispersion  Characteristic  and  b)  Characteristic 

Impedance  of  a  Shielded  Slot  Lino  vs.  Normalized 
Frequency  with  a/D  =  5.5,  h,/D  :  h^/D    8.0  and 


=  20 


h1/D  =  h2/D 
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It  is  observed  that  the  difference  between  the  parameters 
of  open  and  shielded  structures  differ  more  for  the  lower 
frequency  range  than  for  the  higher  one,  where  the  parameters 
of  both  lines  merge  together.   This  is  no  surprise  if  the 
following  computation  is  performed.   According  to  Cohn 
[Ref.  3]  the  ratio  of  the  voltage  V(r)  along  a  semicircular 
path  of  constant  radius  r  in  region  1  divided  by  the  voltage 
V  directly  across  the  slot  can  be  given  as 


V(r)  =  ;kr|H{l)(kr)|  (3-35) 


V     2   c" l"l   v  c 


where  h|    is  the  Hankel  function  of  first  kind  and  first 
order  and 


V  =  '2"is   V(rr)2  -  1  <3-36' 


The  Hankel  function  h!   (jx)  asymptotically  approaches 
e   /  x   for  x  -*■  °°   so  that  an  approximate  result  for  the 
voltage  ratio  in  equation  (3-35)  may  be  easily  obtained. 
In  the  above  stated  example  the  waveguide  walls  are  at 
r/D  =  5.5  in  the  x  direction.   Then  for  W/D  =  0.2  a  voltage 
ratio  of  0.34  is  computed  for  D/A  =  0.02  and  of  0.0  5  for 
D/X  =  0.04.   Similarly,  the  voltage  ratio  for  W/D  =  1.0  is 
0.39  and  0.07  for  the  normalized  frequencies  0.02  and  0.04 
respectively  at  the  position  of  the  sidewalls.   One  may  now 
conjecture  that  at  the  lower  frequencies  a  considerable 
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amount  of  surface  current  flows  on  the  sidewalls  of  the 
shields  at  | x/D |  =  5.5.   This  in  turn  means  that  due  to 
these  surface  currents  more  power  flow  is  distributed  in 
the  air  regions  1  and  3  than  is  the  case  for  the  open 
boundary  line.   This  conjecture  is  verified  and  can  be  seen 
in  Figure  3.2/1  where  the  shielded  slot  line  shows  a  larger 
wavelength  ratio  and  a  larger  impedance  for  the  lower 
frequency  range  due  to  this  larger  amount  of  power  flow  in 
air.   The  above  calculation  shows  furthermore  that  the  shields 
will  effect  a  wider  slot-line  slightly  more  than  a  narrow  one, 
due  to  the  closer  confinement  to  the  slot  region  of  the  latter. 
Again,  this  is  verified  by  investigating  Figure  3.2/1. 

As  the  frequency  increases,  the  shields  have  less  effect 
and  the  line  parameters  approach  those  for  the  open  boundary 
line.   However,  the  upper  limit  on  the  frequency  of  operation 
is  the  onset  of  a  rectangular  TE, 0  mode  in  the  lower  waveguide 
formed  by  regions  2  and  3. 

3.3   END  EFFECT  OF  SHORTED  SLOT  LINE 

With  the  results  of  the  preceding  sections  the  end  effect 
of  a  short  circuited  slot-line  can  now  be  computed  analytically 
A  convenient  description  of  this  effect  is  based  upon  trans- 
mission line  theory  and  can  be  performed  in  terms  of  a 
normalized  reactance.   This  method  was  used  by  Knorr  and 
Saenz  [Ref.  24]  in  their  experimental  investigation  of  this 
effect  and  is  accepted  here  to  provide  for  a  comparison 
between  theory  and  experiment. 


Consider  a  slot-line  which  is  shorted  at  one  end  by 
ending  the  slot  in  the  conductor  plane.   The  slot  voltage 
would  not  only  be  zero  at  this  short  but  also  at  multiples 
of  A'/2  along  the  line  if  the  short  were  perfect,  thus 
forming  a  standing  wave  pattern.   Any  section  on  this  line 
between  two  voltage  nulls  could  be  considered  then  as  being 
a  half  wave  resonator  described  and  analyzed  in  Section  3.1. 
If  the  half  wavelength  A'/2  of  a  wave  on  this  line  is  compared 
with  the  resonating  length  L  of  the  resonator  for  comparable 
physical  dimensions  at  the  same  frequency,  however,  it  is 
found  that  A • /2  >  L.   This  could  be  expected  since  merely 
ending  the  slot  in  the  plane  of  the  conductor  does  not 
provide  a  perfect  short  circuit.   Surface  currents  on  the 
conductor  beyond  the  extension  of  the  slot-line  increase  the 
electrical  length  of  the  line  so  that  the  virtual  short 
circuit  appears  at  a  distance  A  * /4  -L/2  behind  the  physical 
short  at  the  end  of  the  line. 

The  reactance  caused  by  this  effect  is  then  computed  as 
follows.   The  input  impedance  of  a  short  circuited  trans- 
mission line  is 

Zi  =  j  ZQ  tan  3£  (3-37) 

where  3  =  2tt/X'  and  I   =  A'/4  -  L/2.   The  normalized  reactance 
can  then  be  given  as 


X,  =  tan  [it  (1/2  -  L/X')]  (3-38) 

s 
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where  the  necessary  data  is  provided  by  the  analytic/numerical 
methods  developed  in  Sections  3.1  and  3.2.   Results  of  this 
calculation  are  shown  in  Figure  3.3/la-b  and  are  compared 
with  the  experimental  data  from  Reference  24.   Although 
actually  the  theoretical  result  for  a  shielded  line  is 
compared  with  the  experiment  from  an  open  boundary  line, 
excellent  agreement  can  be  seen.   It  should  be  noted,  however, 
that  the  data  presented  was  obtained  by  assuming  that  the 
z-directed  electric  field  component  could  be  neglected.   It 
was  found  in  Section  2.1.1  that  this  component  becomes 
significant  for  normalized  frequencies  above  0.03.   This  fact 
may  explain  why  there  is  a  slight  discrepancy  between  theory 
and  experiment  in  Figure  3.3/lb  for  W/D  =  0.541.   A  solution 
which  accounts  for  the  z-directed  field  component  can  be 
easily  established  since  all  necessary  computations  are 
given  in  Appendix  E.   This  work  was  not  performed  here, 
since  the  amount  of  computation  time  on  the  digital  computer 
increases  substantially.   The  good  agreement  between  theory 
and  experiment  adequately  verifies  the  method. 

The  resulting  data  from  this  section  are  the  basis  for 
the  theoretical  design  of  microstrip  to  slot-line  transitions 
which  are  based  so  far  on  the  experimental  data  from  Knorr 
and  Saenz  [Ref.  24].   Furthermore  these  results  are  used  in 
the  analysis  of  the  following  section. 
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W/D=  0.562 


W/D  =  0.2:2! 
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0.03         0.04 


W/D=  0.541 
W/D  =  0.207 


b) 

0  0.01  0.02         0.03  0.04 

NORMALIZED   FREQUENCY     D/X 

Figure    3.3/1      Short-Circuit  Reactance   Referred   to   a  Plane 
Coincident  with   End   of   Slot   Line    for 
a)     £r   =    12,    b)     cr   =    20 

with    a/D  =5.5    and  h1/D=h2/D   =8.0 
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3.4   END  COUPLED  SLOT-LINES 

An  important  results  of  Section  2  is  that  once  the 
Helmholtz's  equations  in  the  various  spatial  regions  are 
solved  by  the  noted  boundary  matching  technique  in  terms  of 
the  surface  currents  and  electric  field  components  in  the 
plane  y  a  D  any  configuration  of  conductors  in  this  plane  may 
be  solved  for  the  line  or  resonant  parameters  by  approxi- 
mating either  one  of  these  field  quantities.   While  in 
Section  2  this  feature  was  demonstrated  with  the  analysis 
of  coupled  slot  lines  and  coplanar  strip  lines,  the  example 
presented  here  shows  how  still  another  slot-line  configuration 
may  be  analyzed  by  a  modification  of  the  theory  derived  in 
Section  3.1. 

The  coupling  effect  between  transmission  lines  and 
resonators  is  used  in  the  design  of  couplers  and  filters  at 
microwave  frequencies.   This  fact  was  mentioned  before  and 
is  shown  theoretically  and  experimentally  in  Section  4  for 
parallel  coupling  between  two  slot-lines.   An  analysis  of 
parallel  coupled  slot  line  resonators  for  instance  would 
basically  follow  the  same  procedure  outlined  in  Section  2.2. 
Another  configuration  which  has  not  yet  been  studied  is 
investigated  here:  end  coupled  slot-lines.   This  end  coupling 
effect  between  the  two  lines  can  be  described  easily  in  terms 
of  a  mutual  inductance  since  it  was  found  in  the  previous 
section  that  a  simple  short  circuited  slot  line  could  be 
represented  by  an  inductance  at  the  end  of  the  line. 
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The  following  analysis  studies  this  coupling  effect 
between  two  end  coupled  slot  resonators  in  terms  of  odd 
and  even  mode  resonant  lengths.   A  simple  application  of 
transmission  line  theory  unveils  the  coupling  effect  in  terms 
of  network  parameters  such  as  end  inductances  and  mutual 
inductance.   A  physical  layout  of  two  end  coupled  slot 
resonators  is  shown  in  Figure  3.4/la  where  the  shielding 
structure  is  not  displayed  since  the  same  configuration  as 
in  Figure  3.1/1  applies.   A  circuit  model  of  the  two  slots 
with  their  mutual  inductance  can  then  be  thought  of  as  given 
in  Figure  3.4/2.   The  analysis  is  then  as  follows. 

The  impedance  of  an  ideal  transmission  at  any  point  on 
this  line  in  terms  of  the  load  impedance  Z_ ,  the  characteristic 
impedance  ZQ  of  the  line  and  the  electrical  length  3£  is 
given  as 


ZL  +  jZQ  tan  6£ 
Z  =  Z0  ZQ  +  jZL  tan  fU  (3-39) 


where  3  is  the  propagation  constant  and  i   is  the  length  of 

the  line  from  the  load  impedance  to  the  point  of  interest. 

In  this  problem  the  load  impedance  ZT  is  Z^  =  jX_,  the  short 

JLi       s       s 

circuit  reactance  due  to  the  end  effect  which  is  given  as  a 
normalized  quantity  in  equation  (3-38).   Normalize  (3-39)  to 
obtain 

jx   +  j  tan  6£ 

z  =  _ (3-40) 

1  -  xg  tan  6£  VJ  *V} 
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a) 


b) 


Figure  3.4/1 


End  Coupled  Slot  Line  Resonators 

a)  Physical  Layout  and  Dimensions 

b)  Electric  Field  Component  in  x-Direction 
vs.  z  for  Even  and  Odd  Mode  of  Excitation 
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Since   x„  =  tan(arctan  x„) 


tan(arctan  x  )  +  tan  3£ 
x  =  1  -  tan(arctan  x)-tan  P  (3-41) 


or  with  the  angle-sum  relation  of  the  tangent  function 


x  =  tan(arctan  x  +  ££)  (3-42) 


The  network  description  of  the  circuit  model  in  Figure  3.4/2 
is  accomplished  as  follows. 

Using  the  two  current  loops  in  Figure  3.4/2  the  corre- 
sponding loop  equations  are 


ZI,  +  jwL  L  ±  jcoMI^  =  0 
1    J      c  1       2 


ZI0    +    jwL    In  ±  jooMI,    =    0 
2         J       c    2       J         1 


(3-43) 


or  after  normalizing  with  respect  to  ZQ,  the  characteristic 
impedance  of  the  line 


(jx  +  Jxc)I1  ±  jm   I2 


±jml1  +    (jx+jxc)I2   =    0 


(3-44) 
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Figure  3.4/2   Circuit  Model  for  End  Coupled  Slot 
Resonators 
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Nontrivial  solutions  to  this  set  of  equations  can  be 
obtained  for 


(x+x  J2  -  m2  =  0 
or 

x  =  ±m  -  x  (3-45) 


Define   x   =  -m  -  x  (3-46a) 

x^  =   m  -  x  (3-46b) 

o         c 


Then  the  inductances  of  each  circuit  loop  and  the  mutual 
inductance  of  the  coupling  region  are 


x_  +  x 

x„  = ~ (3-4/a) 

C  2. 


x   -  X 
and     m  =      0  (3-47b) 


Equations  (3-47)  could  be  solved  now  if  x   and  x   are  known. 
These  can  be  computed  with  the  aid  of  equation  (3-42)  using 
the  procedure  outlined  in  Section  3.1  as  follows. 
Let  the  even  mode  of  operation  be  defined  as 
E  (x,z)  =  E  (x,-z)   and  similarly  the  odd  mode  as 
E  (x,z)  =  -E  (x,-z) .   A  plot  of  these  functions  is  shown  in 

X  X 

Figure  3.4/lb.  Given  a  set  of  physical  parameters  and  a 
frequency  of  operation  one  computes  the  resonant  lengths 
L   and  L   for  the  even  and  odd  mode  of  operation  respectively 
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where  the  same  method  for  obtaining  these  applies  as  outlined 
in  Section  3.1  with  equation  (3-2  8)  in  general.   However,  due 
to  the  different  geometry  of  this  problem  the  approximating 
transforms  of  the  electric  field  components  in  equation  (3-2 7b) 
are  now  changed  to 

,~x        77  COS  — y-  T+L 

EKy    =  2  -* 7   •  cos(£  —-—•)  (3-46a) 


(0.  it  cos  — ^  T+L 

<f;  =  j  2  -, 1 .  sin(?  — y^)  (3-46b) 


for  the  even  and  odd  mode  respectively  while  equation  (3-2 7a) 

remains  unchanged.   The  shifting  property  of  the  Fourier 

transform  was  used  as  previously  described  in  detail  in 

Section  2.2. 

Once  the  two  quantities  L   and  L.   are  solved  for  by  a 
a  e      o  J 

numerical  root  seeking  procedure  the  even  and  odd  mode 
reactances  are  then 


x  =  tan(27T  r-r  +  arctan  x  )  (3-47a) 

e  a  s 


xQ  =  tan(2TT  y^  +  arctan  x  )  (3-47b) 


The  final  computation  for  the  reactances  in  the  coupling 
region  with  equations  (3-47)  completes  this  analysis. 
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A  representative  set  of  curves  resulting  from  this 
calculation  is  shown  in  Figure  3.4/3  where  the  following 
dimensions  of  the  shielding  structure  were  assumed: 
a/D  =  5.5,  h1/D  =  h2/D  =  8.0,  e   =  20  and  the  normalized 
frequency  of  operation  is  D/A  =  0.0159,  which  corresponds 
for  the  commercially  available  thickness  of  the  dielectric 
D  =  1/16"  to  a  frequency  of  3  GHz.   As  one  expects,  the 
normalized  reactance  in  the  coupling  region  is  essentially 
that  computed  for  the  end  effect  of  a  single  short  circuited 
slot  line,  while  the  mutual  inductance  describing  the  coupling 
effect  rapidly  decreases  as  the  separation  T  between  the 
ends  of  the  two  lines  increases. 

Although  these  results  are  not  confirmed  yet  by  experiment, 
it  is  believed  that  sufficient  experimental  verification 
exists  in  the  comparisons  of  Section  3.3.   The  analysis  of 
this  section  is  nothing  but  a  simple  extension  of  the  same 
method  and  the  results  should  again  be  accurate. 

3 . 5   SUMMARY 

The  theory  developed  in  Section  2  has  been  extended  for 
the  analysis  of  shielded  coplanar  transmission  lines  and 
resonators.   An  exact  formulation  of  the  problem  via  a  full 
or  hybrid  mode  analysis  is  carried  out  in  the  Fourier 
transform  domain  where  the  transforms  of  dyadic  Green's 
function  components  are  derived.   With  the  application  of 
Galerkin's  method  in  the  transform  domain  and  a  finite 
approximation  of  the  electric  fields  in  the  conductor  plane 
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Figure  3.4/3 


Reactances  x   and  m  for  Two  End  Coupled  Slot 
c 


Lines  vs.  Varying  Separation  with 
a/D  =5.5,  h,/D  =  h2/D  =  8.0,  W/D 

D/X  -  0.0159  and  e   =  20 


0.5, 
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above  the  dielectric  substrate  a  numerical  solution  to  the 
slot  line  resonator  problem  is  obtained.   Although  the 
primary  intention  was  the  analysis  of  a  halfwave  resonator  an 
easy  modification  of  the  theory  unveiled  the  line  parameters 
of  a  shielded  slot-line.   Furthermore,  with  this  data 
available  the  end  effect  of  a  single  short  circuited  slot- 
line  is  theoretically  determined. 

Once  the  Helmholtz's  equations  of  the  two  scalar  potential 
functions  are  solved  for  this  inhomogeneous  boundary  value 
problem  through  a  matching  of  boundary  conditions  in  the 
various  spatial  domains  in  terms  of  the  transforms  of  the 
current  density  or  electric  fields  at  the  conductor  plane 
above  the  dielectric  any  other  configuration  may  be  analyzed. 
The  geometry  of  the  problem  encountered  enters  into  the 
calculation  only  through  the  functional  behavior  of  the 
assumed  currents  or  fields. 

As  an  example  of  this  result  the  coupling  effect  between 
two  end  coupled  slot  resonators  was  analyzed.   Furthermore 
it  is  by  now  clear  that  any  coplanar  transmission  line  within 
the  shielding  structure  may  be  investigated  by  simple 
modifications  of  the  assumed  distributions  similar  to  the 
methods  employed  in  Section  2  for  coupled  slot-lines  and 
coplanar  strip  lines. 

The  theoretical  results  of  this  section  are  new  and 
provide  important  design  information.   An  experimental 
verification  is  always  desirable.   Therefore,  experimental 
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data  for  the  end  effect  of  a  shorted  open  boundary  slot  line 
[Ref.  24]  is  used  in  a  comparison  with  the  theory  developed 
here  for  a  closed  boundary  structure.   If  the  shielding 
walls  are  moved  far  enough  from  the  slot  it  is  intuitively 
clear  that  this  comparison  is  valid.   Excellent  agreement  is 
found  for  various  parameters.   Since  the  resonating  length  of 
the  halfwave  resonator  as  well  as  the  wavelength  ratio  of  the 
transmission  line  are  involved  in  this  computation,  this 
experimental  verification  of  the  theory  is  quite  thorough  . 
Furthermore,  shielded  transmission  line  parameters  tend  to 
approach  the  values  of  the  open  boundary  structure  if  the 
shields  are  further  removed  or  similarly  if  the  frequency  of 
operation  is  increased.   This  also  justifies  the  validity  of 
the  comparison. 

Since  modification  of  this  theory  for  the  analysis  of  the 
end  coupling  effect  between  two  slot  lines  is  of  a  simple 
nature,  further  experimental  verifications  were  not  believed 
to  be  necessary. 

Unfortunately  design  curves  for  any  of  the  configurations 
discussed  in  this  section  cannot  be  given  for  various  reasons, 
Depending  on  a  specific  application  the  dimensions  of  the 
shielding  structure  can  vary  and  hence  will  influence  the 
line  parameters.   Line  parameters  are  not  only  frequency 
dependent  but  depend  on  various  considerations  in  specific 
applications,  for  instance  in  transitions  to  microstrip  lines, 
So  the  parameters  D/X,  W/D,  T/D  and  c   are  determined  by 
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other  constraints.   Hence,  the  computer  program  proves  to 
be  more  versatile  than  design  curves  and  can  be  used  in  any 
specific  application. 
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4 .   SLOT  LINE  DIRECTIONAL  COUPLER  DESIGN 

The  design  of  passive  components  such  as  filters  and 
couplers  in  the  upper  UHF  and  SHF  frequency  bands  is  based 
on  electrical  coupling  between  two  or  more  transmission 
lines.   (Some  examples  are  bandpass  filters,  comb-line  filters, 
interdigital  filters  and  directional  couplers.)   The  following 
analysis  and  experiments  will  be  confined  to  the  study  of 
directional  couplers,  which  are  used  as  attenuators,  power 
splitters  or  hybrid  junctions.   However,  most  commonly 
directional  couplers  are  used  as  a  sampling  device  to  measure 
separately  the  forward  and  backward  waves  on  a  transmission 
line,  where  a  fraction  of  the  energy  traveling  in  either 
direction  is  coupled  out  such  that  the  energy  flow  on  the 
main  line  is  not  disturbed.   Measurements  of  the  reflection 
can  be  used  in  matching  loads  to  a  transmission  line  while 
the  differences  between  forward  and  backward  sampled  power 
indicates  the  net  power  transfer.   The  most  often  used 
structure  for  the  realization  of  a  directional  coupler  is 
a  four-part  network  formed  by  placing  two  transmission  lines 
close  enough  to  each  other  so  that  their  fields  interact 
[Ref .  18] .   Numerous  authors  investigated  these  parallel 
line  couplers  with  strip  or  microstrip  lines.   Although 
only  recently  the  interest  in  the  studies  of  these  shifted 
from  TEM  analyses  to  frequency  dependent  solutions,  one  may 
say  that  the  design  and  performance  is  well  understood  and 
known. 
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Since  the  introduction  of  the  slot  line  as  an  alternative 
transmission  line  in  1969  by  Cohn  [Ref.  3]  there  has  appeared 
only  one  publication  in  the  literature  about  a  slot  line 
coupler.   An  experimental  investigation  by  Mariani  and  Agrios 
[Ref.  23]  points  out  that  co-directional  coupling  can  be 
achieved  with  two  parallel  slot  lines.   This  fact  is  no 
longer  surprising  since  Oliver  in  Reference  2  6  showed  that 
contra-directional  coupling  is  possible  only  for  TEM  lines, 
but  co-directional  coupling  occurs  if  the  phase  constants 
for  the  two  natural  modes  is  different.   That  this  is  the 
case  for  the  slot  line  is  clearly  seen  now  in  Section  2.2. 

Since  this  thesis  provides  for  the  first  time  the 
analytical  data  necessary  to  describe  the  coupling  effect 
between  these  two  lines  the  objective  of  this  section  is  not 
only  a  verification  of  the  thoeretical  results  obtained 
in  Section  2.2  by  experiment  but  is  also  the  first  presenta- 
tion of  a  co-directional  as  well  as  a  contra-directional 
coupler  using  slot  transmission  lines. 

4.1   THEORETICAL  METHODS  FOR  ANALYSIS 

A  schematic  diagram  of  the  directional  coupler  to  be 
analyzed  is  presented  in  Figure  4.1/1.   Two  of  the  terminals 
of  this  four-port  network  are  in  general  connected  to  a  main 
transmission  line  while  the  other  two  provide  the  sampled 
output. 

This  four  port  may  be  described  mathematically  by  either 
the  scattering  matrix  or  the  impedance  matrix  where  the 
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Figure  4.1/1   Schematic  of  Directional  Coupler 
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elements  of  the  matrices  may  be  found  from  basic  trans- 
mission line  theory.   Both  matrix  formulations  will  require 
in  general  that  the  line  parameters  such  as  the  character- 
istic impedances  and  the  phase  velocities  for  the  two  natural 
modes  be  known.   Since  in  the  design  of  directional  couplers 
the  ratios  of  power  at  the  various  ports  are  important  to 
determine  coupler  parameters  such  as  coupling  coefficient  and 
directivity,  a  more  explicit  formulation  in  terms  of  the 
ratios  of  the  voltages  at  the  various  ports  seems  desirable. 
All  four  ports  in  general  will  be  matched  by  the  characteristic 
impedances  of  the  connecting  transmission  lines,  thus  the 
load  impedances  for  the  four-port  network  are  given  by  design 
considerations  and  the  ratios  of  power  can  be  calculated  from 
the  voltage  ratios.   In  the  following  sections  two  mathematical 
methods  are  derived,  an  impedance  matrix  method  and  a  method 
which  explicitly  provides  the  voltage  ratios  at  the  various 
parts  as  well  as  the  input  impedance  at  port  1.   The  latter 
is  used  for  computing  the  theoretical  response  of  the  coupler 
while  the  first  method  has  been  used  as  a  numerical  check  by 
an  alternative  method. 

4.1.1   Impedance  Matrix  Formulation 

According  to  Zysman  and  Johnson  [Ref.  2  7]  the  impedance 
matrix  description  of  the  four-port  network  is 

[V]  =  [Z][I]  (4-1) 


124 


where  [Z]  is  a  4x4  matrix  whose  elements  are  given  as 


Zll  "  Z22  =  Z33  "  Z44 


'*  I  <z0e  cot  6e  +  Z0o  cot  eo} 


Z12  '  Z21  "  Z34  "  Z43 


-j  77  (zn  cot  e  -  z_  cot  e  ) 

J  2    Oe      e    Oo      o 


'13     31     24     42 


•j  77  (Z„   esc  6   -  Z_   esc  6  ) 
J  2    Oe      e    Oo      o 


Z14  "  Z41    Z23    Z32 


j  ■=-  (Z„   esc  6   +  Z_   esc  9  ) 
J    2         Oe      e    Oo      o 


where  6   and  6   are  the  even  and  odd  mode  electrical  lengths 
e      o  3 

of  the  coupling  region. 

Assume  now  that  ports  2,    3  and  4  are  terminated  with 

some  known  impedance  Z   and  a  known  test  voltage  V,  is  supplied 

to  port  1 .   Let  the  input  admittance  at  port  1  be  denoted  by 

Y.  .   Then  a  substitution  of  these  values  in  (4-1)  and  a 
in 

reordering  yields 
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or  in  a  shorthand  notation 

[A] [x]  =  [y]  (4-3) 

The  unknown  vector  [x]  can  now  be  found  by  standard  means 
for  solving  a  set  of  simultaneous  linear  equations.   It  is 
seen  that  the  matrices  are  in  general  complex.   To  avoid 
this  complex  arithmetic  an  isomorphism  between  this  system 
of  four  complex  unknowns  and  a  system  of  eight  real  unknowns 
given  by  Zurmuhl  in  Reference  2  8  was  used.   The  imaginary 
quantities  of  (4-2)  are  then  treated  as  real  numbers. 

However,  this  method  of  describing  the  directional 
coupler  is  not  very  useful  for  design  purposes,  since  not 
enough  explicit  information  about  the  voltage  ratios  or  the 
input  impedance  at  port  1  is  provided  by  analytical  means. 
A  more  suitable  formulation  is  given  in  the  next  section. 
4.1.2   Explicit  Coupler  Equations 

An  explicit  form  of  coupler  equations  was  given  by 
Jones  and  Bolljahn  [Ref .  29] ,  where  the  superposition  prin- 
cipal was  used.   The  input  at  port  1  is  considered  as  the 
superposition  of  the  two  modes  of  excitation,  odd  and  even. 
The  analysis  presented  in  Reference  29  is  restricted  in  that 
it  applies  only  for  equal  phase  velocities  of  the  two  natural 
modes.   By  some  algebraic  manipulation  it  is  possible  to 
derive  a  more  general  form  for  unequal  phase  velocities.   This 
calculation  follows  the  general  approach  outlined  in  [Ref.  29] 
and  details  are  not  presented  here. 
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Assume   that   ports    2    to   4    are   terminated  with 

Z0  =     yZ0e'Z0o  (4-4) 

Then  the  ratio  of  the  three  voltages  at  these  ports  with 
respect  to  port  1  are  given  as 


V~    cos  fl   +  j  —  sin  0      cos  9   +  j  R  sin  6 

_2  =      e    J  R      e   _       o    o    (4  -. 

V1  DE  DO  v    ; 


V3  1  J_ 

V,  DE  DO 


V4 


(4-6) 


vl        --W-.  +        -55-  (4"7) 


where   R  =  \\Zn    /Z 
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DE  =  2cos  9e  +  j (R  +  1/R)  sin  6 


DO  =  2cos  9Q  +  j(R  +  1/R)  sin  0O 


The  input  impedance  at  port  1  is  given  by 


1  -  A  tan  6e  tan  6q  +  j (B  tan  9e  +  C  tan  6  ) 

Zin    Z0  1  -  A  tan  6   tan  6   +  j (C  tan  6   +  B  tan  6  )' 
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where 
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The  necessary  data  for  an  evaluation  of  these  equations  is 
provided  by  the  analysis  in  Section  2.2.   The  investigation 
of  several  configurations  of  directional  couplers  is  now 
possible. 

4.2   THEORETICAL  RESPONSE  OF  SEVERAL  COUPLERS 

Some  basic  design  considerations  and  the  theoretical 

performance  of  several  slot  line  couplers  designed  using 

data  from  Section  2.2  are  discussed  in  this  section.   The 

analysis  of  the  first  coupler  follows  a  similar  investigation 

by  Dalley  [Ref.  30]  for  strip-line  directional  couplers. 

It  is  seen  from  equation  (4-8)  that  a  perfect  impedance 

match  at  port  1  will  be  obtained  if  Z.   =  Z  .   If  any 

in    o 

coupling  between  the  lines  is  assumed  then  Zrt   ^  Z~   and 
^  Oe    Oo 

6^9   and  an  impedance  match  is  obtained  only  for 

tan  6   =  tan  G    or 
e        o 


l   +  n-rr     ,   n  =  1,2  ,  3 ,  .  .  . 
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A  realistic  example  shows,  however,  that  this  design  would 

be  almost  impractical  even  for  n  =  1  since  a  10  dB  coupler 

at  3  GHz  should  require  a  coupling  length  of  approximately 

5"  for  e   =16.   Such  a  coupler  shows  the  interesting 

property  of  power  flowing  mainly  from  port  1  to  port  3  while 

a  tenth  of  the  input  power  is  diverted  to  port  2;  port  4  is 

essentially  isolated  at  the  center  frequency.   This 

theoretical  response  is  shown  in  Figure  4.2/1.   It  can  also 

be  seen  that  this  coupler  is  much  narrower  in  its  bandwidth 

response  than  a  conventional  TEM  line  coupler. 

Another  coupler  design  could  be  based  on  the  fact  that 

the  phase  velocities  of  the  two  natural  modes  do  not  differ 

very  much.   An  examination  of  equations  (4-5)  to  (4-8)  with 

the  assumption  G   +  6_  =  it   for  8  z'  B      unveils  that  one 
c  e    o  e    o 

could  indeed  expect  a  similar  performanc  for  this  coupler 
as  for  a  contra-directional  TEM-line  coupler,  though  not  an 
infinite  directivity. 

A  digital  computer  simulation  was  performed  and  the 
characteristics  of  this  simulated  coupler  are  shown  in 
Figure  4.2/2.   This  coupler  has  a  wider  bandwidth  character- 
istic than  the  previous  one  and  exhibits  directivity  of  better 
than  6  dB.   Due  to  the  loose  coupling  the  assumption 


Zq  =  V/Zq  '?Jn0        ^s  still  valid  and  the  mismatch  at  port  1 
causes  a  VSWR  of  less  than  1.08. 

Since  the  objective  of  this  study  is  an  experimental 
verification  rather  than  an  optimum  coupler  design  no  attempt 
to  obtain  a  perfect  match  was  made. 
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Figure  4.2/1   Theoretical  Response  of  Coupler  with 

6   =  G   +  it,  W/D  =  0.5,  S/D  =  0.3  and  e   =  16 
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Figure    4.2/2      Theoretical    Response   of   Contra-Directional 

Coupler  with   S/D   =0.3,    W/D   =0.5    and   er   =    16 
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A  third  possible  choice  of  design  is  one  where 
0  +  6   =  2tt.   This  coupler  will  perform  as  a  co-directional 
coupler  with  high  directivity  as  can  be  seen  if  the 
following  approximations  are  made.   Assume  loose  coupling, 

then   (Z0e+Z0o»/\/Z0e-Z0o'  «  2   with   Z0e  !  Z0o' 

0  +  0_  =  2 7T   and   0   is  slightly  larger  than  ti  while  0 
e    o  o  e 

is  slightly  less  than  it.   Let   0   -  Tr+a   and   0   =  Tr-a 
3   u  o  e 

where  a  is  the  difference  between  the  electrical  lengths  for 
either  even  or  odd  mode  and  the  electrical  length  of  a 
single  slot  line.   The  physical  length  of  the  coupling  region 
is  still  computed  from   0   +  0   =  2tt. 


Equation  (4-6)  can  then  be  simplified  to 


V3 


Vl 


~  a   and  equation  (4-5)  to 


V2 


Vl 


s  0  .   The  theoretical 


performance  of  this  coupler  is  shown  in  Figure  4.3/2. 

4.3   DESIGN  OF  TWO  COUPLERS 

A  comparison  between  the  previously  derived  theory  of 
a  dispersive  line  coupler  and  its  actua]  performance  was 
made  by  designing  two  slot  line  couplers.   A  general  layout 
for  a  coupler  is  shown  in  Figure  4.3/1.   It  was  decided  to 
use  a  microstrip  to  slot  line  transition  to  couple  energy 
into  and  out  of  the  coupler.   The  theory  of  this  transition 
is  well  understood  because  Knorr  in  Reference  5  derived  a 
circuit  model  for  this  transition  and  developed  a  computer 
program  by  which  an  interative  scheme  yields  the  values  of 
the  various  parameters  required  to  achieve  minimum  VSWR. 
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Figure  4.3/1   Geometry  of  Experimental  Slot  Line  Coupler 
with  Microstrip-Slot  Transitions 
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This  program  was  readily  available  and  used  in  this  study. 
The  theory  is  given  in  detail  in  Reference  5  and  therefore 
not  reiterated  here. 

The  following  design  steps  were  taken  before  the  actual 
physical  construction  could  be  made. 

a)  It  was  decided  to  use  standard  miniature  coaxial 
cable  with  an  impedance  of  50  Ohm  to  connect  the  coupler 
to  the  test  equipment.   This  implied  the  use  of  a  50  Ohm 
microstrip  line,  thus  the  width  of  the  microstrip  was 
determined. 

b)  The  iterative  procedure  described  in  Reference  5  was 
used  to  determine  the  physical  dimensions  of  the  slot  and 
microstrip  line  stubs  and  the  width  of  the  slot  line  at  the 
desired  operating  frequency  of  3  GHz  for  lowest  VSWR. 

c)  With  the  known  dimension  of  the  slot  line  the  digital 
computer  program  described  in  Section  5.2.1  computed  the  line 
characteristics  of  two  coupled  slot  lines  for  varying  spacing 
(varying  coupling)  between  the  two  lines. 

d)  The  data  from  step  c)  were  used  in  another  computer 
program  which  evaluated  the  performance  of  the  slot  line 
coupler  but  did  not  include  the  microstrip  to  slot  line 
transitions.   The  coupling  length  could  be  varied.   Both 
mathematical  formulations  of  Sections  4.1.1  and  4.1.2  were 
used  in  this  design  step. 

e)  After  the  desired  theoretical  performance  was  obtained 
the  lines  were  photo-etched  onto  the  dielectric  substrate. 
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Several  iterations  through  steps  c)  and  d)  are  in  general 
necessary  until  the  requirements  for  the  coupler  are  met. 

In  the  two  experiments  the  coupled  slots  were  etched 
onto  one  side  of  a  1  oz.  copper  surface  on  a  dielectric 
substrate  with  thickness  D  =  0.125"  and  permittivity  e   =16, 
dimensions  4.5"  x  6.5",  while  the  four  microstrip  lines  were 
etched  onto  the  opposite  side,  to  make  all  four  ports 
accessible.   OSM  24  8-4  adapters  performed  the  transition  from 
microstrip  to  the  miniature  coaxial  cable  of  the  test 
equipment.   A  network  analyzer  was  used  for  the  recording  of 
the  coupler  characteristic. 

The  first  coupler  designed  was  the  co-directional  coupler 
shown  in  Figure  4.3/2.   The  discrete  points  in  this  graph 
indicate  the  experimentally  measured  coupler  response.   The 
slot  width  was  0.0625"  and  the  separation  between  the  parallel 
slots  0.135".   It  was  originally  thought  that  the  coupling 
length  could  be  defined  to  be  the  length  between  the  points 
A'  and  A"  shown  in  Figure  4.3/1.   This  length  was  0.375"  and 
was  found  not  to  match  the  experimental  with  the  theoretical 
results.   It  was  learned,  however,  that  good  agreement  could 
be  obtained  if  the  coupling  length  was  taken  between  points 
B1  and  B" ,  in  this  case  0.75"  which  corresponds  to 
8  +  6   -  2tt.   The  waves  on  the  two  transmission  lines  do 
obviously  not  decouple  abruptly  in  this  particular  design 
so  that  the  coupling  region  is  somewhat  extended  to  an 
"effective"  coupling  length.   The  sidearms  to  ports  2  and  3, 
however,  could  not  be  bent  further  or  the  waves  on  these  arms 
would  interfere  with  each  other. 
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Figure  4.3/2   Theoretical  and  Experimental  Response  of 
Co-directional  Coupler  with  S/D  =  1.08, 
W/D  =0.5  and  e   =  16 
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The  second  experiment  was  the  design  of  a  contra- 
directional  10  dB  coupler.   The  slot  width  was  again 
0.0625",  but  the  separation  between  the  slots  was  decreased 
to  0.0375".   The  coupling  length  was  taken  between  points  B' 
and  B"  to  be  0.371".   The  result  of  this  experiment,  shown 
by  the  discrete  points  in  Figure  4.3/3  revealed  some 
disagreement  with  the  theoretical  performance  in  Figure 
4.2/2.   Although  the  experimental  curves  include  the  micro- 
strip  to  slot  line  transition,  these  transitions  had  a 
return  loss  of  better  than  10  dB  over  most  of  the  frequency 
band  and  could  not  have  caused  this  disagreement.   Again  it 
was  found  that  the  critical  design  parameter  was  the  coupling 
length.   Due  to  the  closer  coupling  in  this  second  coupler 
design,  the  coupling  length  had  to  be  extended  by  2  0%  beyond 
the  length  B1  -  B"  to  match  experiment  and  theory.   Figure 
4.3/3  shows  the  comparison  between  the  experimental  and  the 
adjusted  theoretical  response,  where  now  better  agreement  is 
achieved. 

Although  initially  intended,  design  curves  for  slot 
line  couplers  cannot  be  presented  since  the  slot  line, 
single  or  coupled  configuration,  is  a  dispersive  transmission 
line  with  frequency  dependent  characteristics.   Thus  any  set 
of  design  parameters  would  be  valid  for  a  narrow  frequency 
band  only. 
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Figure  4.3/3   Theoretical  and  Experimental  Response 

of  a  Coupler  with  S/D  =0.3,  W/D  =0.5 

and  c   =16 
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4 . 4   SUMMARY 

An  analytical  evaluation  of  slot  line  directional 
couplers  has  been  presented.   Due  to  the  different  phase 
velocities  of  the  even  and  odd  modes,  co-  or  contra- 
directional  coupling  can  be  achieved.   Theoretical  results 
were  compared  with  experimental  data  by  the  design  of  two 
couplers.   It  was  found  that  the  coupling  length  is  a  critical 
parameter  in  this  design.   Coupled  lines  were  assumed  to 
extend  infinitely  in  the  computation  of  the  line  parameters 
in  Section  2.2.   However,  in  a  coupler,  abrupt  decoupling 
of  the  two  lines  does  not  occur  due  to  surface  currents  on 
the  metal  and  this  tends  to  increase  the  effective  coupling 
length.   This  result  is  also  new  and  will  be  important  in 
future  designs. 

With  an  adjustment  of  the  coupling  length  good  agreement 
between  the  theoretical  and  experimental  results  further 
verifies  the  spectral  domain  technique  used  in  Section  2  in 
general  and  the  results  of  Section  2.2  in  particular. 
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5 .   COMPUTER  PROGRAMMING 

In  Sections  2  and  3  the  theoretical  background  of  the 
analyses  were  provided  together  with  some  examples  of  the 
results  obtained.   By  no  means  do  these  examples  furnish 
enough  information  for  any  design  purposes.   Since  this 
study  presents  the  first  extensive  investigation  of  coplanar 
transmission  lines,  data  aiding  in  the  design  are  needed. 
Extensive  design  curves  would  not  only  be  beyond  the  scope 
of  this  thesis  but  are  even  sometimes  impossible  to  provide 
as  stated  before.   For  instance,  for  parallel  and  end  coupled 
slot  lines,  not  all  applications  can  be  foreseen  at  this 
point.   Hence,  a  great  effort  was  devoted  to  the  development 
of  two  versatile  computer  programs  which  offer  a  variety  of 
choices  for  the  analysis  of  open  and  shielded  boundary  co- 
planar  transmission  lines.   It  is  believed  that  these  will 
be  of  more  value  to  the  microwave  industry  than  any  design 
curves  could  be.   Before  the  two  programs  are  described  some 
general  considerations  which  led  to  their  development  are 
necessary  for  the  understanding  of  the  numerical  methods 
used. 

5.1   REMARKS  ON  NUMERICAL  ANALYSIS 

The  computation  of  the  characteristic  impedance  and  the 
perturbational  analysis  in  Section  2  is  based  upon  the  solution 
to  the  dispersion  characteristic  problem  for  the  transmission 
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line  under  consideration.   In  other  words,  the  wave  propaga- 
tion constant  or  wavelength  ratio  X'/A  must  be  known  before 
any  other  investigations  can  be  started  since  only  in  this 
case  are  the  scalar  potential  functions  in  the  transform 
domain  known.   These  are  necessary  for  the  evaluation  of 
the  power  distribution  and  for  the  perturbational  analysis 
for  a  ferrite  substrate.   Hence,  in  the  analysis  of  slot  lines 
the  starting  point  is  a  solution  for  the  propagation  constant 
3  to  equation  (2-19) ,  while  for  coplanar  strip  lines  equation 
(2-32)  applies. 

In  general  a  8  is  sought  such  that 

2 

/  G(a,8)  |T(a) |   da  =  0  (5-1) 

—  oo 

where  G(a,8)  is  a  component  of  the  transform  of  the  dyadic 
Green's  function  N..  (a,B)  or  M.(a,8)  with  equations  (2-19) 
and  (2-32)  respectively  and  T(a)  is  the  Fourier  transform 
of  the  electric  field  or  current  density  assumed  in  the 
plane  y  =  D  for  slot  or  coplanar  strip  line  respectively. 
A  first  order  approximation  is  assumed  here;  if  a  second  or 
higher  order  approximation  is  desired  the  zero  of  the  deter- 
minant of  the  appropriate  coefficient  matrix  must  be  found. 
It  is  obvious  that  an  analytic  solution  in  closed  form 
for  6  cannot  be  obtained  due  to  the  algebraic  complexity 
leading  to  the  quantity  G(a,3)  as  shown  in  Appendix  A.   Thus 
a  numerical  solution  is  mandatory  and  must  be  obtained  by 
programming  a  digital  computer. 
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Several  observations  must  be  made  to  adapt  the  theoret- 
ical formulations  in  Section  2  to  this  numerical  solution 
via  a  digital  computer. 

It  is  shown  in  Appendix  A  that  the  transforms  of  the 
dyadic  Green's  function  components  M,  -  M.  or  N,  -  N^  are 
given  in  terms  of  coefficients  F,  through  L.  which  are 
obtained  by  definitions  during  the  application  of  the 
boundary  matching  technique.   These  coefficients  are  either 
real  or  imaginary  so  that  in  general  a  complex  computation 
on  the  computer  seems  to  be  necessary.   Furthermore,  it  is 
shown  in  Section  2.1.1  that  for  all  values  of  the  transform 


V"?  9^ 

k2  -  3    the  nature  of  the  solutions  in 

equations  (2-12)  and  (2-13)  changes  from  hyperbolic  to 

trigonometric  functions.   For  example,  sinh  y^D      changes 

to   j  sin  y-D  .   Since  these  functions  enter  into  the 

computation  of  the  quantities  F-.  to  N.  it  is  found  that  these 

quantities  are  not  only  pure  real  or  pure  imaginary  but 

change  from  real  to  imaginary  or  vice  versa  for  the  two 

2  2 

cases   y9   >  0   and   y2   <  0  . 

This  fact  furthermore  seems  to  suggest  a  complex  analysis 

on  the  computer.   However,  on  the  other  hand,  the  disadvantages 

of  a  complex  arithmetic  with  digital  computation  are  well 

known.   Not  only  is  the  required  amount  of  memory  storage 

twice  as  large  as  for  a  real  computation;  a  fundamental 

setback  is  that  the  computation  times  for  instructions  like 

multiplications,  divisions,  and  evaluations  of  trigonometric 

and  hyperbolic  functions  with  complex  arguments  are  on  the 
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average  four  to  five  times  longer  than  in  real  arithmetic. 
An  avoidance  of  digital  complex  computation  is  therefore 
highly  desirable  and  can  be  obtained  if  some  algebraic 
manipulations  with  the  quantities  F,  to  N.  of  Appendix  A  are 
performed  prior  to  the  computer  programming. 

Before  the  details  of  this  manipulation  are  shown  a 
remark  with  regard  to  a  scaling  technique  implied  in  the 
following  derivation  must  be  made.   Dielectric  substrates 
are  commercially  available  for  various  thicknesses  D  so  that 
in  general  the  line  parameters  should  be  computed  for  each 
thickness  and  varying  frequency.   A  more  generalized  method 
of  analysis  is  simply  to  scale  the  various  geometric 
quantities  involved  with  respect  to  this  thickness  to  obtain 
more  versatile  parameters  for  a  varying  normalized  frequency 
D/X.   The  resulting  design  curves  are  then  applicable  for  any 
geometry  independent  of  the  specific  dielectric  thickness. 
A  manipulation  of  equations  (A-13)  to  (A-24)  including  this 
scaling  technique  is  given  in  Appendix  G,  where  it  is  shown 
that  the  avoidance  of  complex  arithmetic  on  the  computer  and 
thus  long  computation  times  can  be  traded  against  a  somewhat 
lengthy  but  straightforward  algebraic  computation  by  hand 
before  the  programming  is  attempted. 

With  this  formulation  of  the  dyadic  Green's  function  the 
further  steps  to  obtain  the  wavelength  ratio  X'/X  are  as 
follows.   For  a  first  order  approximation  equation  (5-1)  has 
to  be  solved  where  now  due  to  the  scaling  technique  the 
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variable  of  integration  is  aD.   The  transform  of  the 

assumed  electric  field  (or  current  density)  for  the  geometry 

of  transmission  line  to  be  analyzed  is  given  in  Section  2. 

Since  in  (5-1)  the  absolute  value  squared  of  these  functions 

is  to  be  taken,  this  part  of  the  integrand  will  be  in 

general  a  real  and  even  function  of  aD.   Preliminary  numerical 

investigations  showed  that  the  components  M,  ,  M« ,  N,  and  N. 

are  also  even  functions  of  aD  so  that  the  whole  integrand 

in  (5-1)  is  even,  a  fact  which  allows  further  savings  in 

computation  time  since  a  numerical  integration  over  a  half 

interval  between  the  integration  limits  is  sufficient.   Two 

representative  examples  of  the  real  parts  of  the  components 

M4  and  N,  and  their  dependence  on  aD  are  displayed  in 

Figures  5.1/1  and  5.1/2  respectively.   Although  these 

components  according  to  the  derivation  in  Appendix  G  are 

imaginary,  this  fact  is  immaterial  if  a  first  order 

approximation  is  used. 

For  a  second  order  approximation,  however,  the  following 

observation  must  be  made  as  explained  here  for  instance  with 

reference  to  equation  (2-19) .   According  to  equations  (2-21) 

or  (2-23)  E  (a)  is  an  even  and  real  function  of  a  while 

E  (a)  is  odd  and  imaginary.   N,  and  N„  are  found  to  be  even 
z  J.        4 

functions  while  N~  and  N~  are  odd  functions;  all  are 
imaginary.   Equation  (2-19)  can  be  rewritten  in  a  simplified 
form  after  the  integration  with  respect  to  a  is  performed  as 
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Figure  5.1/1   M^(aD,X'/A)  vs.  aD  with   er  =  12  and  D/X  =0.03 
for  Five  Arbitrary  Values  X ' /A 
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Figure  5.1/2   N£(aD,A!/\)  vs.  aD  with  er  =  16  and  D/X 
For  Five  Arbitrary  Values  \* /\ 


=  0.03 
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janu  -  bn12  =  0 

(5-2) 
jan21  -  bn22  =  0 

where  the  r\..    are  real  quantities.   A  solution  for  the  two 
coefficients  a  and  b  of  the  basis  functions  (other  than  the 
trivial  solution)  does  exist  if  rinHnn  =  ^i?*1"1?!   ^n  wn;*-cn 


case   ja  =  b 


n12 
nll 


If  b  is  assumed  to  be  real  this  clearly  indicates  that 
a  must  be  imaginary  or  vice  versa,  in  other  words,  the  electric 
field  components  in  the  x  and  z  directions  are  9  0°  out  of 
phase.   Although  this  result  is  automatically  implied  through 
the  notation  used  in  the  dispersion  characteristic  calculations 
it  should  receive  special  consideration  in  the  power 
calculation.   Since  only  a  first  order  approximation  for 
this  computation  is  used  in  the  present  work,  this  phase 
shift  was  of  no  concern. 

A  desired  solution  for  equation  (5-1)  can  now  be  found 
by  assigning  an  arbitrary  numerical  value  for  the  wavelength 
ratio  such  that   1  <  A/A1  <   e    and  then  iteratively 
changing  this  until  the  integration  yields  zero  or  more 
precisely  a  small  enough  value  so  that  a  prescribed  accuracy 
for  the  wavelength  ratio  is  realized. 

Several  representative  graphs  of  integrands  for  equation 
(5-1)  are  shown  in  Figures  5.1/3-6  where  the  coefficients 
displayed  in  Figures  5.1/1-2  are  used.   Figure  5.1/3  shows 
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Figure  5.1/3   N" (aD, X ' /X) | E    (aD)  |   vs.  aD  with   W/D  -  0.5, 
D/X  =  0.0  3  and  er  =  16 
several  values  of  X'/X 


D/X  =  0.03  and  e   =  16  (single  slot)  for 
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Figure  5.1/4   NJ  (aD,  X  '  /X)  |  E^  (o:D)  |  2  vs.  aD  with  W/D  =  0.5, 
D/X  =  0.03,  S/D  =1.0  and  er  =  16  (coupled 
slots,  even  mode)  for  several  values  of  X'/X 
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Figure  5.1/5   Ni'  (aD,  X1  /A)  |  E  (aD)  |   vs.  aD  with  W/D  =  0.5, 
D/X  =  0.03,  S/D  =1.0   and  er  =  16  (coupled 
slots,  odd  mode)  for  several  values  of  X'/X 
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Figure  5.1/6   MJ (aD, X'/*) | J- (aD) |   vs.  aD  with 


W/D  =  S/D  =1.0,  D/X  =0.03  and  e   =  12 
(coplanar  strip  line)  for  several  values 
of  A'/A 
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the  integrand  for  a  single  slot  line,  5.1/4-5  two  parallel 
coupled  slot  lines  in  the  even  and  odd  mode  respectively, 
and  5.1/6  shows  the  integrand  for  a  coplanar  strip  configura- 
tion.  In  all  cases  the  zero  order  Bessel  function  was  used 
as  the  basis  for  the  transform  of  the  assumed  electric 
field  (current  density)  distribution. 

The  next  task  is  the  preparation  of  the  appropriate 
equations  for  the  numerical  evaluation  of  the  time  average 
power  flow.   Here  again  complex  arithmetic  expressions  can 
be  avoided  if  manipulations  similar  to  those  in  Appendix  G 
are  performed.   General  expressions  for  the  time  average  power 
flow  are  given  in  Appendix  B.   Appendix  H  shows,  how  a  real 
arithmetic  evaluation  of  these  expressions  can  be  accomplished 
for  the  open  boundary  slot  and  coplanar  strip  lines. 

The  normalization  of  equations  with  respect  to  the 
thickness  D  of  the  dielectric  (ferrite)  substrate  and  the 
technique  used  to  avoid  complex  arithmetic  are  given  by  now 
with  sufficient  explanation  and  detail  so  that  no  further 
expansion  seems  necessary.   The  perturbational  expression 
for  the  ferrite  analysis  in  Appendix  D,  for  instance,  need 
only  be  normalized.   Appropriate  remarks  have  already  been 
made  with  regard  to  real  and  imaginary  quantities  involved. 
Hence  equation  (D-2)  will  require  a  real  arithmetic  only. 

Furthermore,  although  so  far  details  are  given  for  open 
boundary  lines  only,  the  manipulations  described  obviously 
apply  for  the  shielded  coplanar  lines  as  well.   The  details 
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of  the  preparation  of  these  equations  for  programming  are 
given  in  Appendix  J.   It  is  observed  that  for  the  dispersion 
characteristic  computations  the  infinite  integration  in 
equation  (5-1)  is  replaced  by  an  infinite  series  due  to  the 
discrete  nature  of  the  transform  variable  a  .   This  remark 
holds  also  for  the  computation  of  time  average  power  flow. 

Preliminary  numerical  investigations  of  the  coefficients 
of  these  series  indicated  an  even  distribution  with  respect 
to  the  variable  a   so  that  computation  time  can  be  saved  by 
summing  over  a  half  interval  only. 

This  result  even  holds  for  computations  concerning  the 
length  of  a  halfwave  resonating  slot,  where  equation  (3-2  8) 
given  in  a  simplified  form  as 


CO         00 

X    /  M1(an,^)  |T(an,U  |   d£  =  0  (5-3) 

n— — °°  — oo 


has  to  be  evaluated.   Instead/  due  to  the  even  behavior  of 
the  integrand  with  respect  to  E,  and  a   the  following 


computation  suffices: 


/   M,  (0,fj  |T(0,U  1 2d£  +  2   E    /   M,  (a.O  |T(a,0  |  2d£  =  0 
0    X  n-1   0    X 

(5-4) 


After  this  preparation  of  the  equations  for  computer 
programming  the  only  remaining  task  is  to  investigate  the 
numerical  integration  with  regard  to  its  limits  and  the 
procedure  to  be  used.   One  of  the  limits  of  the  integrals 
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discussed  so  far  extends  to  infinity  so  that  due  to  the 
finite  presentation  of  numbers  on  a  digital  computer  a  finite 
approximation  has  to  be  used  instead.   From  the  representative 
examples  shown  in  Figures  5.1/1-6  it  is  seen  that  the  inte- 
grands will  decay  very  rapidly  in  general.   Although  no 
analytic  method  can  be  applied  to  estimate  the  order  of 
magnitude  of  a  finite  and  adequate  upper  limit,  a  trial  and 
error  procedure  by  numerical  means  showed  that  for  the  open 
boundary  lines  a  limit  of  ctD  =  55  and  for  the  shielded  lines 
£D  =  9  is  sufficiently  large  to  be  considered  infinite  for 

all  practical  purposes  since  the  integrands  are  at  least  down 

-5 
by  a  factor  of  10    relative  to  their  largest  value  close  to 

the  origin.   Several  different  investigations  showed  that 

these  limits  guarantee  at  least  a  three  digit  accuracy  for 

any  practical  set  of  line  parameters. 

After  establishing  appropriate  limits  of  integration 

the  remaining  problem  of  which  numerical  integration  scheme 

to  use  must  be  considered.   Many  schemes  are  supplied  through 

the  standard  literature  on  numerical  analysis.   In  this  study 

only  three  methods  were  investigated  seriously.   Two 

integration  subroutines  are  supplied  as  part  of  ' IMSL * 

Library  [Ref.  31]  with  names  'DRMBIU'  and  'DCADRE'  where 

both  use  a  repeated  interval-halving  technique,  also  referred 

to  as  a  Romberg  integration.   'DCADRE'  is  an  improved  version 

of  'DRMBIU1  since  it  provides  for  an  adaptive  method  and 

thus  may  save  computation  time.   This  routine  is  quite  new 

and  became  available  only  very  recently.   It  is  used  for 

comparative  purposes  in  this  discussion. 
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An  alternate  integration  scheme  is  the  Gaussian- 
quadrature  method.   According  to  Reference  32  the  Gauss- 
Legendre  formulas  are  significantly  more  accurate  than  the 
equal-interval  formulas  such  as  the  Romberg-method  for  a 
comparable  number  of  sub-intervals  between  the  integration 
limits  with  the  condition  that  the  magnitudes  of  higher  order 
derivatives  of  the  integrand  decrease  or  at  least  do  not 
increase.   In  other  words,  an  equal-interval  formula  will 
in  general  need  more  divisions  into  sub-intervals  to  acquire 
the  same  accuracy  which  results  in  a  greater  computation 
time.   Since  the  computer  programs  presented  here  are  intended 
to  be  used  very  often  for  design  purposes,  the  computation 
time  is  an  important  factor.   Considerable  effort  was  devoted 
therefore  to  finding  the  optimum  numerical  integration  scheme 
where  "optimum"  means  the  shortest  possible  computation  time 
for  a  prescribed  accuracy  of  the  result. 

Due  to  the  algebraic  complexity  of  the  integrands  the 
behavior  of  its  higher  order  derivatives  cannot  be  studied 
on  an  analytical  basis.   It  is  observed,  however,  that  these 
integrands  will  partly  consist  of  the  transforms  of  the 
dyadic  Green's  function  components.   The  examples  shown  in 
Figures  5.1/1  and  5.1/2  indicate  a  smooth  and  continuous 
behavior  with  at  most  one  zero  for  aD  in  the  interval  [0,  ) . 
The  other  part  of  the  integrand  will  consist  of  the  transform 
of  the  assumed  field  (current  density)  distribution.   If, 
for  instance,  coupled  slots  in  the  odd  mode  are  considered, 
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this  function  will  be  in  the  form  of 

.  ,S/D  +  W/D,   ^   D  /W/D  •  otD,     .        _ ,  . 
sin(- 2 — )  aD)  *  Bo* — 2 ^  *       y  =   ^  wnere  y   1S 

the  largest  integer  less  than  or  equal  to  x.   Then  the  first 

55 
term  above  has  E  (-^—  (S/D  +  W/D)  )  zeros  for  aD  in  the  interval 

55 
(0,55].   The  second  term  has  approximately  E  (-p—  W/D)  zeros 

since  for  large  arguments  s,  BQ (x)  ~  cos(x).   These  zeros 

are  not  necessarily  coincident  with  any  of  the  first  ones 

so  that  the  total  number  of  zeros  of  the  integrand  is  the 

sum  of  both. 

For  resonant  slot  length  computations  that  part  of  the 
integrand  involving  the  transform  variable  £D  is  proportional 
to   cos(- — *    )  (see  equation  (3-27b)),  a  function  which 
will  have   E (  -—4 -)  zeros  for  £D  in  the  interval  (0,9]. 

Since  in  general  the  integrands  will  be  basically  in  the 
form  of  a  product  of  a  trigonometric  and  a  zero  order  Bessel 
function  the  following  test  was  performed  to  find  an  efficient 
numerical  integration  scheme  for  the  specific  needs  in  this 
analysis. 

A  test  function  to  be  integrated  in  the  interval  [0,30] 
was  chosen  to  be   f(x)  =  sin (5x) 'BQ  (2x)   which  has  68  zeros 
within  this  interval.   If  a  three  digit  accuracy  is  desired 
a  numerical  Romberg  integration  with  'DRMBIU'  needs 
approximately  6  seconds  computation  time  while  the  cautious 
adaptive  Romberg  method  in  'DCADRE'  needs  only  4  seconds. 
On  the  other  hand,  a  Gaussian-quadrature  80  points  integration 
requires  only  0.5  seconds  to  guarantee  the  same  result. 
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The  advantage  of  this  integration  scheme  becomes  obvious 
then.   However,  in  contrast  to  the  Romberg  integrations, 
where  the  user  can  simply  specify  the  accuracy  of  the  result, 
this  is  not  possible  with  a  Gaussian-quadrature  routine. 
Here  the  user  must  have  advance  knowledge  of  the  general 
behavior  of  the  integrand  to  select  the  number  of  quadrature 
points.   Still,  no  estimate  of  the  integration  error  can  be 
obtained. 

But  from  the  numerical  example  above  it  was  conjectured, 
that  in  the  present  applications  a  quadrature  integration  of 
sufficient  accuracy  can  be  achieved,  if  the  number  of  quadra- 
ture points  is  at  least  twice  the  number  of  zeros  of  the 
integrand  within  the  limits  of  integration. 

This  number  depends  on  the  parameters  such  as  slot  width, 
separation  between  slots  or  strips  and  for  the  resonator 
analysis  on  the  frequency  of  operation  and  will  usually  vary 
over  a  wide  range.   A  basic  feature  in  the  first  program 
'COPLAN'  is  therefore  a  short  computation  which  estimates 
the  number  of  zeros  of  the  integrand  according  to  which  the 
selection  of  an  appropriate  quadrature  routine  is  made.   It 
was  found  by  several  test  runs  and  comparisons  that  a  sub- 
stantial amount  of  computation  time  could  be  saved  without 
any  sacrifice  in  accuracy. 

5.2   COMPUTER  PROGRAM  'COPLAN'  FOR  OPEN  BOUNDARY  COPLANAR 
TRANSMISSION  LINES 

A  numerical  analysis  of  open  boundary  coplanar  trans- 
mission lines  is  provided  by  the  digital  computer  program 
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'COPLAN'  and  its  accompanying  subroutines  'COMBO'  and  'GAUSQ', 
all  shown  in  Appendix  M,  written  in  the  FORTRAN  IV  programming 
language.   The  Gaussian-quadrature  integration  routines  are 
not  shown. 

Although  during  the  development  of  the  theory  in  Section 
2  of  this  thesis  separate  programs  were  written  for  different 
line  configurations,  this  program  'COPLAN'  is  a  combination 
of  all  of  them  and  offers  a  high  degree  in  versatility.   Dis- 
persion characteristic  and  characteristic  impedance  may  be 
computed  for  single  or  coupled  slot  lines,  for  coplanar 
waveguides  or  for  coplanar  strip  lines .   Except  for  the 
latter  the  perturbational  analysis  for  phase  shift  and  attenu- 
ation may.be  applied  as  well  for  ferrite  substrates. 

The  necessary  input  data  and  some  further  information 
about  the  use  of  this  program  are  described  in  the  preceeding 
commentary  section  of  the  main  program.   Input  and  output 
formats  are  chosen  such  that  design  curves  may  be  easily 
generated  if  so  desired. 

After  the  input  data  is  read  in  and  printed  in  an  echo 
check,  preliminary  computations  and  initializations  of  the 
various  computation  loops  are  performed.   The  first  computa- 
tion inside  the  inner  loop  is  a  check  whether  the  operating 
frequency  is  above  the  cut-off-frequency  of  the  lowest  order 
surface  wave  mode;  if  so  the  job  is  terminated.   The  next 
step  is  to  obtain  an  estimate  of  the  number  of  zeros  to  be 
encountered  in  the  interval  of  integration,  according  to 
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which  a  numerical  value  is  assigned  to  an  integer  variable 
to  be  used  in  the  selection  of  the  appropriate  Gaussian- 
quadrature  integration  routine.   This  selection  and  the 
actual  integration  are  performed  in  the  subfunction  'GAUSQ1 
whose  use  is  self-evident. 

The  program  proceeds  with  a  zero  finding  procedure  to 
find  the  wavelength  ratio  A '/A,  according  to  equations  (2-19) 
and  (2-31) ,  where  1  <  A/A'  <    e   .   Integration  limits  are 
computed  to  determine  for  which  values  of  ccD  the  trigono- 
metric and  hyperbolic  functions  in  the  solutions  to  Helm- 
holtz's  equation  apply  (Section  2.1.1).   For  hyperbolic 
functions  (k  ~D  <  | aD |  <  55)  the  appropriate  numerical  inte- 
gration  is  performed  via  the  subroutine  'GAUSQ'  while  for 
0  <_  |  aD  |  <  k  ?D  it  was  found  that  a  24  point  Gaussian-quadra- 
ture  routine  is  accurate  enough  for  all  practical  parameters. 
The  subroutine  COMBO  provides  the  appropriate  integrads  via 
several  identifiers  or  flags. 

After  the  integration  is  carried  out  the  result  is  used 
to  select  a  new  value  for  A' /A  through  a  bisectional  zero 
finding  to  improve  the  result  and  the  procedure  starts  over 
again.   The  computation  is  terminated  when  the  absolute  value 
of  the  difference  between  the  new  and  the  old  value  for 

-A 

A/A'  is  less  than  5*10    and  thus  a  three  digit  accuracy 
after  the  decimal  point  is  guaranteed. 

Preliminary  computations  follow  for  the  determination 
of  the  characteristic  impedance.  Depending  on  the  trans- 
mission line  configuration,  mode  of  operation  and  the  spatial 
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region  the  appropriate  integrands  are  provided  again  through 
the  subroutine  ' COMBO '  in  accordance  to  the  equations  of 
Section  2.   The  numerical  integration  of  the  power  density 
in  the  spectral  domain  follows  the  same  considerations  given 
above  for  the  dispersion  calculations.   The  proper  definition 
of  the  characteristic  impedance  for  the  transmission  line 
to  be  analyzed  is  then  used  to  evaluate  this  line  parameter. 
A  three  digit  accuracy  is  guaranteed  in  this  computation. 

The  two  line  parameters  for  a  dielectric  substrate  are 
thus  determined  and  are  printed  into  a  tabulated  output.   If 
a  perturbational  analysis  for  a  ferrite  substrate  is  desired, 
however,  the  computer  program  proceeds  with  preliminary  com- 
putations to  evaluate  the  external  susceptibility  tensor 
components .   These  are  then  provided  to  the  subroutine 
'COMBO'  which  computes  the  appropriate  integrands  for  the 
perturbational  analysis  of  the  attenuation  and  the  phase 
shift  for  waves  travelling  in  either  direction.   After  the 
numerical  integration  is  performed  some  further  manipulations 
are  necessary  to  allow  the  attenuation  to  be  printed  in  the 
unit  [dB/cm]  and  the  phase  shift  as  a  ratio  of  the  ferrite 
phase  constant  over  the  dielectric  phase  constant.   These 
results  are  printed  together  with  the  results  for  a  dielectric 
substrate.   The  procedure  is  then  reiterated  for  the  next 
set  of  parameters. 

The  subroutine  'COMBO'  provides  the  integrands  for  the 
dispersion  characteristic  and  power  density  computations  as 


160 


well  as  for  the  perturbation  analysis.   The  entry  variable 
of  this  routine  is  the  normalized  transform  variable  ocD 
while  other  coefficients  are  exchanged  in  common  storage 
areas  (COMMON  statement)  with  the  main  program.   Depending 
on  the  transmission  line  configuration  the  appropriate  trans- 
form of  the  assumed  electric  field  (current  density)  is 
computed  by  the  use  of  a  Bessel-f unction  subroutine.   The 
further  evaluations  of  the  respective  integrands  follows  the 
general  theoretical  outline  of  Section  2.   Variable  names 
are  chosen  such  that  these  equations  can  be  easily  compared 
with  the  derivations  of  equations  in  the  appropriate  appen- 
dices so  that  no  further  explanations  are  necessary. 

Representative  computer  outputs  are  given  in  Tables  I 
through  VII.   Waves  travelling  on  a  ferrite  substrate  shew 
a  nonreciprocal  behavior.   It  is  therefore  pertinent  to  cite 
the  orientation  of  the  coordinate  system  and  the  biasing 
field.   Following  the  convention  of  Section  2.4  the  direction 
of  the  coordinates  are  given  in  Figure  1/2  and  the  biasing 
magnetic  field  is  in  the  positive  x-direction.   A  forward 
wave  then  travels  into  the  positive  z  direction. 

Finally  some  statistics  about  this  program  seem  to  be 
apt.   The  total  program  including  all  subroutines  requires 
182  K  bytes  storage  space  and  78  K  bytes  for  the  execution 
step.   To  obtain  a  fast  execution  time  the  program  was  com- 
piled with  the  FORTRAN  H  compiler  which  applies  an  optimiza- 
tion code.   The  computation  time  for  execution  vary  greatly 
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OISP.  AND  CHARACT. 

MPEDANCE  OF  SINGLE  SLOT  LINE 

W/D    1  EPSILON  |    D/L       LP/L 

ZO 

0.100 

11.0 

0.020 

0.4572       57.6 

0.  100 

11.0 

0.030 

0.4446 

59.9 

0.  100 

11. 0 

0.040 

0.4347 

ol.3 

0.100 

11.0 

0.050 

0.4264 

61.9 

0.  100 

11.0 

0.060 

0.4191 

61.  7 

2.000 

11.0 

0.020 

0.5599 

153.6 

2.000 

11.0 

0.030 

0.52  44 

173.6 

;.ooo 

11.0 

0.040 

0.4959 

187.6 

2.000 

11.0 

0.0  50 

0.4718 

195.7 

2.000 

11.0 

0.060 

0.4512 

197.9 

TABLE  I.   Single  Slot  Line  on  a  Dielectric  Substrate 


COUPLED  SLOTS  ON  D 

ELECTRIC  SUBSTRATE 

EVEN 

MODE 

ODD 

MODE 

W/D    |    S/D    !  EPSILON  |    D/L 

LP/L    1    ZO        LP/L        ZO 

0.400 

0.300 

16.0 

0.010 

0.4754 

44  .2 

0.3496 

111.2 

0.400 

0.300 

16.0 

0.030 

0.4065 

55.6 

0.3458 

111.7 

0.500 

0.300 

16.0 

0.010 

0.4861 

47.3 

0.3519 

118.  1 

0.500 

0.300 

16.0 

0.030 

0.4125 

50.6 

0.3  4  72 

118.9 

0.400 

0.600 

16.0 

0.010 

0.h844 

47.5 

0,3537 

93.3 

0.400 

0.600 

16.0 

0.030 

0.4113 

61  .0 

0. 3484 

9'..  7 

0.500 

0.600 

16.0 

0.010 

0.4951 

50.6 

0.3565 

100.5 

0.500 

0.600 

16.0 

0.030 

0.4172 

56.1 

0.3501 

101.9 

TABLE  II.   Coupled  Slot  Lines  on  a  Dielectric  Substrate 
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DI 


SPFRSION    AND    CHARACT . IMPEDANCE    OF    CD-PLANAR    WAVEG 

+ + + + +• 

/D  |  S/D  |     EPSILON     I  C/L  I        L  P/ L  | 


JIDE        | 

b8.5 
7  1.1 
72  .7 
72.8 
61.2 
64.1 
66.1 
66.4 


_W 

1. 
1. 
1. 
1. 
1. 
1. 
1. 
1. 


500 

0.500 

20.0 

500 

0.500 

20.0 

500 

0.500 

20.0 

500 

0.500 

20.0 

500 

i.000 

20.0 

500 

1.000 

20.0 

500 

1.000 

20.0 

500 

1.000 

20.0 

0.C10 

0.3497 

0.020 

0.S412 

0.030 

0.3313 

0.040 

0.3208 

0.010 

0.3599 

0.020 

0.3495 

0.030 

0.3  3  76 

0.040 

0.3252 

TABLE  Hi.   Coplanar  Waveguide  on  a  Dielectric  Substrate 


1   DISPERSION  AND  CHARACT . IMPEDANCE  OF  CO-PLANAR  STRIPS 

i    W/D    1     S/D    1  EPSILON  I    D/L    1   LP/L 

10 

1.500 

0.200 

20.0 

0.010 

0.3582 

58.4 

1.500 

0.200 

20.0 

0.030 

0.33^2 

49.0 

2.000 

0.200 

20.0 

0.010 

0.3794 

59.8 

2.000 

0.2C0 

20.0 

0.030 

0.3455 

47.1 

1.500 

C.400 

20.0 

0.010 

0.3626 

65  .4 

1.500 

0.400 

20.0 

0.030 

0.3359 

54.2 

2.000 

0.400 

20.0 

0.010 

0.3835 

65.8 

2.000 

0.400 

20.0 

0.030 

0.3479 

51.3 

TABLE  IV.   Coplanar  Strip  Line  on  a  Dielectric  Substrate 
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SINGLE     SLOT    LINE    ON    FERRITE     SUBSTRATE 

THICKNESS    OF    DIELECTRIC    =    O.IOO    CM         W/D    =    0.500  S/0    =0.0 

PERMITTIVITY    =    20.0 
SATURATION    MAGNETISATION    =    2300.00    GAUSS 

MAGNETIC    BIASING    FIELD    =       683.60    OERSTED 
LINEWIDTH    =       100.00    OERSTED  LANDE-G    FACTOR    =    2.0 


GENERAL    DATA    FOR 
DIELECTRIC     SUBSTRATE 


FERRITE    SUBSTRATE 
WAVE     TRAVELING 
FORWARD  I  REVERSE 


s  I  l-^b-da-  i  mmE  I 


C2MBDA" 

+ * 

O.C083 
0.0090 
0.0097 
0.0103 
0.0110 
0.01  17 
0.0123 
0.0130 
0.0137 
0.0143 
0.0150 
0.0157 
0.0163 
0.0170 
0.0177 
0.0183 


ZO 


tkm  i  mm  I  &m  i  mm 


2.500 
2.700 
2.900 
3.100 
3.300 
3.500 
3.700 
3.900 
4.100 
4.300 
4.500 
4.700 
4.900 
5.  100 
5.300 
5.500 


0.4043 

61.9 

0.122 

1.136 

0.182 

1.202 

0.4011 

62.6 

0.  180 

1.157 

0.265 

1.232 

0. 3981 

63.4 

0.278 

1.  186 

0.402 

1.271 

0.3952 

64.1 

0.454 

1.226 

0.649 

1.326 

0.3925 

64.7 

0.814 

1.288 

1.144 

1.409 

0.3900 

65.4 

1.688 

1.391 

2.336 

1.547 

0.3875 

66.0 

4.551 

1.582 

6.197 

1.  799 

0.3852 

66.6 

18.403 

1.7  29 

24.646 

1.984 

0.3830 

67.  1 

19.693 

0.2  93 

25.956 

0.077 

0.3309 

67.  7 

5.780 

0.  419 

7.494 

0.256 

0.3788 

68.2 

2.549 

0.602 

3.254 

0.503 

0.3769 

68.7 

1.458 

0.703 

1.832 

0.639 

0.3748 

69.  1 

0.966 

0.  763 

1.197 

0.721 

0.3731 

69.6 

0.701 

0.803 

0.855 

0.775 

0.3713 

70.0 

0.542 

0.331 

0.651 

0.813 

0.3696 

70.4 

0.437 

0.852 

C.518 

0.842 

TABLE  V.   Single  Slot  Line  on  a  Ferrite  Substrate 
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COUPLED    SLOT    LINES    ON    FERRITE    SUBSTRATE 

THICKNESS    OF    DIELECTRIC    =    O.IOO    CM         W/D    =    1.000  S/D    =0.500 

PERMITTIVITY    =    20.0 
SATURATION    MAGNETISATION    =    2300.00    GAUSS 

MAGNETIC    BIASING    FIELD    =       683.60    OERSTED 
LINEWIDTH    =       100.00    OERSTED  LANDE-G    FACTOR    =    2.0 


GENERAL    DATA    FOR 
DIELECTRIC    SUBSTRATE 


FERRITE     SUBSTRATE 
WAVE    TRAVELING 
FORWARD  |  REVERSE 


m  i  usBbds  i  mm1 


zo 


&m  i  mm  I  ?>m  I  mm 


EVEN 

2.900 

ODD 

MODE 

0.0097  j 
MODE 

0.4864 
0.3339 

59.0   j    0.185  j   1.151   j    0.251 
118.3       0.423  j   1.226   1    0.555 

1.209 
1.324 

EVEN 

3.100 

ODD 

MODE 

0.0103 
MODE 

0.4817 
0.3336 

60.1 
118.5 

0.307 
0.682 

1.186 
1.275 

0.410 
0.866 

1.254 
1.378 

EVEN 

3.300 

ODD 

MODE 

0.0110 

MODE 

0.4770 
0.3332 

61.2 
118.6 

0.557 
1.204 

1.2  40 

1.3  50 

0.734 
1  .485 

1.322 
1.462 

EVEN 

3.500 

ODD 

MODE 

0.0117 
MODE 

0.4726 
0.3329 

62.3 
118.8 

1.170 
2.467 

1.329 
1.476 

1.520 
2.951 

1.434 
1.601 

EVEN 

3.700 

ODD 

MODE 

0.0123 

MODE 

0.4683 
0.3326 

63.3 
119.0 

3.193 
6.565 

1.495 
1.708 

4.089 
7.640 

1.641 

1.858 

EVEN 

3.900 

ODD 

MODE 

0.0130 
MODE 

0.4642 
0.3322 

64.3 
119.2 

13.064 
26.226 

1.625 
1.882 

16.491 
29.736 

1.797 
2.036 

EVEN 

4.100 

ODD 

MODE 

0.0137 
MODE 

0.4602 
0.3319 

65.3 
119.4 

14.136 
27.744 

0.392 
0.  123 

17.609 
30.656 

0.252 

0.068 

EVEN 

4.300 

ODD 

MODE 

0.0143 
MODE 

0.4563 
0.3315 

66.3 
119.5 

4.193 

8.052 

0.497 
0.281 

5.160 
8.693 

0.391 
0.262 

EVEN 

4.500 

ODD 

MODE 

0.0150 
MODE 

0.4526 
0.3312 

67.2 
119.7 

1.868 
3.513 

0.654 
0.  506 

2.270 
3.710 

0.591 
0.518 

EVEN 

4.700 

ODD 

MODE 

0.0157 
MODE 

0.4490 
0.3308 

68.  1 
119.9 

1.079 
1.989 

0.741 
0.629 

1.296 
2.058 

0.701 
0.658 

EVEN 

4.900 

ODD 

MODE 

0.0163 
MODE 

0.4455 
0.3304 

69.0 
120.0 

0.722 
1.305 

0.793 
0.703 

0.858 
1.324 

0.767 
0.741 

EVEN 

5.100 

ODD 

MODE 
1   0.0170 
MODE 

0.4421 
0.3300 

69.8 
120.2 

0.529 
0.938 

0.828 
0.752 

0.621 
0.936 

0.812 
0.796 

TABLE  VI.   Coupled  Slot  Lines  on  a  Ferrite  Substrate 
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CO-PLANAR  WAVEGUIDE  ON  FERRITE  SUBSTRATE 


THICKNESS  OF  DIELECTRIC  =  O.IOO  CM    W/D  =  1.000 
PERMITTIVITY  =  20.0 
SATURATION  MAGNETISATION  =  2300.00  GAUSS 

MAGNETIC  BIASING  FIELD  =   683.60  OERSTED 
LINEWIDTH  =   100.00  OERSTED      LANDE-G  FACTOR  =  2.0 


S/D  =0.500 


GENERAL  DATA  FOR 
DIELECTRIC  SUBSTRATE 


FERRITE 
WAVE  T 


SUB 
RAVE 


STRATE 
LING 


FORWARD 

REVERSE 

TO 

LSM655       1     LAMBDA 

20 

mi  i  mim 

»tftt  !  IIH=6f 

2.500 

0.C083 

0.3345 

58.9 

0.  192 

1.  167 

0.269 

1.257 

2.  700 

0.0090 

0.3342 

59.0 

0.279 

1.  192 

0.378 

1.285 

2.900 

0.0097 

0.^339 

59.1 

0.423 

1.226 

0.555 

1.324 

3.100 

0.0103 

0.3336 

59.2 

0.682 

1.275 

0.866 

1.378 

3.300 

0.01  10 

0.3332- 

59.3 

1.204 

1.  350 

1.485 

1.462 

3.500 

0.01  17 

0.3329 

59.4 

2.467 

1.47b 

2.951 

1.601 

3.700 

0.0123 

0.3326 

59.5 

6.565 

1.708 

7.640 

1.858 

3.900 

0.0130 

0. 3322 

59.6 

26.226 

1.882 

29. 736 

2.036 

4.100 

0.0137 

0.3319 

59.7 

2  7. 744 

0.  123 

30.656 

0.068 

4.300 

0.0143 

0.3315 

59.8 

8.052 

0.281 

8.693 

0.262 

4.500 

0.0150 

0.3312 

59.9 

3.513 

0.506 

3.710 

0.518 

4.700 

0.0157 

0.3308 

59.9 

1.989 

0.629 

2.058 

0.658 

4.900 

0.0163 

0.3304 

60.0 

1.305 

0.  703 

I  .324 

0.741 

5.  100 

0.0170 

0.3300 

60.  1 

0.  938 

0.  752 

0.936 

0.796 

5.300 

0.0177 

0.3296 

60.2 

0.718 

0.  786 

0.705 

0.835 

5.500 

0.0183 

0. 3292 

60.3 

0.574 

0.8  11 

0.555 

0.363 

TABLE  VI.   Coplanar  Waveguide  on  a  Ferrite  Substrate 
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with  the  line  configuration  and  the  physical  parameters  so 
that  only  some  typical  examples  are  given  here.   To  generate 
the  Tables  I  through  VII  the  respective  computation  times 
are  60,  65,  95,  45,  55,  300  and  145  seconds. 

For  a  comparison  of  accuracy,  storage  and  computation 
time  the  Gaussian-quadrature  integration  was  replaced  by  the 
adaptive  Romberg  integration  routine  ' DCADRE ' .   The  total 
program  length  was  found  to  be  215  K  bytes  while  the  computa- 
tion time  increased  at  least  fivefold  to  guarantee  the  same 
numerical  accuracy.   The  detailed  and  involved  investigations 
of  the  numerical  methods  are  thus  justified  since  with  the 
here  presented  program  an  efficient  numerical  application  of 
the  theory  of  Section  2  is  presented  and  can  be  used  by 
future  users . 

5.3   COMPUTER  PROGRAM  • RESO '  FOR  RESONATING  SLOT  LINE  AND 

TRANSMISSION  LINE  PARAMETERS  FOR  SHIELDED  CONFIGURATIONS 

A  computer  program  which  provides  a  numerical  analysis 
of  shielded  slot  lines  and  halfwave  slot  resonators  is  shown 
in  Appendix  N  with  the  main  program  called  'RESO'  and  its 
accompanying  subroutines  ' COMBC '  and  'DET',  all  written  in 
the  FORTRAN  IV  programming  language.   The  preceeding  commen- 
tary section  of  the  main  program  describes  the  use,  the  neces- 
sary input  data  and  explains  some  of  the  variable  names  which 
are  important  for  input/output  operations . 

After  tne  input  data  is  read  the  program  computes  the 
lowest  order  mode  cut-off  frequency  for  a  rectangular  wave- 
guide mode  in  the  guide  formed  by  the  shielding  walls  and 
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the  metal  above  the  dielectric  substrate,  thus  including  the 
latter.   The  determining  equation  is  given  by  (3-1)  and  can 
be  easily  programmed,  shown  in  the  subroutine  'DET'.   The 
result  is  printed  and  is  subsequently  compared  with  the 
operating  frequency  where  any  further  computation  is  terminated 
if  the  operating  frequency  can  cause  an  onset  of  a  rectangular 
waveguide  mode.   Initialization  of  computation  loops  for 
varying  normalized  frequency  and  slot  width  follow. 

The  first  major  computation  is  to  find  the  wavelength 
ratio  Xl/X.      The  necessary  coefficients  for  the  determining 
series  are  provided  by  the  subroutine  'COMBC'.   Similar  to 
the  procedure  described  in  Section  5.2  a  bisectional  search 
to  find  the  zero  of  the  determining  functional  is  used  to 
improve  subsequently  the  accuracy  of  the  result  through 
several  iterations .   The  computation  stops  when  an  accuracy 
of  three  digits  after  the  decimal  point  is  guaranteed  for 
the  ratio  A/A' .   The  program  proceeds  with  the  computation 
of  the  characteristic  impedance.   The  coefficients  for  this 
computation  (equation  3-34)  are  provided  again  by  the 
subroutine  'COMBC ' . 

Finally  the  computation  of  a  halfwave  slot  resonator  is 
started.   Similar  to  the  iterative  procedure  for  the  disper- 
sion characteristic,  the  halfwave  resonating  length  is  found 
by  a  bisectional  zero  finding  procedure.   The  result  from 
the  dispersion  characteristic  is  used  to  ensure  that  the 
correct  zero  is  found  since  the  slot  resonator  would  not 


1G8 


only  resonate  at  a  halfwave  length  but  also  at  multiples 

thereof.   For  a  halfwave  length  resonator  the  normalized 

A  '  /A 
resonating  length  L/D  <  yn/<    •   A  first  guess  is  subsequently 

improved  until  a  one  digit  accuracy  after  the  decimal  point 
for  L/D  is  obtained.   The  necessary  integrands  for  equation 
(3-28)  are  provided  by  the  subroutine  ' COMBC ' .   It  was  men- 
tioned in  the  general  considerations  to  the  numerical  methods 
employed,  that  an  estimate  of  the  number  of  zeros  of  the 
integrands  could  improve  the  computational  efficiency.   By 
investigating  various  choices  of  input  parameters  it  was 
found,  however,  that  the  following  number  of  points  of 
Gaussian  quadrature  integration  schemes  guarantee  a  three 
digit  accuracy  for  the  integration: 


2  2 

y,   <  0,   Yo  <   °     12  points 


2  2 

Y-,   >  0,   y2  <   °     24  points 


2         2 

Y-i   >  0 ,   Yo   >  0     5(5  points 


After  finding  the  resonating  length  the  program  proceeds 
with  the  computation  of  the  end  effect  of  a  simple  short 
circuited  slot  line  by  evaluating  the  equations  for  the 
normalized  short  circuit  reactance.   All  results  are  printed 
together  with  some  of  the  input  data.   If  no  analysis  of  end 
coupled  slot  lines  is  desired  the  total  procedure  is  reiterated 
for  the  next  set  of  parameters  of  frequency  and  slot  width. 
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The  analysis  for  end  coupled  slot  lines  is  performed  in 
the  same  manner  as  the  single  slot  resonator  length  is 
computed.   Normalized  reactance  and  mutual  reactance  of  the 
coupling  region  are  derived  from  the  resonating  length  in 
the  odd  and  even  excitation  of  two  resonating  slots.   For  a 
given  operating  frequency  and  slot  width  these  parameters 
are  computed  for  the  desired  number  of  different  separations 
between  slots  and  are  printed.   The  program  starts  for 
another  set  of  input  data  from  the  beginning. 

The  subroutine  ' COMBC ■  provides  the  coefficients  and 
integrands  for  the  various  computations .   The  entry  variable 
is  3D  for  dispersion  and  impedance  computation  and  is  £D 
for  the  resonator  analysis.   All  variable  names  have  been 
chosen  such  that  the  coded  statements  can  be  readily  compared 
with  the  respective  equations  of  the  appropriate  appendices . 
The  transform  of  the  assumed  electric  field  distribution  is 
chosen  as  follows.   It  was  learned,  that  for  the  dispersion 
characteristic  and  resonator  analysis  a  distribution 


e  (1)   = 
xl 


l  +  |-4t|     |x|  <  w/2 

i w/2 '       ii—' 


(5-5) 

elsewhere 


which  is  mathematically  similar  to  (2-35a)  gives  almost  iden- 
tical results  as  if  a  distribution  according  to  (3-26a)  is 
chosen;  both  analyses  seem  to  be  relatively  insensitive  to 
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the  choice  of  basis  functions  while  the  characteristic 
impedance  computation  is  more  sensitive.   Furthermore, 
observe  that  the  largest  computation  time  in  this  program  is 
needed  for  the  iterative  procedures  for  dispersion  and 
resonator  analyses.   The  numerical  evaluation  of  the  trans- 
form of  (3-26a)  given  in  (3-27a)  requires  the  use  of 
another  subroutine  —  a  Bessel  function  routine  —  while  the 
transform  of  (5-5)  simply  involves  algebraic  manipulations 
of  trigonometric  functions  as  can  be  seen  in  equation  (2-35b) 
and  thus  is  faster  on  the  computer.   This  aspect  is  observed 
in  the  subroutine  'COMBC'  which  uses  the  Bessel  function 
for  the  power  density  computation  only. 

The  total  program  length  including  all  subroutines  is 
54  K  bytes  in  storage.   The  program  was  compiled  with  the 
FORTRAN  H  compiler.   A  sample  output  is  shown  in  Table  VIII, 
where  2  3  minutes  computation  was  needed. 
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SLOT-TRANSMISSION    LINE     AND    HALFWAVE    SLOT-RESONATOR 
ON    A    DIELECTRIC    SUBSTRATE    WITH    SHIELDING    WALLS 

DIMENSIONS:  A/D=    5.50  Hl/D=    8.00  H2/D=    8.00  EPSR=20.0 

THE    NORMALIZED    LOWEST    CUT-OFF    FREQUENCY    D/LC=0.0421 

W/D=0.200  D/L  =  C0200  LP/L=0.3549  Z0=    59.0  L/D=  8.1  XSC=0.129 

T/D=0.100  L-OOD/D=    8.5  L-EVEN/D=    7.8  XS=0.134  MS=0.113 

T/D=0.200  L-ODD/D=    8.4  L-EVEN/D=    7.9  XS=0.133  MS=0.100 

W/D=0.200  D/L=C0300  LP/L  =  0.3366  Z0=-    59.1  L/D=  5.0  XSC=0.173 

T/D=0.100  L-ODD/D=    5.3  L-EVEN/D=    4.7  XS=0.182  MS=0.16l 

T/D=0.200  L-OOD/D=    5.2  L-EVEN/D=    4.8  XS=0.179  MS=0.137 

W/D=0.200  D/L=0.0400  LP/L=0,3245  Z0  =    58.5  L/D=  3.4  XSC=0.239 

T/D=0.100  L-ODD/D=    3.7  L-EVEN/D=    3.2  XS=0.248  MS=0.234 

T/0=0.200  L-ODD/D=    3.7  L-EVEN/D=    3.2  XS=0.245  MS=0.200 

W/D=0.500  D/L=0.0200  LP/L=0.3684  Z0=    77.0  L/D=  8.4  XSC=0.148 

T/D=0.100  L-ODD/D=    8.7  L-EVEN/D=    7.9  XS=0.154  MS=0.142 

T/D=0.200  L-ODD/D=    8.7  L-EVEN/D=    8.0  XS=0.151  MS=0.128 

W/D=0.500  D/L=C0300  LP/L=0.3454  Z0=    77.2  L/D=  5.0  XSC  =  0.209 

T/D=0.100  L-ODD/D=    5.3  L~EVEN/D=    <t,6  XS=0.225  MS=0.209 

T/D--0.200  L-ODD/D=    5.3  L-EVEN/D=    4.7  XS=0.219  MS  =  0.182 

W/D=0.500  D/L=0.0400  LP/L=0.3299  Z0=    76.1  L/D=  3.4  XSC=0.295 

T/0=0.100  L~CDD/D=    3.7  L-EVEN/D=    3.0  XS^0.317  MS  =  C306 

T/D=0.200  L-ODO/D=    3.7  L-EVEN/D=    3.1  XS=0.311  MS=0.270 


TABLE  VIII.   Shielded  Slot  Line,  Half wave  Resonator  and 

End  Coupled  Slot  Lines  on  a  Dielectric  Substrate 
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6.   SUGGESTIONS  FOR  FURTHER  RESEARCH 

Although  an  extensive  analysis  of  slot  lines  and  coplanar 
strip  lines  has  been  presented,  several  problems  still  remain 
which  could  be  the  subject  of  further  investigation.   The 
intention  at  this  point  of  the  thesis  is  to  simply  list  a 
few  of  these  and  to  provide  some  guidelines  for  their 
possible  analysis. 

The  transforms  of  the  dyadic  Green's  function  components 
which  were  derived  here  are  general  and  could  be  used  in  the 
analysis  of  different  coplanar  configurations.   For  instance, 
Mariani  and  Agrios  presented  in  Reference  23  some  preliminary 
experimental  results  for  quarter  wave  coupled  slot  line 
resonators  in  a  shielding  box.   This  configuration  could  be 
useful  in  the  design  of  microwave  filters.   To  obtain 
numerical  results  for  this  type  of  line  all  one  has  to  do  is 
to  approximate  the  anticipated  electric  field  distribution.. 
This  approximation  has  to  be  Fourier  transformed  and  is  then 
applied  in  equation  (3-28)  where  it  replaces  the  transform 
of  the  assumed  electric  field  for  a  single  halfwave  slot 
line  resonator. 

Coplanar  transmission  lines  on  a  ferrite  substrate  have 
been  analyzed  in  this  thesis  by  a  perturbational  expression 
for  the  propagation  constant  only.   These  results  still  need 
an  experimental  verification.   Furthermore,  an  exact  analysis 
for  these  configurations  seems  highly  desirable  and  is  still 
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missing  in  the  literature.   It  seems  that  this  formulation 
can  be  easily  obtained  in  the  following  manner.   The 
transforms  of  the  dyadic  Green's  function  components  should 
be  derived  for  a  ferrite  substrate.   The  derivation  would 
follow  exactly  the  general  concept  shown  in  Appendix  A  where 
the  permeability  tensor  characteristics  would  be  implied  in 
a  general  formulation  from  the  start.   Although  this  method 
would  result  in  a  higher  degree  of  algebraic  complexity  it 
should  be  by  now  a  straightforward  procedure.   This  work 
would  be  more  general  and  complete  than  any  of  the  two 
analyses  reported  so  far  [Ref.  11,  this  thesis]. 

Circular  apertures  are  used  in  antenna  engineering  and 
their  analyses  are  described  by  various  authors  in  the 
literature.   But  nothing  has  been  reported  so  far  about  a 
circular  aperture  on  a  dielectric  substrate  and  hence 
specific  properties  are  unknown.   An  attempt  was  made  to 
apply  the  Fourier  transform  method  derived  in  this  thesis  to 
an  analysis  of  this  configuration.   The  steps  of  this  analysis 
are  briefly  outlined  in  Appendix  K  and  indicate  that  various 
and  difficult  problems  are  encountered.   It  cannot  be 
predicted  whether  this  formulation  is  feasible  at  all. 

A  final  suggestion  is  the  analysis  of  sandwich  slot  or 
coplanar  strip  lines.   These  line  configurations  will  be 
similar  to  the  one  described  in  this  thesis  where  another 
dielectric  slab  is  placed  on  top  of  the  conductor.   This 
configuration  was  first  suggested  by  Mariani  and  Agrios 
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[Ref .  23]  since  it  offers  the  advantages  of  shorter  wavelength 
and  greater  confinement  of  the  electromagnetic  fields  and 
thus  can  reduce  radiation  losses  or  perturbing  effects  from 
nearby  obstacles.   A  possible  analysis  would  be  started  by 
developing  the  transforms  of  the  dyadic  Green's  function 
components  for  this  geometry  where  perhaps  some  simplifica- 
tions could  be  made  due  to  the  symmetry  with  respect  to  the 
conductor  plane. 
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7.   CONCLUSION 

A  complete  hybrid  mode  formulation  has  been  presented 
for  the  analysis  of  coplanar  transmission  lines.   The 
analytical  methods  are  applied  in  the  Fourier  transform 
domain  where  dyadic  Green's  function  components  have  been 
derived  for  open  boundary  and  shielded  transmission  lines. 
These  components  are  very  general  and  do  not  depend  on  the 
specific  geometry  of  any  particular  line  configuration. 
Using  these  in  connection  with  Galerkin's  method,  exact 
equations  are  obtained  from  which  line  parameters  such  as 
the  length  of  a  resonator  and  dispersion  characteristics 
of  transmission  lines  can  be  computed.   Approximations  by 
the  use  of  a  finite  set  of  basis  functions  yield  good 
numerical  results. 

Using  the  results  of  the  dispersion  problem  the  theory 
was  extended  to  a  formulation  of  the  characteristic  impedance 
in  the  Fourier  transform  domain.   In  still  another  extension 
a  perturbational  expression  in  the  spectral  domain  was  applied 
for  the  analysis  of  the  propagation  coefficient  of  coplanar 
transmission  lines  on  a  ferrite  substrate.   Both  of  these 
applications  show  the  flexibility  of  the  formulation  in  the 
Fourier  domain  since  all  necessary  computations  can  be 
performed  in  this  domain.   Inversion  of  the  transforms,  a 
task  which  quite  often  is  difficult  to  accomplish,  is  thus 
avoided. 
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Some  design  curves  and  two  computer  programs  have  been 
derived  and  are  provided  within  this  thesis.   Various  new 
information  for  the  design  of  coplanar  transmission  lines 
is  therefore  available. 

Where  applicable,  results  were  compared  with  other 
existing  theories  and  good  agreement  was  found  in  all  cases 

Comparisons  with  experiments  were  also  performed.   In 
particular  the  numerical  line  parameters  for  coupled  slot 
lines  were  used  in  the  design  of  slot  line  directional 
couplers  of  which  several  different  types  were  investigated 
theoretically.   Due  to  the  dispersive  character  of  a  slot 
line  co-  and  contra-directional  coupler  design  is  feasible. 
Two  directional  couplers  were  built  and  their  actual 
performances  measured.   Good  agreement  between  theory  and 
experiment  further  verified  the  validity  of  the  analyses. 

In  this  thesis  particular  interest  was  devoted  to  the 
analysis  of  specific  coplanar  transmission  lines  such  as 
slot  lines  and  coplanar  strip  lines.   However,  the  method 
employed  is  completely  versatile  since  it  is  a  general 
formulation  of  an  inhcmogeneous  transmission  line  problem 
where  metal  conductors  of  any  shape  are  placed  upon  one 
side  of  a  dielectric  substrate. 
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APPENDIX  A 

Derivation  of  Spectral  Domain  Dyadic  Green's 
Function  Components  (Open  Boundary) 


Equations  (2-3)  and  (2-5)  are  substituted  into  the 
continuity  conditions  (2-14)  and  (2-15)  and  the  result  is 
transformed  into  the  Fourier  domain  which  yields 

at  y  =  0 


kj  »2(of0)  -  kj  #3(Of0) 

9  9 

•jay<!>2  (a,0)    -    joou  ^—  vi/2(af0)    =   -  jay<I>3  (a,  0)    -    jooy  ^--   $3(a,0) 

kc2    ^2(a'0)    =   kc0    ^3{ci'0)  (A_1) 

•jay4/2(a/0)    +    jwe2   t~   $2(a,0)    =   -jaY^3(a,0)    +   jwe0   |~-  $3(ar0) 


at  y   =   D 


kc0    *l(a'D)    "   kc2    *2<a'D> 

9  9 

jay^Cc^D)    -    jwu   -|-   ^1(a,D)    =   -jaY$2(a,D)    -    jwy   -^-  H'2(a,D) 

kc0    $i(a'D)    =    Ez(a) 
jay$1(a,D)    -    juy   ~  H'1(a,D)    =    Ex(a)  (A- 2) 

kcO¥l(a'D)    ~    kc2H,2(a'D)       =   Jx(a) 
■jaY^COfD)    +    ju)c0  |-   $1(a,D)    +   jcty^  (a,D)    -    jwe2  Jy  $2(a,D)    =   Jz 
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where  y-,  =  y?  =  y_  =  y  /   and   e,  =  e„  =  e~   so  that 


cl     c3     cO 


The  solutions  to  the  two  Helmholtz  equations  given  by 

equations   (2-12)  and  (2--13)  are  now  substituted;  the 

e  h 

dependence  of  the  coefficients  A   through  D   on  a  is  under- 
stood and  from  here  on  is  not  shown  explicitly. 


kc2  ^  -  kc0  »*  (A"3) 


jayC      -    jwyY2B      =   -JayD     -   jcoyY3D  (A-4) 


i     2    _h        ,     2    _  h  , ,.    r  v 

k    0    C      =   k    n    D  (A-5) 

Q.2.  CO 


h  eh  e 

jayC      +    jcoe2Y2B      =   -jayD      +    ju)e0Y3D  (A-6) 


_?   Ae   =   k_jBe  sinh  y2d  +   k    \^   cosh   y2d  (A-7) 


e  h  e  e 

jayA      +    jwyy^A      =   -jayB      sinh    y„D   -    DcofC      cosh   y2d 

-jcoyy2B      cosh    y2&   -    jwyY2C      sinh    y2& 

(A- 8) 

k   I   Ae   =    E„(a)  (A- 9) 

cO  z 


jaYAe   +    jwpYjA11   =    Ex(a)  (A-10) 
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k  2Ah  -  k  ?>Bh  sinh  y0D  -  k  2Ch  cosh  y0D  -  J  (a)    (A-ll) 
cO       c2         2      c2         '2      x 


jayA  -  jweny-,A   +  jayB   sinh  y2D  +  jayC   cosh  y2D 

6  e 

-joj£:2Y2B   cosh  y^D  -  jwe2Y2C  sinh  y2D  =  J  (a)     (A-12) 


The  objective  at  this  point  is  to  replace  the  unknown 

e  h 

coefficients  A  through  D   to  obtain  a  direct  relation  between 

the  current  density  and  electric  field  components. 

From  equation  (A-3) 

De  =  £c2, 2  ce 

kc0 


and  (A-5) 


Dh  =  (j^-)2  Ch 
cO 


Substitute  these  in  (A-4)  and  (A-6) .   Define  the  coefficients 
such  that 


F,  =  -^~    [(r^2-)2  -  1]  (A-13a) 

1    o)yy2    kcQ 


Y3    r»2  2 
F2  =  Z1    (v2^  (A-13b) 

^    Y2   Kc3 


then    Bh  =  F1Ce  +  F2Ch  (A-14a) 


E 

2 


and     Be  =  -5-  F0CG  -  -H-  F.Ch  (A-14b) 

e~   2      £_   1 
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Substitute  (A-14)  in  (A-7)  and  define 


G,  =  ix—-)       (—  F0  sinh  y9D  +  cosh  y0D)        (A-15a) 
1     Kc0     e2 

G9  =  -  (~i)  2  -H-  F,  sinh  y9D  (A-15b) 

-s       Kcl    e2   X       ^ 

then    Ae  =  G^6  +  G2Ch  (A-16a) 


Likewise  (A-8)  can  be  rewritten  if  the  following  definition 
is  used. 


G 


c0 


then 


_  =   aY   (G,  -  —  F~  sinh  y0D  -  cosh  \nD) 
3    wyy,    1    e2   2       '2         '2 

Y2 
-  —  F-,  cosh  y9D  (A-15c) 


G.    =  -2£-  (G   +  _LL  F   Sinh  YoD) 
4    wyY!    2    e2 

Y2 

-  —  (F9  cosh  y9D  +  sinh  YoD)  (A-15d) 

Y]_    £  £  £ 


Ah  =  G3Ce  +  G4Ch  (A- 16b) 


Now  solve  for  Be  through  C   in  terms  of  Ae  and  A  . 
Inversion  of  (A--16)  yields 

Ce  =  JxAe  +  J2Ah  (A-17a) 

Ch  =  J3Ae  +  J4Ah  (A-17b) 
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where 


G4  "G2 


1    det  G  u2    det  G 


-G  G 


3    det  G  4    det  G 


with  det  G  =  G1G4  -  G2G- 

Substitute  (A-17)  in  (A-14)  and  define 


1    e-   2  1    £         1  3 


G0 
2    e2   2  2    e    1  4 


K3  =  F^  +  F2J3 


K4  =  F  J2  +  F2J4 


then 


Be  =  KnAe  +  K2Ah  (A-19a) 


Bh  =  K3Ae  +  K4Ah  (A-19b) 


The  next  substitution  is  (A-17)  and  (A-19)  in  (A-ll)  and 
(A-12) . 
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Define 


L1  =  -(K_  sinh  y2D  +  J-  cosh  y2d) 

kc2  2 
L2  =  1  -  (^ — )  (K4  sinh  y2d  +  J.    cosh  y2D) 


OiY 

L3  =  T^l    (K3  Slnh  Y2D  +  J3  COsh  Y2D) 


E2Y2 

■g  (K1   cosh    Y2D  +   Jx    sin    Y2D)    -    1 


OtY 
L4    =  HT^TJ    (K4    sinh    Y2D  +    J4    COSh    Y2D   "    1} 


-    K2    cosh    y?D   -    J2    sinh    Y2D 


then 


kc2LlAe  +  kcOL2Ah  =  Jx(a)  (A-20a) 


jwe:0Y1L3Ae  +  jue2Y2L4Ah  =  Jz(a)  (A-20b) 


e      h 
Equations  (A-20)  may  be  solved  for  A   and  A   in  terms  of 

the  current  density  components.   The  result  substituted  in 

(A-9)  and  (A-10)  will  eliminate  all  unknown  coefficients  and 

express  the  electric  field  components  in  terms  of  the  current 

density  components  with  the  following  definition. 
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a,-  (  «,a  L  L  -  jLi  l2l3 

cl  2  2 


cO       cO£2Y2 


M2  =  (^?-L2  +  ^  <ft)2  V  m 


2'2 


2' 2    cl 


(A-21) 


M„  = 


DM 


M„  =  J 


kc0   L2 


we  y   DM 


Then 


Ml      M2 


M3      M4 


X 


X 


(A-22) 


If  the  M-matrix  is  inverted  another  set  of  coupled 
equations  is  obtained  in  the  form 


Nl      N2 


N3     N4 


z  J 


(A-23) 
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where  with 


det  M  =  M,M.  -  M0M. 
14     2  3 


M„  -M 


o 


M 

M    =       Z  N    =   - 

1    det  M  2    det  M 


-M  M, 


(A-24) 


3    det  M  4    det  M 


Note  that  the  use  of  hyperbolic  functions  applies  only 

2  2  .         . 

if   y2   >  0.   For   y«   <  0   the  trigonometric  solutions  for 

I — o T1 

$2  and  V-    aPPly  since   y~  =    jyj   for   Y2'  =\/k  ?  ~  a 
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APPENDIX  B 

Time  Average  Power  Flow  in  the  Spectral  Domain 

(Open  Boundary) 


The  following  expressions  for  the  time  average  power  flow 
in  the  three  spatial  regions  are  derived  from  equation  (2-26) 
The  solutions  to  Helmholtz's  equations  given  in  (2-12)  and 
(2-13)  are  substituted  in  the  transformed  scalar  potentials 
of  (2-27) .   Although  the  detailed  derivation  would  be  quite 
lengthy  and  therefore  prohibits  a  presentation  here,  the 
algebraic  manipulation  and  the  analytic  integration  with 

respect  to  y  is  routine.   In  the  following  derivation  the 

e  h 

dependence  of  the  coefficients  A  through  D   on  the  variable 

a  is  understood.   Note  also  that  y,  =  y~    since  regions  1  and 

3  are  assumed  air  filled. 

In  region  1 


2,   2 

1       °°      a  +Y1  ■  e.2       hi  2 

P    =  -  i-  Re  /   [3oo  —  (en  Aer  +  y  AT) 

avg     8  it        c     Y-i  0  '   '     K '   ' 


+  j2a(kQ2Ae*Ah  -  32AeAh*)]  da  (B-l) 
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2 

In   region    2    for    y2       >   0    ,      y2    real 


2+      2 

P  =   -     1    Re     /°°{{^   {sinh   2y9D[a    J2    (e,  (  1  Be  |  2+  \  Ce  |  2) 

avg  4tt  _oq       2  2  2y?         2    '       '        ' 

+  y(|Bh|2+|ch|2)) ] 


+  D(a2-Y22)[e2(|ce|2-|Be|2)  +  ]X  (  |  Ch  |  2-  |  Bh  |  2)  ] 

«2+Y92 

+  (cosh  2y9D  -  1) — o-=-(e 9  (BeCe*+Be  Ce) 
z        ^y2    z 

+  y(BhCh*+Bh*Ch))} 


+  j  |   sinh  2y2D[B2(BeCh*+Bh*Ce) 

-  k22(Be*Ch+BhCe*)] 

+  (cosh  2y2D-  1)  [B2(BeBh*+CeCh*) 

-k22(Be*Bh+Ce*Ch)]}}  da 

(B-2) 
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2 

In    region    2    for   y2      <    °    /      Y2    =    -3Y2    '      ^2    real 


00  v^      -a 


,„2         2 

Pavg  "   "  TF  Re  _{,   {1T  sin   2Yp  [-l^U,  ( |  Be  |  2"  i  Ce  |  2) 

+    p(|Bh|2-|Ch|2))] 
+    D(a2+y^2)  U2(|Be|2+|ce|2)    +    y(|Bh|2+|ch|2)] 


„2       2 

—ex 

+    2(1  -cos    2y^D)       2       „       [e2    ImfB^C")    +    y    ImfB^C11*)  ]  } 


Y      -a 
2  r        T    /C)e_e*x  T    /r,h_h*. 


2 


+   f   sin    2Y2,D[32(Bh*Ce-BeCh*)   -  k,,2  (Be*Ch~BhCe*)  ] 


+    j    |    (1-cos    2Yp)[32(BeBh*-CeCh*)   -k22(Be*Bh~Ce*Ch)]}    da 

(B-3) 


where  Im(  )  ■*  "imaginary  part  of  (  )" 


In  region  3 

2^   2 

1      °°     a  +Y1     ,  e.2     ,  hi  2 


+  j2a(B2DeDh*  -  kn2De*Dh)]  da       (B-4) 


The  minus  sign  in  front  of  all  these  integrals  appears  since 
a  wave  traveling  into  the  -z-direction  was  assumed. 


188 


For  slot  transmission  lines 

e  h 

Observe  that  the  coefficients  B   through  D   are  expressed 

e      h 
in  terms  of  A   and  A  with  equations  (A-19) ,  (A-17) ,  (A-5) 

and  (A- 3) . 

e      h 
But  A  and  A  are  expressed  with  (A-9)  and  (A-10)  in 

terms  of  the  transforms  of  the  electric  field  components  and 

are  therefore  known. 


For  copianar  strip  transmission  line 

e  h 

As  for  the  slot  line  the  coefficients  B   through  D   are 

e      h 
expressed  in  terms  of  A  and  A  .  (See  above.) 

Since  for  this  line  the  current  densities  are  approximated, 
an  explicit  dependence  of  these  two  coefficients  on  the 
transform  of  the  current  density  basis  functions  is  required. 

This  means  equations  (A-20)  in  Appendix  A  have  to  be 
inverted  which  yields 

Ae  =  -4-  J  (a)  +  — i  J  (a)  (B-5) 

k  ?   x      k?   z 

cl  cl 


L                     k  2    L 
A  =  -  — ^ J  (a)  -  :  ^ J z  (B-6) 


kcl  c2y2  DM  ue2y2  kcl 


DM 


where  all  coefficients  are  as  given  in  Appendix  A. 

With  this  modification  the  expressions  for  the  time 
average  power  (B-l)  to  (B-4)  apply  for  the  copianar  strip 
line  as  well  as  the  slot  line. 
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APPENDIX  C 
Perturbation  Theory  for  Open  Boundary  Waveguides 

In  this  appendix  a  perturbational  expression  for  the 
propagation  constant  of  a  wave  traveling  on  an  open  boundary 
waveguide  structure  is  obtained.   It  is  assumed  the  wave 
travels  into  the  z  direction  and  that  the  perturbation  is  a 
small  change  in  the  type  of  material  within  this  structure. 

-y  -y  , 

Let  E  and  H  represent  the  fields  in  the  unperturbed 

—  YZ 

structure  where  the  z  dependence  as  e    is  implied;  then 

v"  x  e  =  -jw[u]  H  (C-l) 

$  x  H  =  ju)[e]  E  (C-2) 

-  Y  '  Z 

and  for  the  perturbed  guide  (z  dependence  e     implied) 

V  x  e'  =  -jw[y+Ay]  H'  (C--3) 

V  x  h'  =  jw[e+Ae]  E1  (C-4) 

where  "[  ]"  denotes  a  tensor  quantity. 

The  last  two  equations  are  scalarly  multiplied  by  H*  and 
E*  while  the  conjugates  of  the  first  two  are  scalarly 
multiplied  by  H1  and  E1  respectively,  which  yields 
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ft'  .  v"  x  ft*  =   jw[y]*  3*  •  ft'  (C-5) 

ft'  .  $  x  ft*  =  -jw[e]*  ft*  •  E'  (C-6) 

ft*  ,  }x!'  =  -jw[y+Ay]  ft1  •  ft*  (C-7) 

g*  .  v  x  5'  =   jo)[e+Ae]  E'  •  E*  (C-8) 

Equations  (C-5)  through  (C-8)  are  now  combined  to  obtain 

$  .  (ft1  x  ft*)  =  ju>[e+Ae]  E»  •  ft*  -  jw[y]*  ft*  •  ft'      (C-9) 
?  •  (ft*  x  ft1)  =-jw[e]*E*  •  ft*    +  jco[p+Ay]  ft1  •  ft*    (C-10) 

where  use  of  a  vector  identity  is  made: 
namely     V*  (X"x  S)=S»?xX-X'?xS 

Now  add  (C-9)  and  (C-10)  and  integrate  over  some  volume  V 


Iff       [V  •  (ft*  x  ft*)  +  v"  •  (ft*  x  ft1)  ]  dv  = 
V 


ju>  ///  [[e+Ae]  ft'  •  ft  *  -    [el*  ft*  ■  ft']  dv 
V 


-jw  ///   [[y]*  ft*  •  ft'  -  [y+Ay]  ft'  •  ft*]'  dv      (C-ll) 
V 


The  restriction  is  now  imposed  that  the  material  of  the 
unperturbed  structure  is  isotropic  so  that 

[e]  =  e[l]  =  [c-]*   and  similarly  for   [y]  =  [y]* 
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With  this  assumption  some  simplifications  are  possible. 
A  further  simplification  of  (C-ll)  is  possible  through  an 
application  of  the  divergence  theorem.   Then 


<6/(H'  x  S*  +  H*  x  e')  •  n  ds  = 
S 


jco  ///   [[Ae]  •  E'  •  E*  +  [Ay]  H'  •  H*]  dv     (C-12) 
V 


This  equation  is  still  exact  and  may  be  applied  to  any 
region  V.  Choose  this  volume  such  that  a  differential 
cylinder  surrounds  it,  where  the  three  sides  are 


S1      z  =  a 


S2       z  =  b  =  a+dz 


S^      wall  of  cylinder  at  r  =  \j  x  +  y' 


=  V  J 


Note  that 


-t .    -»- .    ->■    -*■ .    ■*■        -*■ 
H'  x  E*  •  n  =  -H'  •  n  x  E' 


and     H*  x  e'  •  n  =  -H*  •  n  x  E' 


For  an  open  boundary  structure  and  r  very  large  the  outward 

-*■        -*■ 
unit  normal  vector   n  -»■  r 

But     lim  |H'  x  e*  +  H*  x  e'  I  •  r  •*•      0    since  the  fields 

finally  vanish  exponentially  for  large  r. 
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This  means  that  the  surface  integral  over  S_  is  equal  to 
zero.   Denote  the  surface  integration  over  S,  as  F(a)  and 
the  surface  integration  over  S~  as  F(b).   Then 


F(b)  =  -[F(a)  +  ±-   f(z) 

a  Z 


dz] 


Since   H1  x  E*  «  e  (V+Y*)z 


and 


H*  x  E'  oc  e 


(y'+y*)z 


then 


9z 


F(z) 


-(Y'+Y*)  F(a) 


z=a 


This  derivation  applies  for  any  a  along  the  z  axis. 
Equation  (C-12)  can  then  be  written  as 

-(Y,+Y*)3z  If       ($'  x  £*  +  ft*  x  FT1)  •  n  ds 
S 

=  jo)   ///   [[Ae]  £'  •  £*  +  [Ay]  ft'  •  ft*]  dv    (C-13) 
V 

where  S  is  the  cross-section  of  the  structure  and  extends 
to  infinity. 

The  right  hand  side  of  equation  (C-13)  remains  essentially 
constant  for  a  small  enough  variation  in  z  so  that  the  volume 
integral  can  be  replaced  by 


ju>  dz  // 
S 


[[Ae]  E"'  •  E*  +  [Ay]  ft'  •  ft*]  dS 
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Noting  that   [Ae]  =  [e1]  -  [e]   and   [Ay]  =  [y«]  -  [y] 


where    [e]  =  eQ([l]  +  lxel) 
and      [y]  =  yQ  (  [1]  +  [x]) 

the  difference  between  the  perturbed  and  unperturbed 
susceptibility  tensors  is  found  to  be 


[Axel  =  Ixll  "  IXe] 

e         t;         e 


[AX]  =  [x'l  -  M      . 


(C-13)  can  then  be  brought  into  the  final  form  of 


ju  //   (en[Axol  *  E«  •  E*  +  yn[AX]  •  H'  •  H*)  ds 

A°; 
Y'+y*  =  — — 

//   (E*  x  H'  +  E*  x  H*)  .""a   ds 
S  Z 


The  susceptibility  tensors  are  internal  tensors  and 
therefore  apply  to  the  internal  fields  only.   As  it  can  be 
seen  in  Section  2.4.2  it  is  more  convenient  to  use  the 
external  fields  and  external  tensors  if  due  to  the  physical  shape 
of  the  ferrite  demagnetization  must  be  taken  into  account. 

[Xe]  E'  -  [X^3  E 
and  [X]  H*  =  txG]  S 
where  the  superscript  "e"  denotes  external. 
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Equation  (2-40)  is  now  obtained.   This  detailed  perturba- 
tional  analysis  of  an  open  boundary  waveguide  could  not  be 
found  in  the  literature  and  hence  was  rigorously  carried 
out  to  prove  its  application  for  the  slot  line. 
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APPENDIX  D 
Perturbational  Expression  for  the  Propagation  Constant  in  the 

Spectral  Domain 

Apply  Parseval's  relation  to  equation  (2-42)  which  is 
then 


avn   —oo  f)     \  J  ■*  -* 


avg   -oo  0 


3  u;  *        3  (jr 


+  XZZkc2l*2|2}   d*  da  (D"1' 


The  transforms  of  the  scalar  potentials  are  known  through 

equations  (2-12)  and  (2-13)  and  the  dependence  of  the 

e  h 

coefficients  B   through  C   on  the  assumed  electric  field 

distribution  is  derived  in  Appendix  A.   The  integration  with 
respect  to  y  may  be  carried  out  analytically.   Since  this 

integration  is  performed  in  the  spatial  region  2,  the  two 

2  2 

different  solutions  in  this  region  for  y9   >  0  and  Y2   <  0 

have  to  be  considered.   The  amount  of  algebraic  manipulation 
is  quite  lengthy  to  obtain  the  final  form  of  the  integrand 
and  will  not  be  shown  here,  because  it  contains  no  new 
information.   Note  that  in  the  lossless  case  the  lefthand 
side  of  equation  (D-l)  is 
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Y«+Y*  =  a1  +  j$'  -  j3  =  a'  +  jA3 

A  shorthand  notation  for  the  final  form  of  equation  (D-l) 
is  then  given  by 


a'*   =  FTP2—       "   {{Cy[(6Y2)2(Bh)2   +    (ace   )2(C2)2 
avg     -oo    v      j  j 

±2we2aBY2BhCe] 


+   2x"    k J[±3Y2BhCh   +    awe2CeCh] 


kyz   c 

-U       -v. 

izz    c2 


+    X"    * o(Ch)2)     •    C 


+    (Xyy[(3Y2)2(Ch)2    +■    (aa)G2)2(Be)2    ±    2a3e2a3Y2BeCh] 


+   4,2Xyzkc2[±BY2Bhch   +    <™e2BeBh]    +    Xzzkc2(Bh)2}    *    Ci 


+    U-x"     [(BY2)2BhCh   +    (aa)£2)2BeCe 


±    u)e2aeY2(BeBh  +  CeCh)  ] 


+    X"    kc2[±PY2((Bh)2    +     4(Ch)2)    +    aoje2(4BeCh+BhCe)] 


+   XSzkjBhCh)      C3    }      d 


a 
(D-2) 


where  the  following  notation  is  used: 
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e  h 

The  coefficients  B   through  C   are  either  pure  real  or  pure 

imaginary;  if  pure  imaginary  take  the  real  part  (see  also 

Section  5) 


for  y2   >  0 


4  =  +1 


sinh  2y2D 
2 


ci =  —wr- +  5 


sinh  2y2D 
7 


C2     -5y7" 


cosh  2y2D  -  1 
:3  =    ZyZ  " 


2 

for  y2   <  0  ,   Y2  "  JY2    (replace  y2  by  y2'  in  (D-2) ) 


4  =  -1 


sin  2y»D    D 

ci    3y5 —  +  2" 


sin  2y"D 

c   =  —  - 

^2    2       Tylj 

1  -  cos  2y^D 

n      =    

^3         2yT| 

furthermore, 


Xjz  =  Re{^z} 
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If  the  expression  for  AB   is  desired,  the  following 
replacements  have  to  be  made: 


replace    xyy  by   xyy 


zz      Azz 


and        Y"    by   X1 
Ayz      Ayz 


The  superscript  "+"  on   a1  and   A£  denotes  that  these 
parameters  refer  to  a  wave  traveling  in  the  positive  z 
direction,  so  that  the  lower  signs  in  the  integrand  have  to 
be  taken  wherever  a  double  sign  appears;  the  superscript  "-" 
refers  similarly  to  a  wave  traveling  in  the  negative  z 
direction  where  upper  signs  in  the  integrand  have  to  be  taken 


199 


APPENDIX  E 

Derivation  of  Spectral  Domain  Dyadic  Green's  Function 
Components  (Closed  Boundary) 


The  continuity  conditions  given  in  equations  (2-14)  and 
(2-15)  are  transformed  via  (3-9)  into  the  two  dimensional 
Fourier  domain.   Transformed  field  quantities  can  be 
expressed  through  the  solutions  to  the  two  Helmholtz's 

equations  in  equations  (3-14) and  (3-15)  in  forms  of  the 

e  h 

coefficients  A  through  D  ,  whose  dependence  on  a  and  E,   is 

understood  from  here  on  and  is  not  shown  explicitly. 


at  y  =  0 


a   £Ce   -    jcoyY2B      =   ~an^De   sinh   Y-in2    ~   JwyYiD*     sinh   Yih2 


(E-l) 


Ce(k    2-£2)  =   De(k   2-£2)sinh   y,h0  (E-2) 


jwg2Y2B      ~    an^C      =    ^coeo^lD      cosh   "Yih2    ~    a    £D      cosh    Yih2 

(E-3) 
Ch(k22-C2)  =   Dh(kQ2-£2)    cosh   Ylh2  (E-4 


at  y  =  D 


e      .  h 

■an^A      sinh   Yt^-i    +   JtoVYiA      sinh   Yini    = 

-    anCB      sinh   y2d    "    an^cG   c°sh   Y2D    "    JWPY2B      cosh   Y2D 

-    ja)V!Y2C      sinh   Y2D  (E-5) 
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Ae(kQ2-£2)  sinh  y^  =  Be(k22-^2)  sinh  y2D  +  Ce(k22-az2)  cosh  y2D 

(E-6) 

-anCAesinh   y^   +    jcoyy-jA  sinh   y1h1   =    E^a^C)  (E-7) 


Ae(kQ2-C2)    sinh    Ylh1  =    Ez(anfO  (E-8 


•JcoEqY^6  cosh   Yihi    ~    an£A     cosh   Yihi   "   Jwe2Y2Be   cosh   ^2D 
Jwe2Y2Ce   sinh   y2D  +    an£B*     sinh   y2d   +    a    E.C     cosh   y2d 


=    V0^'^  (E~9) 


Ah(kQ2-£2)    cosh   Y^    ~   Bh(k22-£2)    sinh   y2D) 

-   Ch(k22-£2)    cosh   y2D  =    Jx(an, FJ  (E-10) 


Define 

2      2 

-r 

(k0"-^)    sinh   Ylh2 
and 


k2   -£' 
Fl    =    „     2      ? (E"lla) 


k   2-£2 

F2    =   2 \ "  (E-llb) 

(kQ    -£    )    cosh   Yj_n2 


then    (E-2)    and    (E-4)    can  be   rewritten   as 

De   =   F^6  (E-12a) 

Dh    =    F2Ch 
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e  h 

The  further  objective  is  to  express  B   through  C   in 

e      h 
terms  of  A   and  A  which  are  related  to  the  transforms  of 

the  electric  field  and  current  density  components  at  y  =  D. 

Substitute  (E-12)  into  (E-l)  and  (E-3)  and  use  the  following 

definition: 


Gl  =  an^Fl  sinh  Ylh2  "  1} 


G2  =  jwyY1F2  sinn  Y1^2 


G3  =  ^WG0Y1F1  COSh  Ylh2 


'4  =  ~an^F2  cosh  Y1n2  ~  U 


(E-13) 


Then 


G1Ce  +  G2Ch  =  ju)uy2Bh 


G3Ce   +   G4Ch   =    jco£2Y2Be 


(E-14) 


(E-14)  can  be  solved  for  Ce  and  C   in  terms  of  Be  and  B 
to  yield 

Ce  =  HxBe  +  H2Bh  (E-16a) 


Ch  =  H3Be  +  H4Bh 
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with 


H,    = 


-ju)G2Y2G2 


1  det   G 


jtoyY2G4 
H2   =      det   G 

jcoe9Y9G 
H3   =      det  g  (E"15) 


H,    = 


-jo)yY2G3 


4  det   G 


det   G   =   G1G4    -    G2G^    . 


With  (E-16)  nov7  substituted  into  (E-5)  and  (E-6)  and  the 
definition  of 

J-,  =  -a  £ (sinh  YoD  +  Hi  cosh  Y?D)  ~  JWVY9H3  s*-nh   ^2D 

J2  =  -jtoyY2^c°sh  Y?0  +  H4  sinn  Y2D)  "  an?H9  c°sn  Y2D 

(E-17) 

2   2 

J-  =  (k9  -E,    )  (sinh  YoD  +  Hi  cosh  y?d) 


2   2 

J4  =  H2 (k2  -£  )  cosh  y2d 


these  are  rewritten  as 

J1Be  +  J2B  =  -Aea  E,    sinh  Y-in-i  +  JA  WPY]_  sinh  y-^1 

J3Be  +  J4Bh  =   Ae(kQ2-^2)  sinh  Yxn! 
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eh  e 

where  now  B   and  B   can  be  solved  for  in  terms  of  A  and 


Ah. 


Define 

[-an£J4  -  (kQ2-C2)J2]  sinh  y1h1 
Kl  =  3eT"J 

jwyY1J4  sinh  y1h]L 
K2  =      det  J 

[0^3  +  (kp2-^2)^]  sinh  y1h1  (E-18) 

3  det  J 

- j cay YyJo sinh  Ynh-i 
K4   =  det   J 


det   J  =  J, J.    -   J2^3 


Then 


Be  =  ^A6  +  K2Ah  (E-19a) 

Bh  =  K3Ae  +  K4Ah   .  (E-19b) 


Now  (E-16)  and  (E-19)  are  substituted  into  (E-9)  and  (E-10) . 
Define 

L1  =  -jtoeQYjL  cos  Y1h1  -•  ju)e2Y2K1  cosh  y2d  _  Jwe2Y2  ^H1K ]+H2K3^  ' 

•sinh  YoD 


+  an^K3  sinh  y2d  +  an^H3Kl+H4K 3^  cosh  Y2D 
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L2    =   -ctnC   cosh    Y1h1    -    jcoe2Y2K2    cosh   Y2D   "    J we2Y2  ^H1K2+H2K4^  sinh   Y2D 
+    anCK4    sinh    y2D   +    an£  (H^+H^)    cosh    y2d 


(E-2  0) 


L3  =  -(k22-^2)K3  sinh  Y2D  "  (k22~^2) (H3K1+H4K3)  cosh  y2d 


L4  =  (k02-?2)cosh  Y1h1  -  (k ^ - E,2 )  K 4  sinh  y2d 
-  (k22-^2) (H3K2+H4K4)  cosh  y2D   . 


Then 


L1aG  +  L2Ah  =  J  (E-21a) 


L3Ae  +  L4Ah  =  Jx   .  (E-21b) 


e       h 
Equations  (E-7)  and  (E-8)  can  be  solved  for  A  and  A 

yielding 


E 
Ae  =  - 2 - —  (E-22a) 

(kQ  -£  )  sinh  y^ 

,  E  a  £E 

An  =  -, — * _  + n   *  , .   (E-22b) 

iuipYn  Sinn  Ynb,    .     ,,  2  r2>   .  .    . 
J   '1       '11    :coyY1(kQ  -5  )  sinh  y1h1 


Substitute  (E-22)  in  (E-21)  and  define 


La 

3   Jwpy][  sinh  y^\ 
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Ll an^L2 

M„    =   ^ « + 


4  2      2  2      2 

(kQ    -£    )    sinh  Y1h1        j(jjyy1(k0    -£    )    sinh  y^^ 


L.  (E-22) 

1        Dooyy,    sinn   Yt"7 


m     =  3 n_    4 

2        Z    2   rs:    ...   .       ~  +  :  #1    2   71 


(kQ   -£    )    sinh   Y1h1        jwyY1(k0    -C    )    sinh   Y]_h1 


Then 


Mt  E      +   MnE      =   J  (E-23a) 

1   x  2    z  x 


M3Ex   +   M4Ez    =    Jz       '  (E-23b) 
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APPENDIX  F 

Time  Average  Power  Flow  in  the  Spectral  Domain 

(Closed  Boundary) 


A  short  hand  notation  for  equation  (3-34)  is 

oo 

P    =  —  Re    I         (P,   +  P0   +  P0  )  (F-l) 

avg    4a         „,    In     2n     3n 

where  the  first  subscript  denotes  the  spatiol  region  1  to  3 . 

These  terms  P.   are  then  found  by  substitution  of  (3-17) 

in  J 

and  (3-18)  into  (3-34)  where  the  variable  £  is  replaced  by 
the  constant  £.   The  integration  with  respect  to  y  can  be 
performed  analytically  since  (3-17)  and  (3-18)  show  a  simple 
dependence  on  this  variable.   Without  showing  the  details 
of  this  rather  lengthy  but  straightforward  algebraic 
manipulation  the  single  terms  for  equation  (F-l)  are  given 
below  with  the  following  implication.   Only  a  first  order 
approximation  will  be  used,  hence  E    [a    )    =0   and  through 
(E-22a)   Ae(an)  =  0  . 
Then 


2    2 
Pln  =  -T-  f  "TyT"  sinh(2y1h1)  +  (a^-y,2)  h,]  |  Ah  (%)  |2 


(F-2) 
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a  2^        2 

P3n  =  -T-  !  -JL^T  Sinh(2Ylh2'  "  <a„  -V»h2]  l^'Vl2 

2    2 

wyn6    a  +y,                   0    9  .      9 

+  _JL_  [  -n2^  -1-   sinh(2Ylh2)  +  (an  "Y!  )h2]  |Dn(an)r 

+  an(kQ2+S2)  De(an)  Dh(an)  sinh(2Y]Lh2)  (F-3) 


For  P0   it  is  not  necessary  to  go  through  this  cumbersome 
algebraic  work  if  one  observes  that  equations  (3-17)  and 
(3-18)  are  written  according  to  the  nomenclature  used  for 
the  solutions  for  the  open  boundary  Helmholtz's  equations  in 
the  spectral  domain  given  in  equations  (2-12)  and  (2-13). 
Therefore  the  derivations  for  region  2  in  Appendix  B  hold  if 
only  the  integrand  is  considered  in  equations  (B-2)  and  (B-3) 

and  a  is  replaced  by  a  . 

2 
So  for  y2  >    0 


P2   given  by  integrand  of  (B-2)  (F-4) 

where  a   replaces  a 

n    c 

2 

and  for   y2   <  0 


P2   given  by  integrand  of  (B-3)  (F-5) 

where  a   replaces  a 
n 


h  h 

The  coefficients  A   through  D*  for  the  above  equations 

are  related  to  the  assumed  transform  of  the  electric  field 

component  E   as  given  in  Appendix  E  with  coefficients  F-,  to 
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K4  and  equation  (E-22) .  All  terms  of  (F-l)  are  then  known 
and  can  be  evaluated  on  the  computer  after  a  normalization 
technique  has  been  applied  which  is  described  in  Section  5 
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APPENDIX  G 
Normalized  Equations  for  Open  Boundary  Lines 

In  this  appendix  the  equations  of  Appendix  A  leading  to 
the  dyadic  Green's  function  components  are  rewritten  such 
that  quantities  are  scaled  with  respect  to  the  thickness  D 

of  the  dielectric  substrate.   The  equations  are  first 

2  2 

written  out  for  the  two  cases  where   y2   >  0   and   y?   <  0 

and  are  then  combined  in  a  form  which  lends  itself  towards 

a  computer  programming  application. 

In  the  following  notation  an  imaginary  quantity  u  is 

written  as  ju"  where  u"  is  then  real.   The  normalization 

of  the  frequency  with  respect  to  D  leads  to 


n.  1  D 
we0D  "  6ff  T 


£r  D 


we:2D  =  601 


wy0D  =  60  j    (2tt)  2 


where  now  D/X  is  the  normalized  frequency  and 


or*    0   D  X 

BD  =  27T  j    p- 
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From  equations    (2-1)    and    (2-2)    then 


(kclD)2    =    -    (2tt   £)2    [(^r)2  -1] 


(kc2D)2   =    (2,  f)2    [er   -    (^)2] 


2  2  ^ 

Y-.D     =  [(aD)   -  (k  ,D)  ]  2  which  is  a  real  quantity 

for  aD  since   1  <  yt  <  e 

—  A   —  *• 


2 
a)       Case  1,  y?   >  0 

According  to  equation  (2-11)  the  above  condition  is 

equivalent  to   (aD) 2  >  (2tt  £) 2  [e„  -  (A-)2]  so  that 

A      r     a 


Y2D  =  [(aD)2  -  (k  2D)2^   is  real. 


The  coefficients  F,  through  N.  from  Appendix  A  are  then 


F,  =  j  F" 


p"  _  _aD_  JD_  rfkc2D  2    ., 
Fl    toyD  y2D  l(kclD;     1J 

YlD   k   D 

r2    YoD  lk  ,D; 
1 2     cl 

kc2D  2   F2 

Gl  =  (i^D)   (f   Sinh  T2D  +  COSh  Y2D) 
cl      r 


G2  =  j  GJ 
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cl       r 


Go  =  J  G  '• 


G3 


-S^r  -££r  (G,  -  -i  sinh  YoD  -  cosh  y.D) 


YoD 


Y 


D  Fl  COSh  ^2D 


G .  -  -  J*  M-    (G;  +  (120ir)   F-  sinh  y9D) 
4     wyD  YiD   2      er    -L       '2 


YoD 


Y 


^=-  (Fn  cosh  YoD  +  sinh  Y2D) 


det  G  =  G1G4  -  G2G3  =  G  G4  +   GJjG" 


Jl  =  G4/det  G 


J9  =  j  J" 


J2  =  ~  G2^det  G 


J3  -  *  ^3 


J3  =  "  G3/det  G 


J.  = 


K,  = 


K„  = 


Gx/det  G 
1 


K2 


(F2J1  +  (120-rr)"  FJ  J£) 


j  K« 

~  (F2J2'  -  (120tt)2  FJ  J4) 


K. 


j  K^ 
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K4   "  "   Fi   J2   +   F2J4 


Ll  ■  i  Li 

Ll  =  ~  *K3  sinh  Y2D  +  J3  cosh  Y2D) 


k  2D 

L2  =  1  ~  j~— 5  (K4  sinh  y2D  +  J4  cosh  y9D) 
cl 


L,  = ~r   -£~  (k»  sinh  y0D  +  J'4  cosh  y~D) 

3  coSqD   y-^D         3  '2  3  '2 

-    er  — -^    (K1   cosh   y2D   +   ^    sinh   y2D)    -    1 


L4    =    J    K 


la   =  -^n  ~^K  <K*  sinh  Y0D  +  J  a    cosh  y0D  -  1) 
4    coe^D  y^D   4       z  4       '2 


-  (K2'  cosh  y2D  +  J"  sinh  y2D) 


kc2D  2  1  Y1D 

DM  =  -  <«    LiLJ-r^L2  L3 

cl  r  '2 


M1  =  j  MJ 

„  _  r  aD  3D    ,  _   "V   ,¥,2      1 

1  "  [  7T^  "4    S  T2D  (1TID)   L3]  DM 

M2  =  j  M2' 

M«  r  aD  3D    t   .  (120tt)2  Y1D  ,kc2D,  2   „  .   1 

M2  l  We2D  y2D  L2  +    ^T~"  y2D  ^k~[DJ    1J  DM 
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M3  =  j  M" 
M'^    =  LJ/DM 

M4  =  j  M4 


(k_nD)2  L. 


„.  _     cl_    2 


M"  = 


4    we  D  Y2D  DM 


det  M  =  M£  M"  -  MJ  M]J 


Nx  =  j  NJ 

NJ  =  MJ/det  M 


N2  =  j  N^ 


N^  =  -M^/det  M 
N,  =  j  N" 


NjJ  =  -M^/det  M 


N,  =  j  N" 


N£  =  Mj/det  M   . 

b)   Case  2,  y22  <  ° 

This  condition  implies  that  the  normalized  transform 

variable   (aD)2<  ( 2 7r  ?)  2  [e   -  (-A-)2]  so  that 

A      e     A 

Y,  =  j  Yo 


Y2  =  C(kc2D)2  ~  {aD)?'^    ■ 


214 


The  coefficients  F,  through  N,  are  then 


_  ocD   BD  r  ,kc2D.  2    n. 


F   =  -i  F»» 
r2    J   2 


Y1D   kc2D  2 

F2  "    yTd  (k~D) 
1 2     cl 

k  9D       F" 
Gl  =  {k£^)   (~  e   sin  Y2D  +  cos  Y2D) 


G0  =  j  G" 


G„=.(W)2U2^PiSinY„D 
cl        r 


G3  =  1  G3 

F"  y"D 

G3  "  ^D  7^5    (G1  +  4  Sin  Y2D  "  COS  Y2D>  "  yjD  Fl  C°S  Y2D 

G   =  -  _^2_  JD        (120tt)2  F   sinY?D) 
4      a)yQD  Y-j_D  '2       er     1      T2 

Y2'D 
+  ~~  (FJ  cos  y2'd  +  sin  Y2D) 


det  G  =  GXG4  +  G^G'^ 


Jl  =  G4/det  G 


J2  "  i   J2 
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J"   =   -GJ/det  G 


J3  "  i   J3 

J"  =  -G"/det  G 

J4  =  G1/det  G 


K,  =  j  K" 


Ki  =  T~  (F2  Ji  "   (1207T>2  Fi  J3> 


K2  =  F"  {"F2  J2  "  (1207r>2  Fi  J4) 
r 


K3  "  F1J1  "  F2J3 


K^  ■  j  K« 


KJ  =  FXJ»  +  F2J4 


Ll  "  3  LJ 


L"  =  -  (K3  sin  Y2D  +  J3'  cos  Y2D) 


k  ?D  ? 
L2  =  1  -  (j—^)   (-KJ  sin  y2'd  +  J4  cos  Y2D) 
cl 


L~  = ^~  -£~    (K,  sin  YoD  +  Jo  cos  Yo^) 

3      ooe^D  y-id    3      '2      3      '2 

Y2D 
+  er  -~  (KJ  cos  y2'd  +  Ji  Sln  Y2D)  "  -1 


L4  -  3ipj  ^5  ("K4  sin  Yp  +  J4  cos  Y-n  -  1) 


-  (K2  cos  y2'd  ~  J2  sin  Y5D) 
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DM  =  j  DM" 


kc2D  2  1  Y1D 

DM«  =  -  (_£i_)2  L,L4  _  A_^L2  L3 


M1    =  j  MJ 


1     <k  .D)2   4    £r  ^2D  <kclD    3'  DM" 

Cl 
M9  =  j  M" 


„  =  .  ap  BD        (120tt)2  Y1D  fkc2D,  2   „ ,   1 
m2    1ojc2D  y^D   2      er     y^D   kc]_D;    1J  DM" 


M3    =    j    M" 


M"    =   L4/DM" 


M,    =    j    M" 

4         ->       4 


(k    ,D)2     L. 


Lcl  2 


M"     = 
4         co£2Dy^D      DM17 

The  equations  for  N,  through  N.  are  identical  with  those 

2 
listed  above  for   y-   >  0   and  need  not  be  reiterated  here 

For  programming  purposes  it  is  now  important  to  notice 

that  indeed  a  complex  arithmetic  on  the  computer  can  be 

avoided  since  only  a  real  arithmetic  is  necessary  for  the 

evaluation  of  the  dyadic  Green's  function  components  M, 

through  N..   A  further  simplification  is  even  possible  by 
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comparing  the  two  sets  of  equations  for  the  two  cases 

2  2 

Y2   >  0   and  y~      <  0   where  it  can  be  seen  that  some  of 

the  equations  differ  only  by  the  hyperbolic  or  trigonometric 
sine  and  cosine  function  and  some  algebraic  signs;  the  double 
prime  superscript  is  immaterial  for  the  further  computation 
since  it  is  only  used  to  distinguish  between  real  and 
imaginary  quantities  in  the  above  derivation.   It  is  then 

possible  to  reduce  the  two  sets  of  equations  above  to  one 

2  • 

set  where  depending  on  whether  y2   1S  greater  or  less  than 

zero  the  proper  signs  and  functions  are  applied. 

The  procedure  is  self  evident  and  needs  no  further 
development  in  this  appendix.   The  result  can  easily  be  seen 
j.n  the  subroutine  named  "COMBO"  of  the  attached  computer 
program   "COPLAN"  where  a  consistent  notation  with  the  one 
used  above  is  found. 
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APPENDIX  H 
Normalized  Power  Flow  Equations 

The  equations  for  the  time  average  power  flow  of  open 
boundary  coplanar  transmission  lines  are  given  in  Appendix  B. 
To  accomplish  a  formulation  which  allows  a  numerical  evalua- 
tion on  the  digital  computer  the  following  steps  are  to  be 

taken. 

e  h 

The  coefficients  A  through  D  depend  on  the  assumed 

field  distribution  at  y  =  D  as  stated  in  Appendix  A.   These 

2 

dependences  are  given  for  the  two  cases  where  y?   ^s  greater 

or  less  than  zero  as  follows  in  normalized  form  for  slot  lines: 


Y22  >  0 

Ah 

h" 
.  An 

Ah" 

E  (aD) 

D2 

wyD  YtD 

Be 

Ah" 
=  "  K2  ^2" 

Bh 

7 

h" 
.  B 

Bh"    ,   Ah" 


k4   2 

D"        D 
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? 


=  -  J 


„  A 


h" 


2  — 5" 


Ch 


.h" 


=  J 


Ah" 
4  "?" 


D 


,kc2D,2  Ce 

(FTd)  IT 

cl     D 


Dh 
7 


.  Dh" 

3  ir 


D 


h" 


(kc2D  2  Ch" 


D 


cl 


D 


and  for  y2   < 


as  above  except 


.  Bv 


D' 


Bc     .  B 
D        D^ 


B* 


=  K 


D 


Ah" 
2"^ 


B 

7 


=  -  K 


„  A 


h" 


4^ 
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With  the  assumption   E  (x)  z    0   it  follows  that   A   =  0. 
The  normalized  equations  for  the  power  flow  then  become. 


For  region  1 

.      -          (aD)2  +  (Y,D)2  h" 

Pl  -  -  h       !  ^   WD  p-^ <V   d(0D)      (H-l) 

-co                      '1  D 


2 

For  region  2 ,  y2   >  0 


k_„D 
P 


1      C    [  BD   Sinh  2y2D   (aD)  +(^2D) 

=  "  **  -k  'd   I  T  {    2     [    V5 


2 

~"kc2 


e      e  h"      h" 

.(we  2D((-T)2  +  (^-)2)  +  o)yD((^-r)2  +  (£-2_)2))] 
D       D  D        D 


e       e  h"      h" 

+  [(aD)2  -  (y9D)2]  [coc9D((^-)2  -  (-T)2)  +  coyD((^-T-)2  -(*   )2)] 
^        ''    D^       D  D       DZ 


(aD)2  +  (y2D)2       Be  ce      Bh"  ch" 

+  (cosh  2y0D  -  1) =r (coe0D  — *-  —*■  +  coyD  — ~ -7-)  ) 

2  y2D  2   D2  D2       D2   D2 


sinh  2y9D   Re   h"    pe   h" 

+  aD[(3D)<f  (k9D)2]  [ --i-  (Sy  --*-  +  —■  S_.) 

^      D^  DZ    D   D 

cosh  2y9D  -  1   Re  nh"    „e  nh"   > 
+  - J- (B   B    +C   C   }]    d(aD)       (H_2a) 

Z        D^  D     DZ  DZ   ) 


221 


Fo 


r   region   2 ,    y„      < 


P2    = 


kc2D      )  f    _      sin    2Yp         (aD)  2- (y2'D)  2 


TtF 


f 
-oo  (kD) 


BD 


{■ 


[ 


Yp" 


e  e"  h 

(we2D((^-)2  -  (^T-)2)    +   a3pD((^-)2 
D  D  D 


D 


e"  e  h  h" 

+    ((aD)2+(Y^D)2)  [o3e0D((^-T-)2  +  (^-)2)    +   coyD  (  (2_)  2  +  (<L_ )  2)  ] 


DZ  D 


DZ  D 


+    (cos    2YoD  -  1) 


(aD)2-(yp)2 


e"      f» 
Be      Ce 


Bh   Ch" 


-^g-  (aje:2D  -^  -^        wyD   -j  -^-J  } 


D       D 


D        D 


2^ 


Bin    2Y2-D        h   ce        Be' 


+   aD[($D)^+(k9D)^]  [- ^~   (Z—  ^  -  i^-  ^-) 

Z  DZ    D^  D         D 


cos    2Y^D  -  1      ne"    _h        ne  rh" 


(^-7-  —  +—^-)]     f      d(aD) 
D       D^         D^     D 


(H-2b) 


For-  region    3 


P3    = 


1      ".kc2D      ,kc2D.2    fBD    <«°>2+<Yi°>2 

_  oo  cl  I.  1  I 


,h" 


e  0h" 


[a)e1D(^T)2    +    a)pD(^-T-)2]    +    aD[(3D)2+(knD)2]    £_  £-_    I    d(aD) 
1      D  D  D        DZ     J 


(H-3) 
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All  coefficients  in  these  integrands  are  real  quantities 
now  and  can  be  easily  programmed  into  a  digital  computer 
language. 

To  accomplish  an  evaluation  of  the  time  average  power 
flow  for  the  coplanar  strip  line  with  real  coefficients, 

equations  (B-5)  and  (B-6)  in  Appendix  B  have  to  be  written 

2  2 

explicitly  for  the  two  cases   y?  >   0   and   y«   <  0. 

Note  that  in  the  equations  given  here   j  (x)  ~  0   is 

assumed.   Furthermore,  from  equations  (2-33)  to  (2-35)  it 

is  seen  that  the  transform  of  the  assumed  current  density 

J  (a)   is  imaginary  so  that 
z 


Jz(a)  =  j  J£(a) 


Following  the  guidelines  above  and  those  in  Appendix  G, 
the  coefficients  Ae  to  D   are  then  derived  in  terms  of  the 


assumed  current  density. 


For  Y,2   > 

0 

Ae 

M" 
4 

Ah 

(kclD)2 
.    Ah" 

J" 

z 
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h"  L"  k    0D    0    J" 

{        _  1 ,    c2    .  2      z 

D  2         2  cl 


e  e  h" 

"T     ~     1  ~T         2  ~~T" 


Bh  .    Bh" 

-j      -    j    -j. 

D  D 


h"  e                 h" 

ff —     Tf"  J-     V         — 

T2T  '    K3  "T         4   TT 

D  D^                   D 


Ce  T      Ae 

~T  1  T 

DZ  X    D 


Ah" 

J2  ~~r 

D 


ch        ,  ch" 


h"  e                  h" 

C  A                        A 

"TT  ~  J3   T       J4 

D  J    DZ                   D 


D        _    ,   c2    ,2   C 
D2  clD        D^ 


Dh 


=  : 


D 


h" 


D 


hM 


,kc2»2    Ch" 
cl  D 
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For  y_   <  0 


Ae  _            MJ 

JZ 

D2              (kclD)2 
Ah          .    Ah" 

h" 


D 


LI'     k  ~D  0  J" 
,  c2  > 2   z 


we2D  Y^D   k  ,D'   DM* 


5* 


B 


e" 


D 


e"  e        h" 

B  A         A 

A    —   X\,-j    X-  T  1\~        tj— 

DZ  X  DZ            Z  D 


h        e       h" 
B        A       A 


and  the  rest  of  the  coefficients  are  as  given  above  for 

Y22  >  0  . 

The  normalized  equations  for  the  power  flow  are  then, 
in  region  1 


Pl  " 


1 

Tjt" 


!°/BD  (oD)  +  (Y-.D)2        e  ,       -h"  , 
--  ^   '   Yid       I^0D(^)2  +  oviDC^)2] 


e  h" 
+  aD[(3D)2  +  (k  D)2]  $j  ^T- 

D   D 


d(aD) 


H-4) 
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For  regions  2  and  3  equations  (H-2)  and  (H-3)  apply 
without  any  further  change  so  that  here  again  a  formulation 
with  real  quantities  only  can  be  achieved. 

The  further  translation  into  a  computer  language  is  then 
routine. 
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APPENDIX  J 
Normalized  Equations  for  Closed  Boundary  Lines 

Similar  to  the  detailed  manipulation  of  the  various 

coefficients  for  the  open  boundary  structures  shown  in 

Appendix  G,  a  complex  arithmetic  for  the  closed  boundary 

coplanar  lines  can  be  avoided.   The  equations  leading  to  the 

transforms  of  the  dyadic  Green's  function  components  are 

derived  in  Appendix  E  in  a  general  form.   It  must  be  observed 

that  the  smallest  value  the  discrete  transform  variable  a 

n 

can  attain  is  zero.   For  the  computation  of  the  dispersion 
characteristic  and  subsequently  for  the  line  average  power 
flow  the  transform  variable  £  is  replaced  by  the  propagation 
constant  3  as  stated  in  Section  3.2  where   k,  <  $  <  k~  , 

so  that  for  computations  concerning  these  two  line  parameters 

2  2  .  . 

y,   will  always  be  positive  but  y?   may  be  positive  or  negative 

depending  on  the  value  of  a  . 

In  the  computation  for  the  resonant  slot  length,  however, 

the  transform  variable  £  extends  such  that   -°°  <  E,    <   °°  , 

thus  in  general  three  possible  solutions  for  the  transformed 

Helmholtz's  equations  (3-15),  (3-16)  must  be  considered, 

namely 

a)  Y]L2  >  0    ,    y22  >  ° 

b)  Yt2  >  0    ,    Yn2  <  0 


x  ,  i2 

2 

i    <  °     '     v 


c)  y,2   <   0   ,   Yo2  <  0 
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where  a)  and  b)  apply  for  the  dispersion  and  characteristic 
impedance  calculations  as  well. 

The  coefficients  in  Appendix  E  for  these  three  possible 
sets  of  solutions  can  then  be  derived  in  normalized  form, 
where  an  imaginary  quantity  is  written  as   u  =  j  u"  .   It  is 
found  that  the  transform  of  the  dyadic  Green's  function 
component  matrix  [M]  is  imaginary  and  that  the  evaluation  of 
each  single  component  can  be  based  upon  an  evaluation  of  real 
quantities  only  for  the  three  cases  mentioned  above. 

Similar  to  the  derivation  indicated  in  Appendix  G  it  is 
now  possible  to  reduce  these  three  sets  of  equations  to  one 
set,  where  the  following  notation  is  used  for  the  three  types 
of  solutions: 

Case  a)   upper  signs  apply 

Case  b)  middle  signs  apply,  'sinh  Y2D'  should  be  replaced  by 
•sin  Yo15'  and  'cosh  Y?1^'  ky  'cos  YoD'  i  furthermore 
fYoD'  replaces  'Y?0'  • 


Case  c)   lower  signs  apply,  replacements  as  indicated  in 

b)  above.   Furthermore  replace  'sinh  Yih'  by  "sin  yjh'  / 
'cosh  Yin'  by  'cos  y^h'  and  lY1DI  by  'yjD'  . 


The  coefficients  are  then 


+     (k9D)2  -  (CD)2 
F   =  + 

1   -  [(k0D)2-(£D)23  sinh  y-^ 
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(k  D)2   -    UD)2 
F„    = 


[(k0D)2-UD)2]    cosh    y1h2 

+ 

G,    =      a   D    £D    (+    F-     sinh    y,  h0    -    1) 
1  n  _      x  '12 

+ 
G2    =   +    coyD   YtD   F2    sinh    Y1h2 

+ 
G,    =   +    tue0D   y,D   F1    cosh    Yin2 


G.    =    -    a   D    £D    (F2    cosh   Ynn2    "    ^ 


det   G   =    G,G4    +    G2G 


H,    =   co£2D   YoD   G2/det   G 


+ 
H2    =   -    wyD   Y2D   G4/det   G 


+ 
H3   =   -    a3e2D   y2d  G-^/det   G 


H4  =    wyD  Y2D  G3/det  G 


J,  =  -  a  D  £D  (sinh  y2D  +  H,  cosh  y2^)  +  wyD  y2&  H3  sinh  y2D 

+ 
J2  =  +  wyD  y2D  (cosh  y2D  -  H4  sinh  y2D)  -  anD  £D  H2  cosh  y2D 


J.  =  [(k0D)  -UD)  1  (sinh  YoD  +  Hi  cosh  ^?D) 
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J4    =   H2     ^k2D)2    "    (£D)2]    cosh    ^2D 


det   J   =   J, J.    -   J0J_ 
14  2    3 


K 


1   -   {+    anD    ^D  J4   +    [(k0D)2-(CD)2]J2|    sinh   y^  /det   J 


+ 
K2    =   +    coyD   YiD   J/    sinh   y-ihj/det   J 


=   +  {    «nD    CD   J3   +     [  (k0D)2-(^D)2]J1|    sinh    y^/det   J 


K3   = 


K.   =  -   wyD   YiD  J.   sinh   Yi^n/det  J 


L-,    =   -    co£0D   y,D    cosh   Ynh,    +    we2D   YoD  Ki    cosh   Y2D 


-    we2D    y2D    (H-.K      -   H2K3)     sinh   Y2D   +    a   D    CD   K_    sinh    Y2D 


+ 
+    aRD    CD    (H3K1   +   H4K3)     cosh   y2D 


L2    =  aRD    CD   cosh    Yini    +    ooe2D    Y2D  k2    cosh    y2d 

+  + 

-  W£2D  y2d  (H!K2  +  H2K4J  Sinh  Y2D  "  anD  ^D  K4  sinh  Y2D 


+  aRD  CD  (H3K2  -  H4K4)  cosh  Y2D 


L3  =  -  [(k2D)2-(CD)2]  K3  sinh  y2D 


2  2  + 

-     [(k2D)     ~(CDr]      (H3K1    +    H4K3)     COSh    Y2D 
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L4  =  [ (kQD)2-(£D)2]  cosh  y1h1  +  [ (k2D) 2- ( CD) 2 ]K4  sinh  y2D 

+        2      2 

-  [(k2D)z-UD)z]  (H3K2  -  H4K4)  cosh  y2D 


M,    =   -   L2/wyD   y,D   sinh   Yih, 

+  2  2 

M2   =  +    [Lx   -    anD    CD  L2/wyD  y^D]  /  [  (kQD)    -  (  CD)  *  ]    sinh   Y1h]_ 


M~    =   -   L./coyD   y,D   sinh   y,  h.. 

T 

+  2      2 

M4  =  +  [L3  -  anD  CD  L4/wyD  y]_D]  /  [  (kQD)  -(CD)  ]  sinh  y^ 


While  the  above  equations  are  applicable  for  the  dispersion 
characteristic  calculation  of  shielded  slot  lines  and  slot  line 
resonator  lengths,  some  further  manipulations  are  necessary 
to  provide  means  for  a  numerical  evaluation  of  the  time 

average  power  flow  computation  for  slot  lines. 

e  h 

The  dependence  of  the  coefficients  A  through  D   on  the 

assumed  electric  field  distribution  is  given  in  Appendix  E. 

Rewriting  these  equations  in  the  normalized  form  and  combining 

the  two  sets  of  equations  for  the  two  possible  solutions  as 

?  2 

above  where  y2 '  >  0  and  y2   <  0  then  leads  to  the  following 

expressions . 


Ae 

— j   =  0   since  in  the  first  order  approximation  it  is  assumed 

that  B      as  0  . 
z 
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A   _ x 

72"      wyD  y-j^D  sinh  YiD  h-^/D 

Be  .  -  „   Ah 
?  "  +   2  D2 

B*  „  +  „   A* 
?  "  "   4  D2 

CS  _  +  „   B6  -     Bh 

?    -.Hi  7  +  2  ? 

rh        Re       h 

~T     H3  ~7   H4  ~T 
D  D        D 

DS  _  p   Ce 

T  "  1  T 

D        D   _ 

Dh  =  p   Ch 

^2"  "   2  IT 

D        D 

2 

where  the  upper  sign  applies  for   y2   >  0   an<^  the  lower  sign 

2 

for   y2   <  0  .   (Note  that  a  quite  similar  notation  was  used 

for  the  expressions  in  Appendix  H  for  the  open  boundary 

2 

transmission  lines,  written  there  separately  for   y?  >    0 

2 

and   Y?   <  0  while  the  above  form  shows  the  advantage  of  a 

more  composite  form. ) 

The  final  step  for  a  preparation  for  numerical  analysis 
is  then  the  normalization  of  the  expressions  for  the  discrete 
Fourier  transformed  power  density  terms  given  in  Appendix  F. 


232 


For    region    1, 

ftn      Ah    9     f                o                o      sinh    2y1Dh1/D 
Pin   =   wyD    W    ,^,2    |[(anD)2+(T1D)2]    ^ 


^    [(anD)2-(YlD)2]} 


and  for  region  3, 


n    BD    ,De.     /.,       ,,,2^,       n,2.     Sinh    2yjDh    /D 
P3n   -   ue0D  —    (^>    \[(anD)    +(YlD>    '    2r[D 

-  -£    [(anD)2-(YlD)2i} 
+   wyD    »    (2^)    {[(anD)2+(YlD)2]    Slnh22^lDh2    D 

+  ^  [(anD)2-(YlD)2]} 

e   h 

+  a  D  [(3D)2+(knD)2]  ^5-  ^r  sinh  2YlDh9/D   . 
n  u      Dz    Dz 


The  two  expressions  for  region  2  for  y2   >  0   and   y?   <  0 
are  readily  available  with  the  integrands  of  equations  (H-2a) 
and  (H-2b)  in  Appendix  H  due  to  the  same  nomenclature  used 
for  the  solutions  to  the  Helmholtz  equations  for  open  and 
closed  boundary  slot  lines  in  equations  (2-12,  2-13)  and 
(3-17,3-18)  where  the  continuous  variable  otD  is  replaced  by 
the  discrete  variable  a  D. 
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The  necessary  translation  into  a  digital  computer  language 
is  then  routine  and  can  be  easily  seen  in  the  subroutine  "COKBC 
of  the  attached  computer  program  "RESO"  which  computes 
parameters  for  shielded  slot  lines  and  resonators.   A 
consistent  nomenclature  with  the  one  used  above  has  been  used. 
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APPENDIX  K 
Circular  Resonating  Aperture  on  a  Dielectric  Substrate 

This  appendix  presents  an  attempt  to  analyze  a  circular 
slot  on  a  dielectric  substrate  enclosed  with  conducting 
cylindrical  walls  of  a  circular  waveguide  closed  at  both 
ends  by  conductors  across  the  cross-sections  of  this  guide. 
A  schematic  of  this  configuration  is  shown  in  Figure  K/l . 

This  geometry  is  considered  to  be  analyzed  if  depending 
on  the  assumed  field  mode  of  operation  for  a  certain  set  of 
geometric  parameters  the  resonant  frequency  can  be  found. 

The  following  derivation  is  not  meant  to  be  complete  and 
rigorous  but  is  merely  intended  to  outline  possible  problems 
an  investigator  may  encounter.   For  analytical  purposes  the 
conductors  and  the  dielectric  are  assumed  to  be  ideal  and 
the  dimensions  of  the  surrounding  walls  are  such  that  a 
regular  circular  cavity  mode  cannot  exist. 

It  is  a  well  known  fact  that  the  scalar  potential  functions 
satisfy  the  Helmholtz  equation  which  is  in  cylindrical 
coordinates 


f  If  <r  §F>  +  \  ?4  +  aA  +  k2*  ■  °    •  »-l» 

r   30     3z 


If  a  general  formulation  of  the  resonating  problem  via  a 
transformation  is  sought  —  following  the  guidelines  in 
Section  2  and  3  —  the  partial  differential  equation  (K-i) 
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/T\ 


>-  r 


Figure  K/l   Circular  Resonant  Aperture  on  a  Dielectric 
in  a  Closed  Cylinder 
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is  to  be  transformed  into  an  ordinary  differential  equation. 
This  suggests  the  use  of  a  two  dimensional  transform  which 
consists  of  a  finite  Fourier  transform  in  the  variable  0  and 
a  finite  Hankel  transformation  in  the  variable  r. 
Thus  define 


H(n,Sm,z)  = 


7T 

/  e 

-7T 


jn6 


/  r  B  U  r)  h(6,r,z)  dr  d9 
0 


(K-2) 


where  B   is  the  Bessel  function  of  the  first  kind  of  order  n 
n 

It  is  easily  shown  that  this  transform  applied  to 
equation  (K-l)  will  transform  the  differential  operators  of 
the  variables  6  and  r  into  algebraic  operations  so  that  an 

ordinary  differential  equation  in  the  variable  z  remains. 

d  2 

Obviously  the  differential  operator  — «-  is  replaced  by  -n 

so  that  equation  (K-l)  via  (K-2)  becomes 


/  r  Bn (£  r) {I  !-  (r  2£)  -  K   h  +  i-£  +  k2h}  dr  =  0   (K-3) 
0      n  m   r  ^    W         r1  dz2 


where  h  is  now  a  function  of  n,  r  and  z. 

Integration  by  parts  of  the  first  two  terms  will  then 
lead  to 


W>  r  w 


-  h  •  r  k 

0  8r 


^B  (Er) 


+  I        7?    (r  ^"ar^^  h  dr  "  ^  ^   F  Bn( V>  h  dr 


(1_  +  k2)  /   r  B(£  r)  h  dr   =   0 
8z^       0 


(K-4) 
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Bessel's  equation  is  known  as 


3B  U  r)       2      9 
|_  (r  -VL-,  =  <^-?ra2r>  Bn(?mr) 


so  that  equation  (K-4)  can  be  written  as 


2 

1-2-  H(n,£m,z)  =  (Km2   -   k2)  H(n,Cm,z)  (K-5) 

3z 


which  is  the  desired  ordinary  differential  equation. 
However,  this  only  holds  if  the  boundary  terms  in  (K-4) 
vanish.   These  terms  must  now  be  investigated  in  more  detail. 

But  first  a  sidestep  is  to  be  taken  to  formulate  the 
behavior  of  the  field  components  and  their  dependence  on  the 
scalar  potential  functions.   The  fields  are  contained  within 
the  circular  waveguide  cavity  so  that  no  energy  is  lost  due 
to  radiation.   This  fact  makes  the  problem  symmetrical  in  the 
8,  r  and  z  coordinates  because  standing  waves  do  exist  in 
these  three  directions.   But  then  the  fields  may  be  expressed 
as  TE  or  TM  to  any  one  of  these  coordinates  where  a  conven- 
tional choice  suggests  the  use  of  the  z  coordinate,  similar 
to  cylindrical  cavities. 

According  to  Collin  [Ref .  14]  one  may  assume  a  magnetic 
Hertzian  potential   iL  =  a   ijj(6,r,z)   and  an  electric 
Hertzian  potential   S  =  a   <j>(6,r,z)   for  a  full  or  hybird 
mode  analysis. 


238 


The  field  quantities  in  the  three  cylindrical  coordinates 
can  then  be  expanded  by  Maxwell's  equations  to  yield 

„         92<f>         .1  dip 
i  s2<t>  aui 

Ez   =    (~T+  k2)* 

o  Z 
H       -     8^       +    W     1     3* 

Hr  -  SrTz"  +   3°J£   F  STT 
He   "  F  T6Tz"  "   DW£  3r 


(K-6) 


Hz  -  (7^+  k^J 
8z 


Bessel  functions  of  the  first  kind  are  finite  at  the 

origin  and  have  finite  derivatives;  similarly  the  scalar 

potential  functions  for  this  configuration  will  be  finite 

with  finite  derivatives  so  that  the  boundary  terms  in  equation 

(K-4)  vanish  at  the  lower  limit.   At  the  upper  limit  r  =  a 

where  E   and  Efi  must  vanish  (tangential  electric  field 

vanishes  at  the  surface  of  a  conductor) ,  one  of  the  boundary 

conditions  for  the  scalar  potential  by  equation  (K-6)  is 

then   <f>(a,0,z)  =  0.   This  requires  then  that  Bn(£ma)  =  0 

so  that  a  discrete  set  of  transform  variables  £   with 

m 

m  =  1,2,3,...  is  obtained  where   £  •  a   is  the  m-th  zero  of 
the  n-th  order  Bessel  function. 
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While  this  derivation  holds  for  the  transformation  of 
the  potential  function  <$>    it  is  not  obvious  that  the  same 
set  of  transform  variables  will  apply  for  the  transformation 
of  the  function  i|;  for  which  the  boundary  condition  is  J^-  =   0. 

A  further  complication  in  this  analysis  can  be  foreseen 
as  follows.   For  the  matching  of  the  boundary  conditions  at 
the  two  interfaces  at  z  =  0  and  z  =  D  the  field  components 
in  the  r  and  the  0  direction  must  be  transformed  where  it  is 
observed  that  the  field  components  in  the  6  direction  consists 
of  partial  derivatives  of  the  scalar  potential  functions  with 
respect  to  the  variable  r.   The  finite  Hankel  transform  of 
a  derivative  of  a  function  is  given  by  Sneddon  [Ref.  25]. 
For  instance,  if  n  ^  0  and  if 


W    =   £     r  Bn(^mr>    9(r>    dr 


then 


/    r  W>  ai  dr  "  TST  Vn-l'V   "  TT  VWV    • 


It  is  seen  that  a  derivative  in  the  space  domain  is  thus 
transformed  into  a  algebraic  operation  in  the  transform 
domain  but  the  order  of  the  transforming  Bessel  function  has 
to  be  changed  which  further  complicates  this  analysis. 

As  one  may  observe  at  this  point,  the  relatively  straight- 
forward transform  method  in  a  rectangular  coordinate  system 
is  anything  but  simple  in  a  problem  of  cylindrical  geometry. 
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Finally,  it  should  be  mentioned  that  if  the  Hertzian 
potentials  are  assumed  to  be  in  the  other  directions  the  set 
of  equations  (K-6)  is  changed  in  such  a  way  that  the  boundary 
conditions  are  even  more  difficult  to  handle. 

This  problem  was  abandoned  here,  since  it  was  felt  that 
if  this  method  of  analysis  is  applicable  at  all  further  study 
would  more  appropriately  be  the  subject  of  new  thesis  work. 
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APPENDIX  L 
Design  Curves 

This  appendix  includes  some  design  curves  for  single 
slot  and  coplanar  strip  transmission  lines.   The  appropriate 
parameters  can  be  found  in  the  figure  captions. 
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Figure  L/l   Dispersion  Characteristic  and  Characteristic 
Impedance  of  Single  Slot  Line  with  er  =  12 
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Figure  L/2   Dispersion  Characteristic  and  Characteristic 
Impedance  of  Single  Slot  Line  with  e   =16 
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Figure  L/3   Dispersion  Characteristic  and  Characteristic 
Impedance  of  Single  Slot  Line  with  e   =20 
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Figure   L/4a      Dispersion   Characteristic    of   Coplanar 
Strip   Line   with  W/D   =   1.5    and    er    ~-    2.5 
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Figure  L/4b   Characteristic  Impedance  of  Coplanar 

Strip  Lire  with  W/D  =  1.5  and  e   =  2.5 
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Figure  L/5a   Dispersion  Characteristic  of  Coplanar 
Strip  Line  with  W/D  =1.5  and  er  =  9.0 
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Figure  L/5b  Characteristic  Impedance  of  Coplanar  Strip  Line 
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Figure  L/6a   Dispersion  Characteristic  for  Coplanar 
Strip  Line  with  W/D  =  1.5  and  e   =  20 
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Figure  L/6b   Characteristic  Impedance  of  Coplanar 
Strip  Line  with  W/D  =1.5  and  er  =  20 
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APPENDIX  M 
Computer  Program  * COPLAN' 
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APPENDIX   N 
Computer   Program    'PESO' 
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